Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 197,032 C→T 9.9% T846T (ACC→ACT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947197,0320CT 9.9% 914.8 / 64.4 359T846T (ACC→ACTPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (170/152);  new base T (23/13);  total (194/165)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 359 total aligned reads displayed.

GCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGT  >  NC_002947/196910‑197178
                                                                                                                          |                                                                                                                                                  
        gCCGGAGTCAGGGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGGCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTaa                                                                                                                 >  2:693828/1‑150 (MQ=255)
          cGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAAtt                                                                                                               >  1:1062972/1‑150 (MQ=37)
          cGGTGACAGTGACCATGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAAtt                                                                                                               <  1:2073898/150‑1 (MQ=37)
            gTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTc                                                                                                             <  1:1377052/150‑1 (MQ=38)
              gACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                                                                                           <  2:1996125/150‑1 (MQ=35)
                cAGTGACCCTGTCCAACGGCCAACCCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGc                                                                                                         >  1:1252058/1‑150 (MQ=21)
                       ccTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAc                                                                                                  <  2:2172778/150‑1 (MQ=17)
                       ccTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGACCACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAc                                                                                                  >  2:73371/1‑150 (MQ=17)
                       ccTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGt                                                                                                                   >  2:356053/1‑133 (MQ=17)
                        cTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAcc                                                                                                 <  2:1082549/150‑1 (MQ=32)
                        cTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAcc                                                                                                 >  2:114937/1‑150 (MQ=32)
                        cTGTCCAACGGCCAAACCATCACCACCGCATCCGGCAGGACACGGGGCCGCGCCCATTTCCAGATCTCGGCGAGTGCCGTCTCCATCAACAGTTAAACCGCCAGCCGCACCATCGAGAACCCCACAAGCGGCAACTTCGAGCAACTGAcc                                                                                                 >  2:180080/1‑150 (MQ=2)
                              aaCGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACTATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAAt                                                                                           >  1:2185756/1‑150 (MQ=18)
                              aaCGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTaa                                                                                                                 <  2:1500054/128‑1 (MQ=32)
                              aaCGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGa                                                                                                   <  1:1051271/142‑1 (MQ=32)
                               aCGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCg                                                                                                            >  2:667603/1‑132 (MQ=32)
                               aCGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATc                                                                                          >  1:691865/1‑150 (MQ=34)
                                 ggCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                                                                                           >  2:2420951/1‑131 (MQ=33)
                                 ggCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                                                                                           <  1:2420951/131‑1 (MQ=33)
                                 ggCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCg                                                                                        >  1:991379/1‑150 (MQ=34)
                                  gCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGa                                                                                       >  1:1681331/1‑150 (MQ=34)
                                  gCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGa                                                                                       <  1:1559640/150‑1 (MQ=35)
                                   ccAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGAc                                                                                      <  2:531978/150‑1 (MQ=35)
                                     aaaCCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGAccc                                                                                    >  1:1497809/1‑150 (MQ=38)
                                       aCCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                                                                                           >  1:397621/1‑125 (MQ=33)
                                       aCCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCg                                                                                  >  1:672472/1‑150 (MQ=39)
                                       aCCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCg                                                                                  >  2:2245313/1‑150 (MQ=39)
                                         cATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCg                                                                                                                     <  1:2258468/113‑1 (MQ=33)
                                         cATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGc                                                                                <  2:792934/150‑1 (MQ=39)
                                          aTTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGTCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCt                                                                               >  2:2363339/1‑150 (MQ=37)
                                          aTTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCt                                                                               >  1:910842/1‑150 (MQ=38)
                                             aCCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGCCCTCGAATTCGACCCTGGCTTAg                                                                            <  2:1971392/150‑1 (MQ=11)
                                                 ttGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGa                                                                           >  2:1917279/1‑147 (MQ=38)
                                                 ttGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGa                                                                           <  1:1917279/147‑1 (MQ=38)
                                                 ttGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCa                                                                        >  1:814676/1‑150 (MQ=38)
                                                 ttGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCa                                                                        <  1:1534226/150‑1 (MQ=38)
                                                    aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGAcgacga                                                                     >  1:1291849/1‑150 (MQ=38)
                                                    aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGa                                                                     <  2:455314/150‑1 (MQ=38)
                                                         ggCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAAc                                                                >  1:877387/1‑150 (MQ=38)
                                                          gCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGAt                                                                    >  1:1220534/1‑145 (MQ=39)
                                                             aGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACt                                                            <  2:38705/150‑1 (MQ=38)
                                                              gACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGACACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAca                                                                                                 >  2:1378789/1‑111 (MQ=18)
                                                                cccAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATc                                                                   >  1:1757490/1‑140 (MQ=37)
                                                                  cAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCTACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTc                                                       <  2:1248453/150‑1 (MQ=25)
                                                                     ggCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCGTCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAt                                                    <  2:2211255/150‑1 (MQ=37)
                                                                     ggCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAt                                                    <  1:486919/150‑1 (MQ=37)
                                                                     ggCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAt                                                    <  1:1514268/150‑1 (MQ=37)
                                                                     ggCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAt                                                    <  2:576195/150‑1 (MQ=35)
                                                                      atctcGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGGCCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGaca                                                   <  1:1138273/146‑3 (MQ=255)
                                                                       cAGCGTCGATTTCCAGCCACCAGCGAATCACCTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGCAAACGCCCCATGCTGTAATTTCCAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGc                                                  <  1:1306491/150‑1 (MQ=11)
                                                                        aGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGcc                                                 <  2:1864305/150‑1 (MQ=39)
                                                                          cGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGccac                                               >  1:1431916/1‑150 (MQ=40)
                                                                           gTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGAc                                                                          <  2:1767880/122‑1 (MQ=37)
                                                                            tCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCaa                                                                                                       >  1:100434/1‑92 (MQ=34)
                                                                             cGATTTCCAGACCCCGGCGAATGACGTCTACAACAGAAGTGCGACCGTCAGGGACAGCATGGAGACGGCCGCATGCGGTAAGGTCGAGCACCTGACATCGAATACGAACCCGGCTCAGACCACGATCAACGACTCGGTCGATGCcaccac                                            >  2:2061898/1‑150 (MQ=11)
                                                                              taattCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCAGAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCGCAGACCACGAGCACCGACTCGGTCGATGCCaccacc                                           >  2:1174853/4‑150 (MQ=11)
                                                                              gATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGa                                                                           >  2:838094/1‑118 (MQ=37)
                                                                                  tCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGa                                       >  1:2117650/1‑150 (MQ=255)
                                                                                   ccAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAAttt                                                                                                              <  2:1441567/78‑1 (MQ=35)
                                                                                   ccAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGAc                                      >  1:352041/1‑150 (MQ=40)
                                                                                    cAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGAcc                                     >  1:439801/1‑150 (MQ=255)
                                                                                     aGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGAccc                                    <  1:725539/150‑1 (MQ=255)
                                                                                      gACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCt                                   <  2:1617943/150‑1 (MQ=255)
                                                                                        ccccGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCt                                   >  1:1180796/1‑148 (MQ=255)
                                                                                         cccGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGaa                                                                                            >  1:364517/1‑90 (MQ=35)
                                                                                         cccGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACCGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACGCGGTCGACACCACTACCGCGACCCTGAC                                <  2:2109557‑M1/150‑1 (MQ=255)
                                                                                         cccGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCAAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGAC                                >  1:550904‑M1/1‑150 (MQ=255)
                                                                                          ccGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACG                               >  1:1658319‑M1/1‑150 (MQ=255)
                                                                                          ccGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACG                               >  1:600475‑M1/1‑150 (MQ=255)
                                                                                          ccGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACG                               >  1:786673‑M1/1‑150 (MQ=255)
                                                                                           cGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGCGCAACTGACCCCGAACCCGACCCCGGCTCACACCACGCTCAACGACTCGGTCGATGCCACCACCGCGACCCCGACGG                              >  2:1115320‑M1/1‑150 (MQ=255)
                                                                                           cGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGG                              >  1:115116‑M1/1‑150 (MQ=255)
                                                                                            ggccagtGACGTCTACAACAACGGTTCGACCCTCAGCGGCACCCTCGAGCATGCCACAGGCGGTAATATCGGGTACCTGACCCCTAATCCGACCCCGGCCCAGACCTCGATCAACGTCTCCGTCGATGCCCCCACCGCGACCCTGACGGc                             <  2:769688/144‑1 (MQ=11)
                                                                                             gCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAt                                                    <  1:2441838/126‑1 (MQ=37)
                                                                                              cGAATGACGTCTACAACAACGTTTAGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGCCCCCGCGCCCGCCCCAGGCTCAGCCCACCCTCAACCACTCGTGCTATGCCACCACCGCGACCCTGACGGCGa                           <  2:1207239/150‑1 (MQ=11)
                                                                                               gAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCGGCGACCCGACCCCGGCTCAGACCACGATCACCGACTCGCTCGATGCCACCGCCGCGACCCTGACGGCGAg                          >  1:1614284/1‑150 (MQ=11)
                                                                                               gAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAg                          <  1:576465/150‑1 (MQ=255)
                                                                                                aaTGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTTCATCCCACCACCGCGACCCTGACGGCGAGc                         <  1:1023866/150‑1 (MQ=11)
                                                                                                aaTGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAt                                                    <  2:681418/123‑1 (MQ=34)
                                                                                                 aTGACGTCTACAGCAACGGTTCCACTGTCAGCGTCACTATCGAGAACGACACAGGCGGTAATTTCAAGCAACTGACCCCGGACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCc                        >  2:269416‑M2/1‑149 (MQ=255)
                                                                                                 aTGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCc                        >  2:1465755‑M2/1‑149 (MQ=255)
                                                                                                  tGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCcc                       >  1:1799469‑M2/1‑148 (MQ=255)
                                                                                                    ccgTCTACAACAACTGTTCGACCGTCAGCGCCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCccgt                     >  1:1931770‑M2/2‑146 (MQ=255)
                                                                                                     cGTCTACAACAACGGTTGGGCTGTCAGCGTCGCCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAGCCCGACCCCGGCCCACACCACGCCCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCccgtc                    >  2:1047752‑M2/1‑145 (MQ=255)
                                                                                                      gTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCg                   >  1:940878/1‑150 (MQ=38)
                                                                                                        cTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTACTTTCGAGCAACTGTCCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGgctcgccgtcagg                 <  2:1305621‑M1/150‑14 (MQ=255)
                                                                                                        cTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGt                 <  2:1549475/150‑1 (MQ=39)
                                                                                                         tACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGgctcgccgtcaggg                <  2:1390341‑M1/150‑15 (MQ=255)
                                                                                                           caacaaCGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGgctcgccgtcagg                 >  1:756663‑M1/1‑134 (MQ=255)
                                                                                                           caacaaCGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTc                                                       <  2:1021986/109‑1 (MQ=33)
                                                                                                           caacaaCGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGG                              <  1:2012758‑M1/134‑1 (MQ=255)
                                                                                                            aacaacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCa               >  1:234399/1‑148 (MQ=39)
                                                                                                            aacaacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCAcc             >  1:1397379/1‑150 (MQ=39)
                                                                                                            aacaacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCAcc             >  1:1885577/1‑150 (MQ=39)
                                                                                                              caacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGa           >  2:1675602/1‑150 (MQ=39)
                                                                                                              caacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGc                                                  >  2:607565/1‑111 (MQ=25)
                                                                                                               aacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAt                                                    <  1:2191763/108‑1 (MQ=35)
                                                                                                                 cGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCACACCGCGATCAACGACTCGGGCCATGCCACCCCCGCGACCCTGCCGCCGAGCCCGCCGGTCACCGCAgg        >  2:914195/1‑150 (MQ=11)
                                                                                                                  ggTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggc       <  1:652214/150‑1 (MQ=38)
                                                                                                                     tCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc    <  2:413431/150‑1 (MQ=38)
                                                                                                                          |                                                                                                                                                  
GCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGT  >  NC_002947/196910‑197178

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