Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,685,967 G→C 6.2% R328G (CGC→GGC) ‡ PP_4988 ← signal transduction histidine kinase CheA/CheY‑like receiver

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,685,9670GC 6.2% 405.0 / 14.2 128R328G (CGC→GGC) ‡PP_4988signal transduction histidine kinase CheA/CheY‑like receiver
Reads supporting (aligned to +/- strand):  ref base G (65/55);  new base C (3/5);  total (68/60)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 128 total aligned reads displayed.

TCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTTCCTCCAAAGGGCC  >  NC_002947/5685822‑5686114
                                                                                                                                                 |                                                                                                                                                   
tcaATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGc                                                                                                                                                 <  1:1362650/150‑1 (MQ=255)
 caATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACCGCCCCGGCAAGGGCTGCAGGCGTTGcc                                                                                                                                                >  2:543796/1‑150 (MQ=255)
  aaTCGGCCGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCCATTTGCGCGCGCAGCCGGCTAATAAGTGGGCCACGGCCGGCAGCATCAACATGGGCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGccc                                                                                                                                               >  2:82109/1‑149 (MQ=255)
   aTCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAGCCCTTGCCgg                                                                                                                                              <  2:2135859/150‑1 (MQ=255)
   aTCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAGCCCTTGCCgg                                                                                                                                              <  1:2184256/150‑1 (MQ=255)
           ctctTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTc                                                                                                                                      <  1:1573213/150‑1 (MQ=255)
            tctTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCg                                                                                                                                     >  1:1675208/1‑150 (MQ=255)
            tctTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAGCCCTTGCCGGGGctgt                                                                                                                                        <  2:1284095/147‑4 (MQ=39)
            tctTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAGCCCTTGCCGGGGctgt                                                                                                                                        >  1:1284095/1‑144 (MQ=39)
            cctTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCg                                                                                                                                     <  2:2430754/149‑1 (MQ=255)
            acttcttCGAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCg                                                                                                                                     >  2:1174661/2‑150 (MQ=255)
              ttCCTCCGGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGa                                                                                                                                   <  2:805883/150‑1 (MQ=255)
              ttCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCGTCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGa                                                                                                                                   >  2:1186158/1‑150 (MQ=255)
              ttCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGCTCTGCAGGCGTTGCCGGGGCCAGTCGGa                                                                                                                                   >  1:2207165/1‑150 (MQ=255)
                cctccAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACt                                                                                                                                 >  1:222730/1‑150 (MQ=255)
                  tccAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTcc                                                                                                                               <  1:906424/150‑1 (MQ=255)
                      ggTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCgggg                                                                                                                                            >  1:1376676/1‑133 (MQ=255)
                      ggTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCgggg                                                                                                                                            <  2:1376676/133‑1 (MQ=255)
                          ctgcAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAg                                                                                                                             <  1:696955/144‑1 (MQ=255)
                          ctgcAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAg                                                                                                                             >  2:696955/1‑144 (MQ=255)
                          ctgcAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCGCGGCCACGGCCTGGAGCATCAACATGGGAAATGAATTAATGCTAAAGATTGACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATc                                                                                                                       >  1:1124788/1‑150 (MQ=255)
                           tgcAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCg                                                                                                                      >  2:730448/1‑150 (MQ=255)
                            gcAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGc                                                                                                                     <  1:2224793/150‑1 (MQ=255)
                              aaccaacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCg                                                                                                                   >  1:1187081/1‑150 (MQ=255)
                               accaacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGa                                                                                                                  <  2:2432417/150‑1 (MQ=255)
                                ccaacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGAt                                                                                                                 >  2:890217/1‑150 (MQ=255)
                                 caacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATc                                                                                                                >  1:2433295/1‑150 (MQ=255)
                                  aacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTcc                                                                                                                               >  2:2384361/1‑134 (MQ=255)
                                  aacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTcc                                                                                                                               <  1:2384361/134‑1 (MQ=255)
                                  aacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCa                                                                                                               >  1:1186535/1‑150 (MQ=255)
                                     gCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGAt                                                                                                            >  2:1915483/1‑150 (MQ=255)
                                       tGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATcc                                                                                                          >  1:182537/1‑150 (MQ=255)
                                        gCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCt                                                                                                         >  1:385485/1‑150 (MQ=255)
                                        gCACATAGCGTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCt                                                                                                         <  1:383915/150‑1 (MQ=255)
                                         cacaTAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAgg                                                                                                                            <  2:567544/130‑1 (MQ=255)
                                         cacaTAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAgg                                                                                                                            >  1:567544/1‑130 (MQ=255)
                                          acaTAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTg                                                                                                       >  1:2064442/1‑150 (MQ=255)
                                            aTAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGc                                                                                                     <  2:1675208/150‑1 (MQ=255)
                                            aTAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGc                                                                                                     <  2:530961/150‑1 (MQ=255)
                                               gCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTc                                                                                                                                      <  2:1716255/114‑1 (MQ=255)
                                               gCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTc                                                                                                                                      >  1:1716255/1‑114 (MQ=255)
                                                cTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTgg                                                                                                      <  1:1735398/145‑1 (MQ=255)
                                                     gcAGCCGACAAACAGCTCGGCCACGGCCTGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCaa                                                                                            >  1:1519632/1‑150 (MQ=255)
                                                         ccGACAAACAGCTCGGCCACGGCCGGTAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTgg                                                                                                      <  2:1360403/136‑1 (MQ=255)
                                                         ccGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACg                                                                                        <  2:451741/150‑1 (MQ=255)
                                                          cGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGCCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGt                                                                                       >  2:2192687/1‑150 (MQ=255)
                                                              aaaCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAg                                                                                   <  1:839686/150‑1 (MQ=255)
                                                               accAGCTCGGCGACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGc                                                                                  <  2:1997596/148‑1 (MQ=255)
                                                                aCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGa                                                                                                                  <  1:606003/117‑1 (MQ=255)
                                                                aCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGa                                                                                                                  >  2:606003/1‑117 (MQ=255)
                                                                aCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCaa                                                                                            >  1:465788/1‑139 (MQ=255)
                                                                aCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCaa                                                                                            <  2:465788/139‑1 (MQ=255)
                                                                aCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCg                                                                                 >  1:302956/1‑150 (MQ=255)
                                                                    cTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGattc                                                                             >  2:1412634/1‑150 (MQ=255)
                                                                     ccGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGattca                                                                            <  2:48068/149‑1 (MQ=255)
                                                                         ccACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGAttcattca                                                                        >  1:1600785/1‑150 (MQ=255)
                                                                           aCGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGc                                                                      >  2:1359018/1‑150 (MQ=255)
                                                                           aCGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGc                                                                      >  2:186431/1‑150 (MQ=255)
                                                                               ccggcAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCa                                                                  <  2:2064442/150‑1 (MQ=255)
                                                                                   cagcaTCAACCTGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTg                                                              <  2:1124788/150‑1 (MQ=255)
                                                                                         cAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATc                                                                                                                       <  1:2014246/87‑1 (MQ=255)
                                                                                             aTGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCg                                                    <  1:1412634/150‑1 (MQ=255)
                                                                                             aTGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCg                                                    <  2:1742185/150‑1 (MQ=255)
                                                                                               ggCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTTCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGc                                                  <  1:510642/150‑1 (MQ=255)
                                                                                                gCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGc                                                                                  <  2:1131176/117‑1 (MQ=255)
                                                                                                 ccAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCa                                                <  1:1328779/150‑1 (MQ=255)
                                                                                                  cAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAg                                               >  1:587528/1‑150 (MQ=255)
                                                                                                  cAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAg                                               >  1:1066470/1‑150 (MQ=255)
                                                                                                     tGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTc                                            <  2:385485/150‑1 (MQ=255)
                                                                                                     tGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTc                                            <  1:1186158/150‑1 (MQ=255)
                                                                                                          cgcgGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGACAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGc                                       >  2:2325064/1‑150 (MQ=255)
                                                                                                              gATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCgg                                                                                                                                              >  1:106155/1‑43 (MQ=255)
                                                                                                                   gctgcAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGc                              <  1:186431/150‑1 (MQ=255)
                                                                                                                    ctgcAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAg                                                                                                                             <  1:204301/54‑1 (MQ=255)
                                                                                                                    ctgcAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAg                                                                                                                             >  2:204301/1‑54 (MQ=255)
                                                                                                                    ctgcAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGa                                         <  1:1311502/138‑1 (MQ=255)
                                                                                                                    ctgcAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGa                                         >  2:1311502/1‑138 (MQ=255)
                                                                                                                     tgcAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTg                            >  2:1042146/1‑150 (MQ=255)
                                                                                                                     tgcAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTg                            >  1:472927/1‑150 (MQ=255)
                                                                                                                        aaCAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCt                         >  1:390666/1‑150 (MQ=255)
                                                                                                                         aCAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCtt                        <  1:730448/150‑1 (MQ=255)
                                                                                                                          aaGCCCCGGCAAGGGCTGCAGACGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCAGTCAGCACCAGCTGGCGCAGCTCTGCCAGGTCCGAGCCCGGCAGGCTGCCTTg                       >  2:1729620/2‑150 (MQ=255)
                                                                                                                           aGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGa                                         <  1:449657/131‑1 (MQ=255)
                                                                                                                           aGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGa                                         >  2:449657/1‑131 (MQ=255)
                                                                                                                           aGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTgg                      <  2:222730/150‑1 (MQ=255)
                                                                                                                            gCCCCGGCAACGCCTGCAGCCCTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGc                     >  2:676138/1‑150 (MQ=255)
                                                                                                                             ccccGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGcc                    >  1:176880/1‑150 (MQ=255)
                                                                                                                              cccGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGc                                                                                  <  1:1027816/87‑1 (MQ=255)
                                                                                                                              cccGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGc                                                                                  >  2:1027816/1‑87 (MQ=255)
                                                                                                                              cccGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCg                   >  1:1555049/1‑150 (MQ=255)
                                                                                                                               ccGGCAACGCCTGCAGCCCTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGt                  >  2:2440202/1‑150 (MQ=255)
                                                                                                                                cGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGtt                 >  2:549208/1‑150 (MQ=255)
                                                                                                                                     aGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTtcct            <  1:1876005/150‑1 (MQ=255)
                                                                                                                                     aGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTtcct            >  1:2036306/1‑150 (MQ=255)
                                                                                                                                        gCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGa                                         >  2:506287/1‑118 (MQ=255)
                                                                                                                                          tGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTTCCTCCaaa       >  2:709593/1‑150 (MQ=255)
                                                                                                                                           gCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTTCCTCCaaa       <  1:475063/149‑1 (MQ=255)
                                                                                                                                           gCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTTCCTCCAAAg      <  1:786532/150‑1 (MQ=255)
                                                                                                                                            cAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGc                                                                                  <  2:1274386/73‑1 (MQ=255)
                                                                                                                                            cAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGc                                                                                  >  1:1274386/1‑73 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTTCCTCCAAAGGGCC  >  NC_002947/5685822‑5686114

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