Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,969,349 C→A 6.8% R78L (CGA→CTA)  cheA ← fused chemotactic sensory histidine kinase in two‑component regulatory system with CheB and CheY: sensory histidine kinase/signal sensing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,969,3490CA 6.8% 1306.9 / 51.0 338R78L (CGA→CTA) cheAfused chemotactic sensory histidine kinase in two‑component regulatory system with CheB and CheY: sensory histidine kinase/signal sensing protein
Reads supporting (aligned to +/- strand):  ref base C (148/167);  new base A (15/8);  total (163/175)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 340 total aligned reads displayed.

TGCCAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCGCTTCCGGCTGCAAAA  >  CP009273/1969201‑1969496
                                                                                                                                                    |                                                                                                                                                   
tGCCAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTc                                                                                                                                                     <  1:2980513/149‑1 (MQ=255)
 gCCAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCg                                                                                                                                                    >  2:4386304/1‑149 (MQ=255)
   cAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTc                                                                                                                                                  <  1:4182479/149‑1 (MQ=255)
   cAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTc                                                                                                                                                  <  1:5133382/149‑1 (MQ=255)
   cAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTc                                                                                                                                                  >  2:1405013/1‑149 (MQ=255)
   cAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTAGTc                                                                                                                                                  <  1:4683943/149‑1 (MQ=255)
    aGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCt                                                                                                                                                 <  1:4283609/149‑1 (MQ=255)
    aGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCt                                                                                                                                                 <  1:1686991/149‑1 (MQ=255)
    aGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCt                                                                                                                                                 <  1:1914025/149‑1 (MQ=255)
       tGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGc                                                                                                                                              <  1:1289649/149‑1 (MQ=255)
              aGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATTATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCg                                                                                                                                       >  2:766955/1‑149 (MQ=255)
              aGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCg                                                                                                                                       >  1:2432838/1‑149 (MQ=255)
               ggCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGa                                                                                                                                      >  2:888508/1‑149 (MQ=255)
               ggCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGa                                                                                                                                      >  1:3971587/1‑149 (MQ=255)
               ggCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGa                                                                                                                                      <  1:3352383/149‑1 (MQ=255)
               ggCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGa                                                                                                                                      >  2:299319/1‑149 (MQ=255)
               ggCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGa                                                                                                                                      >  1:1374306/1‑149 (MQ=255)
                gCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAg                                                                                                                                     >  2:3442539/1‑149 (MQ=255)
                  cTGGCAGATATAATCGACGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCa                                                                                                                                   <  2:4953890/149‑1 (MQ=255)
                  cTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCa                                                                                                                                   >  2:644352/1‑149 (MQ=255)
                  cTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCa                                                                                                                                   >  2:3951052/1‑149 (MQ=255)
                   tGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAg                                                                                                                                  <  2:3769023/149‑1 (MQ=255)
                   tGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAg                                                                                                                                  <  1:1339521/149‑1 (MQ=255)
                        gATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGtttt                                                                                                                             >  1:1095363/1‑149 (MQ=255)
                        gATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGtttt                                                                                                                             >  2:2690791/1‑149 (MQ=255)
                        gATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGtttt                                                                                                                             >  2:1608194/1‑149 (MQ=255)
                            taATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGGGCAGGTTTTCCAt                                                                                                                         >  2:2104306/1‑149 (MQ=255)
                            taATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCAt                                                                                                                         >  1:582613/1‑149 (MQ=255)
                            taATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCAt                                                                                                                         <  2:3528057/149‑1 (MQ=255)
                                cGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGa                                                                                                                     <  1:1620631/149‑1 (MQ=255)
                                cGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGa                                                                                                                     <  2:2583261/149‑1 (MQ=255)
                                 gAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGAt                                                                                                                    <  1:1939830/149‑1 (MQ=255)
                                  aaGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATg                                                                                                                   <  2:4076761/149‑1 (MQ=255)
                                  aaGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGCTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATg                                                                                                                   <  1:2694748/149‑1 (MQ=255)
                                       ggcggcATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTgg                                                                                                              <  1:794481/149‑1 (MQ=255)
                                        gcggcATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGt                                                                                                             <  2:1985528/149‑1 (MQ=255)
                                          ggcATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTAGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGttt                                                                                                           <  1:2597808/149‑1 (MQ=255)
                                           gcATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTc                                                                                                          >  2:552887/1‑149 (MQ=255)
                                            cATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTcc                                                                                                         <  2:1241247/149‑1 (MQ=255)
                                            cATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTcc                                                                                                         <  1:5014225/149‑1 (MQ=255)
                                             aTCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCt                                                                                                        <  1:5156292/149‑1 (MQ=255)
                                              tCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTg                                                                                                       >  1:2872296/1‑149 (MQ=255)
                                                cGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCa                                                                                                     >  1:2391812/1‑149 (MQ=255)
                                                  gCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCaaa                                                                                                   >  1:2483966/1‑149 (MQ=255)
                                                    tctTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAAc                                                                                                 >  2:1184725/1‑149 (MQ=255)
                                                    tctTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAAc                                                                                                 <  1:1036398/149‑1 (MQ=255)
                                                        gCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTg                                                                                             <  2:1671188/149‑1 (MQ=255)
                                                           aCTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAg                                                                                          <  1:3808459/149‑1 (MQ=255)
                                                            cTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGc                                                                                         <  2:4740938/149‑1 (MQ=255)
                                                            cTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGc                                                                                         >  1:3538224/1‑149 (MQ=255)
                                                            cTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGc                                                                                         <  2:1331030/149‑1 (MQ=255)
                                                            cTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGc                                                                                         <  2:3928664/149‑1 (MQ=255)
                                                            cTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTAGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGc                                                                                         >  2:3252710/1‑149 (MQ=255)
                                                             tGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGcc                                                                                        <  2:5176277/149‑1 (MQ=255)
                                                              gTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCa                                                                                       <  2:1104063/149‑1 (MQ=255)
                                                               tttATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCaa                                                                                      <  2:5122476/149‑1 (MQ=255)
                                                                   taAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTAGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGt                                                                                  >  1:4795657/1‑149 (MQ=255)
                                                                           gAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGccc                                                                          >  1:4506599/1‑149 (MQ=255)
                                                                            aGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGcccc                                                                         >  1:4096050/1‑149 (MQ=255)
                                                                              cTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCcctc                                                                       <  1:2253046/149‑1 (MQ=255)
                                                                              cTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCcctc                                                                       <  1:4857605/149‑1 (MQ=255)
                                                                                gTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCcctcct                                                                     <  2:1010642/149‑1 (MQ=255)
                                                                                gTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCcctcct                                                                     >  2:1821359/1‑149 (MQ=255)
                                                                                   cgtGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTg                                                                  <  2:4630745/147‑1 (MQ=255)
                                                                                    ttGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGa                                                                 <  2:2043873/149‑1 (MQ=255)
                                                                                    ttGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCACCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGa                                                                 <  2:1911950/149‑1 (MQ=255)
                                                                                     tGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGAt                                                                <  2:4680325/149‑1 (MQ=255)
                                                                                      gCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATc                                                               <  1:3606538/149‑1 (MQ=255)
                                                                                      gCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATc                                                               <  1:1009536/149‑1 (MQ=255)
                                                                                      gCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATc                                                               >  1:2755372/1‑149 (MQ=255)
                                                                                      gCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATc                                                               <  1:4341089/149‑1 (MQ=255)
                                                                                         tgatgTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAg                                                            >  1:1998152/1‑149 (MQ=255)
                                                                                            tgTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTgg                                                         <  1:532065/149‑1 (MQ=255)
                                                                                               ccTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCa                                                      <  1:1405013/149‑1 (MQ=255)
                                                                                                    gTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCg                                                 >  1:3637977/1‑149 (MQ=255)
                                                                                                    gTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCg                                                 >  1:753735/1‑149 (MQ=255)
                                                                                                    gTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTAGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCg                                                 >  2:3626553/1‑149 (MQ=255)
                                                                                                      ttCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGaa                                               <  2:3983500/149‑1 (MQ=255)
                                                                                                        ccAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAg                                             >  1:860676/1‑149 (MQ=255)
                                                                                                                 gATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTc                                    >  1:3244513/1‑149 (MQ=255)
                                                                                                                     gATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAAtt                                <  2:4832528/149‑1 (MQ=255)
                                                                                                                     gATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAAtt                                <  2:4555698/149‑1 (MQ=255)
                                                                                                                       tAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGt                              >  2:4775916/1‑149 (MQ=255)
                                                                                                                       tAATGTCGGTGTTGAGTTGCATCTCACCTAGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGt                              <  1:3922770/149‑1 (MQ=255)
                                                                                                                         aTGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTc                            <  1:1561036/149‑1 (MQ=255)
                                                                                                                          tGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCg                           >  1:916205/1‑149 (MQ=255)
                                                                                                                              ggTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCAt                       >  2:3853049/1‑149 (MQ=255)
                                                                                                                                  ttGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTgg                   >  2:1562568/1‑149 (MQ=255)
                                                                                                                                   tGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGgc                  <  2:4810951/149‑1 (MQ=255)
                                                                                                                                      gTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCCAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATGCAATTGTTCGGCATCTGGCGCt               >  2:2792932/1‑149 (MQ=255)
                                                                                                                                         gCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCGCTTcc            <  2:31001/149‑1 (MQ=255)
                                                                                                                                           aTCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCGCTTCCgg          <  1:2630575/149‑1 (MQ=255)
                                                                                                                                           aTCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCGCTTCCgg          >  1:4612371/1‑149 (MQ=255)
                                                                                                                                            tctcACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCGCTTCCGGc         <  2:2506953/149‑1 (MQ=255)
                                                                                                                                            tctcACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCGCTTCCGGc         <  2:1728249/149‑1 (MQ=255)
                                                                                                                                            tctcACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTGTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCGCTTCCGGc         <  2:4136558/149‑1 (MQ=255)
                                                                                                                                            tctcACCTAGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCGCTTCCGGc         >  1:2891545/1‑149 (MQ=255)
                                                                                                                                              tcACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCGCTTCCGGCTg       >  1:3390456/1‑149 (MQ=255)
                                                                                                                                              tcACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTACTGCCCCACCTTTGATAGAGTGGGCAGCCCGAAATAGGGCATTCAATTGTTCGGCATCTGGCGCTTCCGGCTg       >  1:310902/1‑149 (MQ=255)
                                                                                                                                                   tCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCGCTTCCGGCTGCaaaa  <  1:4979304/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGCCAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCGCTTCCGGCTGCAAAA  >  CP009273/1969201‑1969496

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