Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,936,634 A→T 7.9% K119N (AAA→AAT rbsB → D‑ribose ABC transporter periplasmic binding protein; ribose chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,936,6340AT 7.9% 186.1 / 4.8 62K119N (AAA→AATrbsBD‑ribose ABC transporter periplasmic binding protein; ribose chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base T (3/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGA  >  NC_000913/3936494‑3936780
                                                                                                                                            |                                                                                                                                                  
gCGACCGTGCAGGACTTACCCGTTCGCGGCACAAAAATTCTGCTGATTTACCCGCCCGCCTCCGACGCAGTGGGTAATGGTGTGAAGATGGCTAACCAGGCTATCATACCGGGTATCTCTCGTGGGCTCCAGGCAAGGAATGGGGAAg                                                                                                                                             >  6:118513/1‑148 (MQ=255)
gCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGa                                                                                                                                               <  3:56646/146‑1 (MQ=255)
 cGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTTGGTAATGCTGTGAAGATTGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAgt                                                                                                                                            <  5:118513/148‑1 (MQ=255)
   aaCGTGCAGGACTTAACCGATCGCGGAACAAAAATACTGCTGATTAACCCGACCGACTCCGACGCAGAGGGAAAAGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCGAGGCAACGAAAGGTGAAgtgg                                                                                                                                          >  5:335375/1‑148 (MQ=255)
           ggACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCcac                                                                                                                                  <  2:110810/148‑1 (MQ=255)
                     gTTCGCGGCACAAAAATTCTGCTGATTAACCCGCCCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCt                                                                                                                         >  4:274493/1‑147 (MQ=255)
                       tCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGa                                                                                                                       >  4:350892/1‑147 (MQ=255)
                         gcgGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATaa                                                                                                                    <  2:302766/148‑1 (MQ=255)
                         gcgGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATa                                                                                                                     <  3:33036/147‑1 (MQ=255)
                                aaaaaTTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGa                                                                                                                                               <  1:142117/114‑1 (MQ=255)
                                aaaaaTTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGa                                                                                                                                               >  2:142117/1‑114 (MQ=255)
                                 aaaaTTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTgg                                                                                                            >  1:168888/1‑148 (MQ=255)
                                 aaaaTTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTgg                                                                                                            >  6:61800/1‑148 (MQ=255)
                                               tAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGTACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCt                                                                                              >  3:290326/1‑148 (MQ=255)
                                                 aCCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTg                                                                                             <  5:61800/147‑1 (MQ=255)
                                                 aCCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGGGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTg                                                                                             <  3:350892/147‑1 (MQ=255)
                                                  cccGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGt                                                                                           <  3:274493/148‑1 (MQ=255)
                                                   ccgaccgaCTCCTACGCAGTGGGTAATGCTGTTAAGATGGCTAACCAGGGGAACATCCCGGTATTCTCTCTTGACCGCCAGGCTACGAAAGGTTAAGTGGTTAGCCACATTGCTTCTGATAACGTACTGGGCGGCATAATCGCTGGTg                                                                                          <  6:335375/148‑1 (MQ=255)
                                                   ccgaccgaCTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTg                                                                                                                        >  6:65898/1‑118 (MQ=255)
                                                   ccgaccgaCTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTg                                                                                                                        <  5:65898/118‑1 (MQ=255)
                                                                       gggTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAg                                                                                                                                      <  2:112099/84‑1 (MQ=255)
                                                                       gggTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAg                                                                                                                                      >  1:112099/1‑84 (MQ=255)
                                                                       gggTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCg                                                                      >  7:21067/1‑148 (MQ=255)
                                                                       gggTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCg                                                                      >  5:250624/1‑148 (MQ=255)
                                                                       gggTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCg                                                                      >  4:194690/1‑148 (MQ=255)
                                                                          tAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGt                                                                   <  8:21067/148‑1 (MQ=255)
                                                                             tGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGaa                                                                <  6:250624/148‑1 (MQ=255)
                                                                               cTGTGAAGATGGCTAACCAGGCGTACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAg                                                               <  4:290326/147‑1 (MQ=255)
                                                                                tgtgAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGt                                                             >  4:337388/1‑148 (MQ=255)
                                                                                        tGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGt                                                     >  1:309979/1‑148 (MQ=255)
                                                                                          gCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTa                                                   >  3:69819/1‑148 (MQ=255)
                                                                                           cTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGa                                                                                         <  5:62875/109‑1 (MQ=255)
                                                                                           cTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGa                                                                                         >  6:62875/1‑109 (MQ=255)
                                                                                              aCCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAATGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCg                                                <  6:385063/147‑1 (MQ=255)
                                                                                                 aGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCt                                            >  1:155550/1‑148 (MQ=255)
                                                                                                           cccGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCAtt                                  >  6:205746/1‑148 (MQ=255)
                                                                                                                   tCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGc                                     >  8:233937/1‑137 (MQ=255)
                                                                                                                   tCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGc                                     <  7:233937/137‑1 (MQ=255)
                                                                                                                    cACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTAc                          >  2:262720/1‑147 (MQ=255)
                                                                                                                        ctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGa                                                                 >  6:179756/1‑104 (MQ=255)
                                                                                                                        ctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGa                                                                 <  5:179756/104‑1 (MQ=255)
                                                                                                                        ctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGa                                               >  1:346532/1‑122 (MQ=255)
                                                                                                                        ctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGa                                               <  2:346532/122‑1 (MQ=255)
                                                                                                                        ctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGt                            >  8:109498/1‑141 (MQ=255)
                                                                                                                        ctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGt                            <  7:109498/141‑1 (MQ=255)
                                                                                                                         ttGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATcc                      <  3:337388/146‑1 (MQ=255)
                                                                                                                            aCCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGc                                     <  8:248659/128‑1 (MQ=255)
                                                                                                                            aCCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGc                                     >  7:248659/1‑128 (MQ=255)
                                                                                                                               gCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCg              <  5:205746/148‑1 (MQ=255)
                                                                                                                                  aGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGa           <  5:206106/148‑1 (MQ=255)
                                                                                                                                  aGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGa           <  4:69819/148‑1 (MQ=255)
                                                                                                                                   ggCAACGAATGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAg                                                               >  8:279740/1‑95 (MQ=255)
                                                                                                                                   ggCAACGAATGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAg                                                               <  7:279740/95‑1 (MQ=255)
                                                                                                                                   ggCAACGAATGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGAGAAGAAAGCGGGTGAAg                                                               >  8:279739/1‑95 (MQ=255)
                                                                                                                                   ggCAACGAATGGTGAAGTGGTGAGACACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAg                                                               <  7:279739/95‑1 (MQ=255)
                                                                                                                                   ggCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCAATGGTCGTACATCCGCAGCCCGTGaa          >  6:206106/1‑148 (MQ=255)
                                                                                                                                    gCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGAAc         <  2:168888/148‑1 (MQ=255)
                                                                                                                                    gCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGAAc         <  3:194690/148‑1 (MQ=255)
                                                                                                                                        cGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGAACGTgg     >  4:150578/1‑148 (MQ=255)
                                                                                                                                         gAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGaa          <  7:102635/142‑1 (MQ=255)
                                                                                                                                         gAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGaa          >  8:102635/1‑142 (MQ=255)
                                                                                                                                         gAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGAACGTGGc    <  2:309979/148‑1 (MQ=255)
                                                                                                                                           aaGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGa  <  3:150578/148‑1 (MQ=255)
                                                                                                                                           aaGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGa  <  2:155550/148‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
GCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGA  >  NC_000913/3936494‑3936780

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: