Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 550,089 Δ1 bp 6.4% coding (761/1194 nt) tufB → elongation factor Tu‑B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947550,0890C. 6.4% 164.1 / 6.4 62coding (761/1194 nt)tufBelongation factor Tu‑B
Reads supporting (aligned to +/- strand):  ref base C (24/33);  new base . (2/2);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GAGCCAGTTCGTGCCATCGACCAGCCGTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCGAGAACTGCGGCGTTCTGCTGCGTGGTACCAAGCGTGACGAC  >  NC_002947/549941‑550195
                                                                                                                                                     |                                                                                                          
gAGCCAGTTCGTGCCATCGACCAGCCGTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGAcacc                                                                                                           <  2:580388/150‑1 (MQ=35)
 aGCCAGTTCGTGCCATCGACCAGCCGTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTCCCGTGTTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGAcacca                                                                                                          >  1:758599/1‑150 (MQ=255)
     aGTTCGTGCCATCGACCAGCCGTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGAcaccaccac                                                                                                      <  2:940919/150‑1 (MQ=35)
          gTGCCATCGACCAGCCGTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGTCCGTA‑CGATCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGAcaccaccaccacca                                                                                                 >  1:1334001/1‑150 (MQ=255)
            gCCATCGACCAGCCGTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGAcaccaccaccacca                                                                                                 <  1:1923462/148‑1 (MQ=35)
            gCCATCGACCAGCCGTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGAcaccaccaccacca                                                                                                 >  2:1923462/1‑148 (MQ=35)
            gCCATCGACCAGCCGTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACaccaccaccacctcc                                                                                               <  2:2064324/150‑1 (MQ=21)
             ccATCGACCAGCCGTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCt                                                                                              >  1:1722092/1‑150 (MQ=34)
             ccATCGACCAGCCGTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCt                                                                                              >  2:384522/1‑150 (MQ=34)
             ccATCGACCAGCCGTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTCAGCGCCCCTTCCACCACCACCt                                                                                              >  2:1637477/1‑150 (MQ=11)
                tCGACCAGCCGTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCa                                                                                           >  1:2085527/1‑150 (MQ=32)
                    ccAGCCGTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGA‑‑‑CACCACCACCACCTGCACCGgtgt                                                                                    >  1:1925780/1‑150 (MQ=255)
                      aGCCGTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGgtg                                                                                     >  1:2316203/1‑150 (MQ=32)
                       gCCGTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGgtgt                                                                                    >  2:1577292/1‑150 (MQ=32)
                        ccGTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGtt                                                                                   <  2:1876488/150‑1 (MQ=32)
                        ccGTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGACCCTGCGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGtt                                                                                   <  1:1636750/150‑1 (MQ=255)
                          gTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAACTCGTTGGTCTGCGTGACACCACCACTTCCACCTGCACCGGTGTTga                                                                                 <  2:289484/150‑1 (MQ=12)
                          gTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGA‑‑‑CACCACCACCACCTGCACCGGTGTTGAGAt                                                                              >  2:1335601/1‑150 (MQ=255)
                              cTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATg                                                                             >  1:2287717/1‑150 (MQ=32)
                                 ttGCCGATCGAAGACGTGTTCTCGCTCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTc                                                                          <  2:1045193/149‑1 (MQ=17)
                                 aTGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCTGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTc                                                                          <  2:203714/150‑1 (MQ=17)
                                 aTGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTc                                                                          <  1:466556/150‑1 (MQ=32)
                                  tGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTcc                                                                         >  2:1082853/1‑150 (MQ=32)
                                  tGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTcc                                                                         <  2:758599/150‑1 (MQ=32)
                                      gATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCaa                                                                     >  1:2214930/1‑150 (MQ=32)
                                      gATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATAGCGGCGCGTCGGGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCaa                                                                     >  2:1244672/1‑150 (MQ=11)
                                       aTCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAg                                                                    >  2:2212978/1‑150 (MQ=32)
                                       aTCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAg                                                                    <  2:2362177/150‑1 (MQ=32)
                                       aTCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAg                                                                    >  2:991743/1‑150 (MQ=32)
                                         cgatctcGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCCCATGCACCGGTGTTGAGATGTTCCGCAAgct                                                                  <  1:1881504/144‑1 (MQ=255)
                                            aGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGA‑‑‑CACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGa                                                            <  1:1335601/150‑1 (MQ=255)
                                            aGACGTGTTCTCGATCTCGGGTCGTGCTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAgctgct                                                               <  2:2400048/150‑1 (MQ=17)
                                               cGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGa                                                            <  1:1696941/150‑1 (MQ=32)
                                               cGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGa                                                            <  2:2046857/150‑1 (MQ=32)
                                                gtgtTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGATCTGACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAgctgagatcg                                                           >  2:1881504/1‑144 (MQ=11)
                                                gtgtTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGAc                                                           <  1:1577292/150‑1 (MQ=32)
                                                 tgttgtCGATCTCGGGTCTTGGTAGCGTTTTTACCGGT‑CGTATCGAGCGTGCCATCGTTCTCGTTCATGATCCGCTGGAAATCGTTGCTCTGCGT‑‑CA‑CACCACCCCCGCCTGCACCGGTGTTGACGTGTTCCGCAAGCTGCTGGACGAAg                                                       <  2:710685/150‑1 (MQ=11)
                                                 tgtTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGAt                                                                              >  2:1980465/1‑130 (MQ=33)
                                                 tgtTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCAGCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGAt                                                                              <  1:1980465/130‑1 (MQ=18)
                                                       cGATCTCGGGTCGTGGTATGCTGGTTACCGGG‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTc                                                    <  1:722220/150‑1 (MQ=11)
                                                       cGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTc                                                    >  1:2204679/1‑150 (MQ=32)
                                                       cGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTc                                                    <  1:384522/150‑1 (MQ=32)
                                                        gATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGa                                                         >  1:2175090/1‑144 (MQ=32)
                                                        gATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGa                                                         <  2:2175090/144‑1 (MQ=32)
                                                           ctcGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGCCACCACCACCACGACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGc                                                <  2:1330725/150‑1 (MQ=17)
                                                              gggTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTgg                                             <  2:1334001/150‑1 (MQ=32)
                                                                 tcatGGTACCGTTGCTACCGGC‑CGTATCGAGCGGGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCga                                          <  2:1716622/147‑1 (MQ=255)
                                                                 tCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCga                                          <  2:281897/150‑1 (MQ=32)
                                                                    tGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCaa                                                                     <  2:1810240/120‑1 (MQ=33)
                                                                    tGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCaa                                                                     >  1:1810240/1‑120 (MQ=33)
                                                                         ccGTTGTGACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCGAGAACTgcg                                  <  1:711857/150‑1 (MQ=17)
                                                                          cGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCga                                          >  1:2419134/1‑141 (MQ=32)
                                                                          cGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCga                                          <  2:2419134/141‑1 (MQ=32)
                                                                              gttACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAgc                                                                   <  1:462828/112‑1 (MQ=33)
                                                                              gttACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAgc                                                                   >  2:462828/1‑112 (MQ=33)
                                                                               ttACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGCGCCTCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCGAGAACTGCGGCGTTc                            <  1:23551/150‑1 (MQ=12)
                                                                                   cGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCGAGAACTGCGGCGTTctgct                        >  1:2412275/1‑150 (MQ=32)
                                                                                      t‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAGGCTGCTGGACGAAGGTCGTGCTGGCGAGAACTGCGGCGTTCTGCTGCg                     <  2:2316203/150‑1 (MQ=17)
                                                                                          tATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCGAGAACTGCGGCGTTCTGCTGCGTgg                  >  2:174140/1‑150 (MQ=32)
                                                                                          tATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCGAGAACTGCGGCGTTCTGCTGCGTgg                  <  2:2204679/150‑1 (MQ=32)
                                                                                           aTCGAGGGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCGAGAACTGCGGCGTTCTGCTGCGTGGt                 <  1:181468/150‑1 (MQ=17)
                                                                                            tCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCGAGAACTGCGGCGTTCTGCTGCGTGGTa                <  2:1584376/150‑1 (MQ=32)
                                                                                               aGCGTGGCGTCGTTCGCGTTCAGGCTCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCGCCGGTGTGGAGATGTTCGGAAAACTGCTGGACGAAGGGCGTGCTGGCGAGAACTGCGGCGCTCTGCTGCCTGGTACCa             >  2:267514/1‑150 (MQ=11)
                                                                                                gCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCGAGAACTGCGGCGTTCTGCTGCGTGGTACCaa            >  2:480427/1‑150 (MQ=32)
                                                                                                  gTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCGAGAACTGCGGCGTTCTGCTGCGTGGTACCAAGc          <  1:88153/150‑1 (MQ=32)
                                                                                                     gCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACAGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCGAGAACTGCGGCTTTCTGCTGCGTGGTACCAAGCGTg       >  2:771472/1‑150 (MQ=17)
                                                                                                      cATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCGAGAACTGCGGCGTTCTGCTGCGTGGTACCAAGCGTga      >  1:1873363/1‑150 (MQ=32)
                                                                                                      cATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCGAGAACTGCGGCGTTCTGCTGCGTGGTACCAAGCGTga      <  2:1873363/150‑1 (MQ=32)
                                                                                                          gTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCGAGAACTGCGGCGTTCTGCTGCGTGGTACCAAGCGTgacgac  <  1:1082853/150‑1 (MQ=32)
                                                                                                                                                     |                                                                                                          
GAGCCAGTTCGTGCCATCGACCAGCCGTTCCTGATGCCGATCGAAGACGTGTTCTCGATCTCGGGTCGTGGTACCGTTGTTACCGGT‑CGTATCGAGCGTGGCATCGTTCGCGTTCAGGATCCGCTGGAAATCGTTGGTCTGCGTGACACCACCACCACCACCTGCACCGGTGTTGAGATGTTCCGCAAGCTGCTGGACGAAGGTCGTGCTGGCGAGAACTGCGGCGTTCTGCTGCGTGGTACCAAGCGTGACGAC  >  NC_002947/549941‑550195

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