Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,424,254 G→T 10.4% intergenic (‑12/+247) cvpA ← / ← dedD colicin V production membrane protein/membrane‑anchored periplasmic protein involved in septation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,424,2540GT10.4% 1082.7 / 63.5 298intergenic (‑12/+247)cvpA/dedDcolicin V production membrane protein/membrane‑anchored periplasmic protein involved in septation
Reads supporting (aligned to +/- strand):  ref base G (137/130);  new base T (17/14);  total (154/144)
Fisher's exact test for biased strand distribution p-value = 8.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 300 total aligned reads displayed.

GTAGTAATGACTGGCAACAAAGAAAGCACAACCCCATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACTTTTTATGCCTTTGCGGCATCGGGCAAT  >  CP009273/2424108‑2424403
                                                                                                                                                  |                                                                                                                                                     
  agtaATGACTGGCAACAAAGAAAGCACAACCCCATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGAtttt                                                                                                                                                  >  2:2319727/1‑150 (MQ=255)
    taATGACTGGCAACAAAGAAAGCACAACCCCATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTAc                                                                                                                                                <  1:1090003/150‑1 (MQ=255)
      aTGACTGGCAACAAAGAAAGCACAACCCCATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGc                                                                                                                                              >  1:356468/1‑150 (MQ=255)
      aTGACTGGCAACAAAGAAAGCACAACCCCATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGc                                                                                                                                              >  1:4260327/1‑150 (MQ=255)
      aTGACTGGCAACAAAGAAAGCACAACCCCATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGc                                                                                                                                              >  1:148511/1‑150 (MQ=255)
      aTGACTGGCAACAAAGAAAGCACAACCCCATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGc                                                                                                                                              >  1:3097802/1‑150 (MQ=255)
             gCAACAAAGAAAGCACAACCCCATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCa                                                                                                                                       >  1:1217959/1‑150 (MQ=255)
              caacaaAGAAAGCACAACCCCATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCAt                                                                                                                                      >  1:2742960/1‑150 (MQ=255)
                acaaAGAAAGCACAACCCCATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATcc                                                                                                                                    <  1:114182/150‑1 (MQ=255)
                acaaAGAAAGCACAACCCCATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATcc                                                                                                                                    >  2:2442607/1‑150 (MQ=255)
                    agaaagCACAACCCCATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGtt                                                                                                                                >  2:1244095/1‑150 (MQ=255)
                         gCACATCCCCATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCATGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTACCACGATTTTACGCCCTGTCATCCTGttagtgg                                                                                                                           <  1:3078381/150‑6 (MQ=255)
                            caACCCCATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGgcgc                                                                                                                        >  2:1665279/1‑150 (MQ=255)
                             aaCCCCATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCa                                                                                                                       >  2:4361710/1‑150 (MQ=255)
                              aCCCCATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCAt                                                                                                                      >  2:925974/1‑150 (MQ=255)
                              aCCCCATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCAt                                                                                                                      >  2:3760429/1‑150 (MQ=255)
                               ccccATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCAtt                                                                                                                     >  2:3588845/1‑150 (MQ=255)
                               ccccATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCAtt                                                                                                                     <  1:17984/150‑1 (MQ=255)
                               ccccATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCAtt                                                                                                                     >  1:33970/1‑150 (MQ=255)
                                 ccATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCt                                                                                                                   <  2:2778290/150‑1 (MQ=255)
                                 ccATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCt                                                                                                                   >  1:3460642/1‑150 (MQ=255)
                                   aTGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACTATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTaa                                                                                                                 >  2:3532545/1‑150 (MQ=255)
                                   aTGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTaa                                                                                                                 >  2:1393400/1‑150 (MQ=255)
                                    tGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACTATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAAc                                                                                                                >  2:3172341/1‑150 (MQ=255)
                                      tCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACTATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACag                                                                                                              >  2:3101483/1‑150 (MQ=255)
                                      tCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACag                                                                                                              <  1:4605949/150‑1 (MQ=255)
                                        aCCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACagaa                                                                                                            <  1:2594721/150‑1 (MQ=255)
                                        aCCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACagaa                                                                                                            <  2:2081400/150‑1 (MQ=255)
                                        aCCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACagaa                                                                                                            <  1:1996536/150‑1 (MQ=255)
                                        aCCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACagaa                                                                                                            <  1:3301013/150‑1 (MQ=255)
                                           aGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACagaaaaa                                                                                                         >  2:1132524/1‑150 (MQ=255)
                                              gATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACTATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACagaaaaagaa                                                                                                      >  2:2236631/1‑150 (MQ=255)
                                               aTAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACTATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACagaaaaagaaa                                                                                                     >  2:3324345/1‑150 (MQ=255)
                                                tAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACagaaaaagaaaa                                                                                                    >  2:4010772/1‑150 (MQ=255)
                                                tAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACagaaaaagaaaa                                                                                                    <  2:4522342/150‑1 (MQ=255)
                                                       tCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGc                                                                                             <  2:4142608/150‑1 (MQ=255)
                                                        cACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCg                                                                                            >  1:3891100/1‑150 (MQ=255)
                                                        cACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCg                                                                                            >  1:3960460/1‑150 (MQ=255)
                                                        cACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCg                                                                                            >  2:1917774/1‑150 (MQ=255)
                                                        cACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCg                                                                                            <  2:846191/150‑1 (MQ=255)
                                                             aCAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAggg                                                                                       <  1:2236561/150‑1 (MQ=255)
                                                             aCAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAggg                                                                                       <  2:766475/150‑1 (MQ=255)
                                                              cAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCTATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGa                                                                                      <  2:1718805/150‑1 (MQ=255)
                                                                 aGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACTATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGAttt                                                                                   <  1:680185/150‑1 (MQ=255)
                                                                 aGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGAttt                                                                                   >  1:4086078/1‑150 (MQ=255)
                                                                   ccGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATttcc                                                                                 <  2:3931443/150‑1 (MQ=255)
                                                                          tCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACTATTTTACGCCCTGTCATCCTGTTCGGGGCGCACTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCcgc                                                                          >  2:4145082/1‑150 (MQ=255)
                                                                             ggCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCa                                                                       <  2:571456/150‑1 (MQ=255)
                                                                              gCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCAt                                                                      <  1:1905535/150‑1 (MQ=255)
                                                                              gCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCAt                                                                      <  2:1233894/150‑1 (MQ=255)
                                                                                tAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCa                                                                    >  2:1868156/1‑150 (MQ=255)
                                                                                    cAGAGAGGAAAAAGCAATCACCGCGATTATTGCGTAATCAATCCAGACCATATGTGTCCCACGCTTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCGTCCCGCGCATCAATaa                                                                <  1:2147349/150‑1 (MQ=255)
                                                                                    cAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATaa                                                                <  2:894759/150‑1 (MQ=255)
                                                                                     agagagGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATaaa                                                               >  2:4605325/1‑150 (MQ=255)
                                                                                       agagGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATaaaaa                                                             >  1:1031651/1‑150 (MQ=255)
                                                                                          ggAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATgg                                                          <  2:26722/150‑1 (MQ=255)
                                                                                           gAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGgc                                                         >  1:3067070/1‑150 (MQ=255)
                                                                                            aaaaaGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGgcg                                                        <  1:4323460/150‑1 (MQ=255)
                                                                                             aaaaGCAATCACCGCGATTATGGCGTAATCAATCCAGACAATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTAGGTGCGCATTCTAACAGAAAAAGACAACTTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGtcgc                                                       >  2:1968764/1‑150 (MQ=255)
                                                                                               aaGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTg                                                     >  2:4573557/1‑150 (MQ=255)
                                                                                                aGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGa                                                    <  1:2354942/150‑1 (MQ=255)
                                                                                                 gCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCACTAAAAATGGCGCTGaa                                                   >  1:2418098/1‑150 (MQ=255)
                                                                                                 gCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGaa                                                   >  1:3407575/1‑150 (MQ=255)
                                                                                                 gCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGaa                                                   >  2:1144845/1‑150 (MQ=255)
                                                                                                   aaTCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGaaaa                                                 >  1:2013392/1‑150 (MQ=255)
                                                                                                    aTCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGaaaaa                                                >  1:4554881/1‑150 (MQ=255)
                                                                                                      cACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAat                                              <  2:4340395/150‑1 (MQ=255)
                                                                                                      cACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAat                                              <  1:2933989/150‑1 (MQ=255)
                                                                                                      cACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAat                                              <  2:3220557/150‑1 (MQ=255)
                                                                                                       aCCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAata                                             <  1:3096405/150‑1 (MQ=255)
                                                                                                        ccGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAatat                                            <  1:358741/150‑1 (MQ=255)
                                                                                                           cgATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCa                                         >  1:785238/1‑150 (MQ=255)
                                                                                                                 tGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCa                                   >  1:3526214/1‑150 (MQ=255)
                                                                                                                    cGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCg                                <  1:373062/150‑1 (MQ=255)
                                                                                                                    cGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCg                                <  2:4541410/150‑1 (MQ=255)
                                                                                                                    cGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCg                                <  2:3939979/150‑1 (MQ=255)
                                                                                                                     gcaatcaatcCAGACCATATGTGTCCCACGATTTTACTCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGa                               >  1:2774568/3‑150 (MQ=255)
                                                                                                                     gTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGa                               <  1:1343745/150‑1 (MQ=255)
                                                                                                                       aatcaatcCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGaat                             >  2:3605564/1‑148 (MQ=255)
                                                                                                                       aatcaatcCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACt                             <  1:2846361/150‑1 (MQ=255)
                                                                                                                       aatcaatcCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACt                             <  1:534275/150‑1 (MQ=255)
                                                                                                                         tcaatcCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACttt                           <  2:181275/150‑1 (MQ=255)
                                                                                                                          caatcCAGACCATATGTGTCCCACTATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACtttt                          >  1:1451719/1‑150 (MQ=255)
                                                                                                                           aatcCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACttttt                         >  1:3204640/1‑150 (MQ=255)
                                                                                                                           aatcCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACttttt                         >  1:607675/1‑150 (MQ=255)
                                                                                                                           aatcCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCACGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACttttt                         >  2:4350444/1‑150 (MQ=255)
                                                                                                                             tcCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACTTTTTAt                       <  2:3983810/150‑1 (MQ=255)
                                                                                                                             tcCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACTTTTTAt                       <  1:3836967/150‑1 (MQ=255)
                                                                                                                                aGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACTTTTTATGcc                    >  2:3491139/1‑150 (MQ=255)
                                                                                                                                  aCCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACTTTTTATGCCtt                  <  2:4200051/150‑1 (MQ=255)
                                                                                                                                   ccATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGTGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACTTTTTATGCCttt                 <  2:1555435/150‑1 (MQ=255)
                                                                                                                                    cATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACTTTTTATGCCTTTg                >  2:2284368/1‑150 (MQ=255)
                                                                                                                                    cATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACTTTTTATGCCTTTg                <  1:3152502/150‑1 (MQ=255)
                                                                                                                                     atatGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACTTTTTATGCCTTTGc               >  2:1013937/1‑150 (MQ=255)
                                                                                                                                     atatGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACTTTTTATGCCTTTGc               <  2:2039833/150‑1 (MQ=255)
                                                                                                                                        tgtgTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACTTTTTATGCCTTTGCGGc            >  1:4344705/1‑150 (MQ=255)
                                                                                                                                         gtgtCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACTTTTTATGCCTTTGCGGCa           >  1:1010900/1‑150 (MQ=255)
                                                                                                                                           gtCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACTTTTTATGCCTTTGCGGCATc         <  2:1081280/150‑1 (MQ=255)
                                                                                                                                               cACTATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACTTTTTATGCCTTTGCGGCATCGGGc     <  2:3003484/150‑1 (MQ=255)
                                                                                                                                                  gATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACTTTTTATGCCTTTGCGGCATCGGGCAAt  <  2:3812825/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTAGTAATGACTGGCAACAAAGAAAGCACAACCCCATGTCACCAGCGATAACGCTTCACGAACAAAGCCGCGGATCAGGCTAACCAGAGAGGAAAAAGCAATCACCGCGATTATGGCGTAATCAATCCAGACCATATGTGTCCCACGATTTTACGCCCTGTCATCCTGTTCGGGGCGCATTCTAACAGAAAAAGAAAACGTTTGCGTAGGGATTTCCTTCCCGCGCATCAATAAAAATGGCGCTGAAAAAATATTCAACGCCATCGACTTTTTATGCCTTTGCGGCATCGGGCAAT  >  CP009273/2424108‑2424403

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: