Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 202,309 C→G 6.4% V2605V (GTC→GTG PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947202,3090CG 6.4% 255.5 / 7.9 93V2605V (GTC→GTGPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (40/47);  new base G (2/4);  total (42/51)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

TGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAAC  >  NC_002947/202165‑202458
                                                                                                                                                  |                                                                                                                                                     
tGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCgc                                                                                                                                                  >  1:967824/1‑149 (MQ=25)
   ccAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTGAccaaggg                                                                                                                                               <  2:647916/150‑6 (MQ=255)
    cAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGCCGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAggg                                                                                                                                              <  1:1616295/150‑2 (MQ=25)
    cAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGc                                                                                                                                              <  2:461183/150‑1 (MQ=35)
    cAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGc                                                                                                                                              <  1:1836282/150‑1 (MQ=35)
    cAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACAACATCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGc                                                                                                                                              >  1:115336/1‑150 (MQ=25)
     aGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCg                                                                                                                                             <  1:627210/150‑1 (MQ=35)
         cAGCCTGGGCATCACGTCCGCTGTGGATTTAGACGGTCCCATAATGGGGAAC‑‑GCCATTTGGCATACCGATGTCTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTga                                                                                                                                         <  1:223704/150‑1 (MQ=2)
           gCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTgaga                                                                                                                                       >  1:731837/1‑150 (MQ=35)
            ccTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAggggtcg                                                                                                                                          >  2:86542/1‑141 (MQ=255)
            ccTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAggggtcg                                                                                                                                          <  1:86542/146‑6 (MQ=255)
             cTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATc                                                                                                                                     <  2:791574/150‑1 (MQ=35)
               gggCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCAc                                                                                                                                   <  1:1246581/150‑1 (MQ=33)
                ggCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGCAGGCGGTCCGATCAcc                                                                                                                                  <  2:967824/150‑1 (MQ=12)
                      aCGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGA‑‑GTGAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACAcc                                                                                                                            >  1:1078179/1‑150 (MQ=12)
                       cGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCa                                                                                                                           <  1:1048869/150‑1 (MQ=32)
                        gTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCAt                                                                                                                          <  2:1305869/150‑1 (MQ=32)
                           cGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCAc                                                                                                                       >  1:44483/1‑150 (MQ=32)
                              tgtgGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTg                                                                                                                                          <  1:1715257/128‑1 (MQ=33)
                              tgtgGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTg                                                                                                                                          >  2:1715257/1‑128 (MQ=33)
                              tgtgGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCt                                                                                                                    <  2:731837/150‑1 (MQ=32)
                               gtgGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTg                                                                                                                   >  2:1376393/1‑150 (MQ=32)
                                tgGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGa                                                                                                                  <  1:1307611/150‑1 (MQ=32)
                                 ggATGTAGACGGTCACACCTTCGAGAACCTGCCCCCGGCC‑‑GGAAACACCTCGCTACCAGTGACCGACAGCTGCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGAc                                                                                                                 <  2:1653539/150‑1 (MQ=11)
                                 ggATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGAc                                                                                                                 <  2:1172518/150‑1 (MQ=32)
                                    tGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCCCCGAAGGCGGTGAGATCACCTACACCATCACCCTGACcaa                                                                                                              <  1:981806/150‑1 (MQ=17)
                                     gTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTTGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGaa                                                                                                                                                 >  1:426401/1‑114 (MQ=17)
                                     gTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGaa                                                                                                                                                 <  2:426401/114‑1 (MQ=33)
                                     gTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACcaac                                                                                                             <  2:1278225/150‑1 (MQ=32)
                                        gACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCACCaaa                                                                                                          >  1:787599/1‑150 (MQ=17)
                                         aCGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACAGAACGCGGTGACAGCACCTACACCATCACCCTGACcaacacaa                                                                                                         >  2:858521/1‑147 (MQ=11)
                                           ggccacaCCATAGAGTACCTGCAACTGGGC‑‑GGGAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGAc                                                                                                       <  2:1052412/147‑1 (MQ=11)
                                               ccacCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTc                                                                                                   <  2:1867773/149‑1 (MQ=25)
                                               acacGTTCGAGATCCTGCCACGGGGC‑‑GGCAACGCTTCGGTTCTAGTGTGCTACACCCTCGACCAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGTCCAACAATGACGGTc                                                                                                   <  2:889227/150‑1 (MQ=11)
                                               acacCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTc                                                                                                   >  2:545496/1‑150 (MQ=35)
                                               acacCTTCGAGAACATGCAACTGGGC‑‑GGCAACGCTTCGGTACAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCCCCTACACCCTCACCCTGACCAACAAAGACGGTc                                                                                                   <  2:116953/150‑1 (MQ=255)
                                                 acCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTg                                                                                                 >  1:658143/1‑150 (MQ=35)
                                                  ccTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAACGCTCTCAGATCACCTACACCAGCACCCTGACCAACAAAGACGGTCTGc                                                                                                <  1:800085/150‑1 (MQ=11)
                                                   cTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGcc                                                                                               >  2:800085/1‑150 (MQ=35)
                                                   cTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGAGGTGAGATCACCTACACCACCACCCTGACCAACAAAGACGGTCTGcc                                                                                               >  2:721861/1‑150 (MQ=21)
                                                   cTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTCGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGcc                                                                                               >  2:1408896/1‑150 (MQ=25)
                                                           aCCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATcaaca                                                                                       >  1:1740905/1‑150 (MQ=35)
                                                            ccTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTTTCGGTGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATcaacaa                                                                                      <  1:1772277/150‑1 (MQ=255)
                                                            ccTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATcaacaa                                                                                      <  2:1740905/150‑1 (MQ=35)
                                                             cTGGAACTGAGC‑‑GGCAACGTGTCAGTTCAAGTGACCGACTCCCTCGACGAGGTGGTGGCGAAGCTGGCCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACCAAGACGGTCTGCCGATCaacaac                                                                                     <  2:1444594/150‑1 (MQ=11)
                                                                 aaCTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACt                                                                                 >  1:1917005/1‑150 (MQ=35)
                                                                 aaCTGGGC‑‑GGCAACGCTGCGGTTCAAGTGACCGACACCCTCGGCGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGACCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACt                                                                                 <  1:875723/150‑1 (MQ=255)
                                                                   cTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCGGACCAACAAGGACGGTCTGCCGATCAACAACCACTCg                                                                               <  2:690580/150‑1 (MQ=21)
                                                                    tGGGC‑‑GGCAACGCTTCGGTTGAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCgg                                                                              >  1:965315/1‑150 (MQ=25)
                                                                    tGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCAcagtg                                                                              >  1:1575983/1‑146 (MQ=25)
                                                                      ggc‑‑ggcAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAg                                                                            <  1:32146/150‑1 (MQ=34)
                                                                      ggc‑‑ggcAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAg                                                                            <  2:186105/150‑1 (MQ=35)
                                                                            gcAACGCATCGGTGCCAGTGAGCAACACCGTCCACGGGGGGGTGGCGAAGCTGACCGCGACCCCTGCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGCCCAACAAAGACGGTCTGCCTATCAACAACCACTCGGAGCTGt                                                                        <  2:115336/150‑1 (MQ=11)
                                                                                    tCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGACGCCCAACCGGAGCTCTGCTTTCACCGAAGGCGGTGATGTCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAg                                                                <  1:1376393/150‑1 (MQ=11)
                                                                                       gTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAgggg                                                                                                                                             <  2:369758/70‑3 (MQ=255)
                                                                                       gTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAgggg                                                                                                                                             >  1:369758/1‑68 (MQ=255)
                                                                                       gTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTg                                                             >  2:392588/1‑150 (MQ=35)
                                                                                       gTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTg                                                             >  2:192146/1‑150 (MQ=35)
                                                                                            aGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCt                                                              <  2:931148/144‑1 (MQ=35)
                                                                                            aGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCt                                                              >  1:931148/1‑144 (MQ=35)
                                                                                            aGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGATGACAGa                                                        >  1:393027/1‑150 (MQ=21)
                                                                                              tGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACg                                                      >  2:124797/1‑150 (MQ=35)
                                                                                                    acacCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGt                                                                                                                                           >  2:465510/1‑59 (MQ=35)
                                                                                                    acacCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGt                                                                                                                                           <  1:465510/59‑1 (MQ=35)
                                                                                                     cacCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAATACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGcaccac                                               <  2:787599/150‑1 (MQ=25)
                                                                                                        ccTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACcgt                                            >  1:283258/1‑150 (MQ=38)
                                                                                                        ccTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTCCCGCTCAACAACCACTCGGAGCTGTACTGCAAGCTGACCGACGGCACCACcgt                                            >  2:1155313/1‑150 (MQ=255)
                                                                                                         cTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACcgtc                                           <  1:392588/150‑1 (MQ=38)
                                                                                                         cTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGCCAAGCCGGTAAACAACCACTCGGGGCTGTACTTCAAGCTGACCGACGGCACCACcgtc                                           >  2:529180/1‑150 (MQ=11)
                                                                                                             acgaGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACcgtcgtcg                                       <  1:1119356/150‑1 (MQ=38)
                                                                                                                aGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTgg                                    <  2:154339/150‑1 (MQ=38)
                                                                                                                  gtggtgGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCa                                  >  2:259586/1‑150 (MQ=39)
                                                                                                                  gtggtgGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCa                                  <  2:658143/150‑1 (MQ=39)
                                                                                                                   tggtggCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAg                                 >  2:604974/1‑150 (MQ=39)
                                                                                                                     gtggCGAAGCTGACCGCGCCCCCTTCGGTCAACGAAGGCGGTGAGATCACTTATACCATCCCGCTGACCAACAAATACGGTCTGCCGATCAACAACCACTCGGCGCTTTACGTCAAGCTGACCGACGGCGCCACCGTCGTCGTGGCAGcc                               <  2:1805697/150‑1 (MQ=11)
                                                                                                                       ggCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACcaac                                                                                                             <  2:1758342/70‑1 (MQ=37)
                                                                                                                       ggCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACcaac                                                                                                             >  1:1758342/1‑70 (MQ=37)
                                                                                                                          gAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAg                          >  1:1032694/1‑150 (MQ=39)
                                                                                                                          gAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAg                          >  2:1181426/1‑150 (MQ=39)
                                                                                                                           aaGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGc                         <  2:1032694/150‑1 (MQ=38)
                                                                                                                             gCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGcac                       >  2:360168/1‑150 (MQ=39)
                                                                                                                                  ccGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCg                  <  1:1181426/150‑1 (MQ=37)
                                                                                                                                   cgcgACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCgg                 >  2:1464106/1‑150 (MQ=35)
                                                                                                                                    gcgACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCTATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGt                >  1:1602117/1‑150 (MQ=25)
                                                                                                                                    gcgACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGt                >  1:1602340/1‑150 (MQ=34)
                                                                                                                                      gACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTc              <  2:237064/150‑1 (MQ=33)
                                                                                                                                      accgccTTCGGTCACCGAAGGGGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTc              >  1:1281560/5‑150 (MQ=17)
                                                                                                                                        ccccTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTgg                                    <  1:216375/126‑1 (MQ=33)
                                                                                                                                        ccccTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGGACTTCAAGCTGACCGACGGCACCACCGTCGTCGTgg                                    >  2:216375/1‑126 (MQ=18)
                                                                                                                                         cccTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGc           <  2:44483/150‑1 (MQ=32)
                                                                                                                                          ccTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCaacaac                                                                                     <  2:1576147/75‑1 (MQ=35)
                                                                                                                                          ccTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCaacaac                                                                                     >  1:1576147/1‑75 (MQ=35)
                                                                                                                                            ttCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCAc        <  1:259586/150‑1 (MQ=32)
                                                                                                                                                  cACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAAc  <  1:1464106/150‑1 (MQ=11)
                                                                                                                                                  |                                                                                                                                                     
TGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGC‑‑GGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAAC  >  NC_002947/202165‑202458

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: