Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 196,795 G→C 8.2% T767T (ACG→ACC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947196,7950GC 8.2% 194.2 / 10.1 77T767T (ACG→ACCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (34/34);  new base C (3/3);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

CGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCT  >  NC_002947/196657‑196935
                                                                                                                                          |                                                                                                                                            
cGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGAc                                                                                                                                   <  1:1051087/150‑1 (MQ=40)
cGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGAc                                                                                                                                   <  2:1797999/150‑1 (MQ=40)
  ttGAAGCCGGCAAGATCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGGCCGCGTCACCATCGAGAACGCAACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCCGCTCAGACCa                                                                                                                                 <  2:1629179/150‑1 (MQ=11)
         cGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCaa                                                                                                                          <  2:1290333/150‑1 (MQ=38)
          ggCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGTCTGCAACTTCTAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAAc                                                                                                                         <  1:1454095/150‑1 (MQ=255)
                  ccAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGGCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGCCGGCCACTTCGAAGAGCTCACTCCGCACCCGACACCGGCTCAGCGCACAATCAACAACTcctg                                                                                                                 >  2:292993/1‑147 (MQ=11)
                  ccAGGGCAGCGGCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGCGCACGGGCCTGGGGGCAAGTGCGAGCAGTTGACCACGAACCCGACGCCGCCTCAGACCACGATCAACGAACCGGt                                                                                                                 >  2:31483/1‑150 (MQ=11)
                            gTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCTTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCaa                                                                                                                          >  1:1918140/1‑131 (MQ=37)
                            gTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCaa                                                                                                                          <  2:1918140/131‑1 (MQ=35)
                            gTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGccacc                                                                                                       >  1:1473646/1‑150 (MQ=35)
                                   tCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGa                                                                                                <  1:1803463/150‑1 (MQ=35)
                                     cAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACcgcgctt                                                                                              <  1:1835096/150‑4 (MQ=25)
                                      aGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACt                                                                                                                     >  1:470606/1‑126 (MQ=37)
                                      aGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACt                                                                                                                     <  2:470606/126‑1 (MQ=35)
                                          cccGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGAc                                                                                         >  1:1518182/1‑150 (MQ=35)
                                                 aaTGACGTGTACAACAAAGCTTCGACCGCAAGCTTCACGGTGTAGAAGGTCGCTGTCTGCAACTTCGAACAGCTGACCCCGATCCCGACGCCGGCTCAGACCACGATCAACGACGCGGTCGATGCCACCACCGCGACCCTGACGGCGAGc                                                                                  <  2:255606/150‑1 (MQ=11)
                                                       gTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCg                                                                            <  1:1378232/150‑1 (MQ=32)
                                                       gTCTACAACAACGGCTCGACCGTCAGCGTCACCAGCGAGAACGCCACTGGCGGCAACTTCGCCCAGCTGCCCCCGAACCCGACGCCGGCACAGACGACGATCAACGACTAGGTCGAGGCCGCCACCGCGACCCTCACGGCGACCccgccg                                                                            >  2:570511/1‑150 (MQ=11)
                                                        tCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCgg                                                                           >  1:2002246/1‑150 (MQ=32)
                                                         cTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGt                                                                          <  1:1270428/150‑1 (MQ=32)
                                                           acaacaACGGTTCGACCGTCAGCGATCGTACGGAGAACGCCACTGGCGGCAACTTCGAACATCTTACCCCGATCCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGAAGGCGAGCCCGTCGGTCa                                                                        >  1:1228705/1‑150 (MQ=11)
                                                            caacaaCGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGAt                                                                                                             <  1:658746/112‑1 (MQ=33)
                                                            caacaaCGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGAt                                                                                                             >  2:658746/1‑112 (MQ=33)
                                                             aacaacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCAcc                                                                      >  1:1367210/1‑150 (MQ=32)
                                                             aacaacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCAcc                                                                      >  1:2025896/1‑150 (MQ=32)
                                                                aacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAAc                                                                                                                         <  2:1044683/96‑1 (MQ=34)
                                                                aacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAAc                                                                                                                         >  1:1044683/1‑96 (MQ=34)
                                                                aacGGTTCGACCGTCAGCGTCACCATCAAGAACGCCACTGGCGGCAACTTCGAACAGTTGACCCCGCACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGaa                                                                   <  2:2227848/150‑1 (MQ=12)
                                                                 acGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAg                                                                  <  2:1367210/150‑1 (MQ=32)
                                                                      tCGACCGTCAGCGTCACCATCGAGCACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGTTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcctg                                                             <  1:2200874/150‑4 (MQ=12)
                                                                         aCCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGgcggcg                                                            >  2:1005497/1‑148 (MQ=32)
                                                                         aCCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGgcggcg                                                            <  1:1005497/148‑1 (MQ=32)
                                                                         aCCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTg                                                          <  2:1473646/150‑1 (MQ=32)
                                                                             tCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGt                                                                                                                 >  2:1267433/1‑91 (MQ=34)
                                                                             tCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGt                                                                                                                 <  1:1267433/91‑1 (MQ=34)
                                                                                 cGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTa                                                  <  1:1042757/150‑1 (MQ=32)
                                                                                    cACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTacac                                               >  1:1607987/1‑150 (MQ=32)
                                                                                      ccATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGACGGCTAAGACGACGATCAACGACTAGGTCGATGCCACCAACGCGACACTGAAGGCGAGCGGGTTGGTCAACGAAGGAGGGGTGGTGAGCTacaagg                                             >  2:1116606/1‑147 (MQ=11)
                                                                                       cATCGAGAACGCCACTGGCGTCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGt                                            <  1:148871/150‑1 (MQ=17)
                                                                                       cATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGt                                            >  1:1865265/1‑150 (MQ=32)
                                                                                        aTCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCaccacc                                                                                                    <  1:1348520/93‑1 (MQ=34)
                                                                                        aTCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCaccacc                                                                                                    >  2:1348520/1‑93 (MQ=34)
                                                                                        aTCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTg                                           >  1:352898/1‑150 (MQ=32)
                                                                                         tCGAGAGCGCCACTGGCGGCAACTTCGAACAGCGGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGACCACCTACACCgggc                                          >  2:14150/1‑147 (MQ=255)
                                                                                         tCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAc                                                     >  1:1146256/1‑139 (MQ=33)
                                                                                         tCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAc                                                     <  2:1146256/139‑1 (MQ=32)
                                                                                              aaCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTg                                     >  1:493680/1‑150 (MQ=35)
                                                                                               aCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGa                                    <  1:1963186/150‑1 (MQ=35)
                                                                                               aCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCATGTCATCGACCCCGACGGAGTCCTATGGCTGCGTGATCACCTACACCGTGACCCTGa                                    <  1:292993/150‑1 (MQ=11)
                                                                                                  ccACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCa                                 >  1:986097/1‑150 (MQ=35)
                                                                                                    aCTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAAt                               <  2:986097/150‑1 (MQ=35)
                                                                                                        gcggcAACTTCGAACAGCTGCCCCCGACCCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTg                           <  1:1076937/150‑1 (MQ=17)
                                                                                                        gcggcAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGCCAGCAGCCCGGGCACCTACACCGTGACCCTGAGCAATcccg                           >  2:776023/1‑148 (MQ=11)
                                                                                                          ggcAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGcc                         <  2:493680/150‑1 (MQ=32)
                                                                                                           gcAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGccc                        >  1:1508080/1‑150 (MQ=32)
                                                                                                              aCTTCGAACAGCTGACCCCGAACCCCCCGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGa                     <  1:31483/150‑1 (MQ=17)
                                                                                                                   gAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCt                            <  1:1971262/138‑1 (MQ=33)
                                                                                                                   gAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCt                            >  2:1971262/1‑138 (MQ=33)
                                                                                                                    aaCAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCgg               >  2:1974296/1‑150 (MQ=32)
                                                                                                                         cTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACa          <  2:2002246/150‑1 (MQ=32)
                                                                                                                          tGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCATGCCCAGACGCCGATGACAg         >  2:1234441/1‑150 (MQ=17)
                                                                                                                           gACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGt        <  2:2033577/150‑1 (MQ=9)
                                                                                                                            aCCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTg       >  1:2095460/1‑150 (MQ=17)
                                                                                                                              cccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATc                              <  1:706055/125‑1 (MQ=33)
                                                                                                                              cccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATc                              >  2:706055/1‑125 (MQ=33)
                                                                                                                              cccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAc     >  1:835391/1‑150 (MQ=32)
                                                                                                                              cccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCACGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAc     >  2:8997/1‑150 (MQ=17)
                                                                                                                              cccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGGCGATGCCACCACCGCGACCCTGACGGCGAGCCCGGCGGTCACCGAAGGCGGCGTGATCACCTCAACAGTGACCCGGAGCAATCCTGCCCAGACGCCGGTGACAGTGAc     >  1:1906224/1‑150 (MQ=11)
                                                                                                                              cccGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAc     <  1:1878798/150‑1 (MQ=12)
                                                                                                                              cccGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCAAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAc     >  1:1837792/1‑150 (MQ=12)
                                                                                                                               cggatcccGACAACGACTCAAACCTCGATCAACTACTCGGTCGATACCACCACCGCGTGCCTGAGGGGGAGCCCGTCGGTGACCGAGGGCGGCGTGATCAGCTACACCGTGACCCTGGGCAAGCCTGCCGAGACGCCGGTGACAGTGAcc    >  2:1840950/6‑150 (MQ=11)
                                                                                                                               ccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAcc    >  2:1387856/1‑150 (MQ=32)
                                                                                                                               ccGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAcc    <  2:1359504/150‑1 (MQ=12)
                                                                                                                               ccGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAcc    >  2:1247281/1‑150 (MQ=9)
                                                                                                                                cGAACCCGACGCCGACTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGc                                  <  2:419024/119‑1 (MQ=17)
                                                                                                                                cGAACCCGACGCCGACTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGc                                  >  1:419024/1‑119 (MQ=17)
                                                                                                                                 aaccccGACTCCGCCTCAGACCACGATCAACGACTCGTTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCt  <  1:1116606/147‑1 (MQ=2)
                                                                                                                                          |                                                                                                                                            
CGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCT  >  NC_002947/196657‑196935

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