Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273748,5810TA32.7% 101.7 / 17.0 55intergenic (‑330/+60)ybgD/gltAfimbrial protein/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (30/7);  new base A (0/18);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 3.34e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAAAATATTGTTAGTTTTATTAAATGCAAAACTAAATTATTGGTATCATGAATTTGTTGTATGATGAATAAAATATAGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACC  >  NZ_CP009273/748336‑748826
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
aaaaaTATTGTTAGTTTTATTAAATGCAAAACTAAATTATTGGTATCATGAATTTGTTGTATGATGAATAAAATATAGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaa                                                                                                                                                                                                                                                    >  1:656858/1‑249 (MQ=255)
aaaaaTATTGTTAGTTTTATTAAATGCAAAACTAAATTATTGGTATCATGAATTTGTTGTATGATGAATAAAATATAGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATAGATTTTATTTTTACAATTTTTTGTGTTTAGGCatatg                                                                                                                                                                                                                                                      <  2:344994/247‑2 (MQ=255)
             gTTTTATTAAATGCAAAACTAAATTATTGGTATCATGAATTTGTTGTATGATGAATAAAATATAGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCa                                                                                                                                                                                                                                       <  1:566516/249‑1 (MQ=255)
               tttATTAAATGCAAAACTAAATTATTGGTATCATGAATTTGTTGTATGATGAATAAAATATAGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCatat                                                                                                                                                                                                                                    >  1:744159/1‑250 (MQ=255)
                                    tttttGGTATAATAAATTTGTTTTATGATGAAAAAAATAGGGGGGGGTATAGATAGACGTCTTTTTAATGGGGTTAAAAATGCGACTACCATGAATTTTTTAATGAAAATTTTTGGGTTGTTGATATTTTGATTTTTTCTTTTTTTTTTAAATATCTATATAGTTTTTTTAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGAAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATAttt                                                                                                                                                                                                               <  2:184025/247‑1 (MQ=255)
                                             taataatttttttGTATGATAAATAAAATATGGGGGGGTATAAATAGACGTCTTTTTCATAGGGTTAAAAATGCGACTACCATGAATTTTTTAATTGAAAGTTTTGGTTTGCTGATAATTTGAGCTGTTCTTTTTTTTTTAAATATCTATATAGGTTTGTTAATGGATTTTATTTTTACATTTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTa                                                                                                                                                                                                       <  1:372867/242‑1 (MQ=255)
                                                    tttGTTGTATGATGAATAAAATATAGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCaa                                                                                                                                                                                               >  1:747627/1‑250 (MQ=255)
                                                      tgttgtATGATGAATAAAATATAGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTTAAATATTTAAAACTTAACAaacaac                                                                                                                                                                                              >  2:787637/1‑249 (MQ=255)
                                                             tgatgaAAAAAAAATAGGGGGGTATAGATAGACGTAATTTTCATAGGGTTATAAATGCGACTACCATGAATTTTTTAATTGAAATTATTGGGTTGCTAATAATTTGAGCTTTTCTATTCTTTTTAAATATTTATATAGTTTTGTAAATGGTTTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACg                                                                                                                                                                                      <  1:799663/250‑1 (MQ=255)
                                                               ataaaaaaaaTATGGGGGGTTATATATAGACGTATTTTTCAGAGGGTTAAAATTGCGACTACCGTAATTTTTTTAATTGAAAGTATTGGTTTACTGATAATTTGATCTTTTATATTTTTTTAAAATTTTTATATGGGTGTGTAAATGTTTTTTATTTTAAAATTTTTTTTTTTTTAGAAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCt                                                                                                                                                                                    <  2:183698/244‑1 (MQ=255)
                                                                tgataaaaaTATGGGGGGGTATAGATATCTTCTATTTTCATTGTGTAATAAAGGCGATTACAATGAATTTTTAAATAGAAAGTATGGGGTTTCTGATATTTTTGTTTTTTTTTTTTTTTTTAAATATCTATATATTTTTGTAAAGGTTTTTTATTTTAAAATTTTTTTGTTTTAAGAAAAATAAAAAACACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCtt                                                                                                                                                                                   <  2:30645/245‑1 (MQ=255)
                                                                     aaaaTATAGGGGGGTATAGATAGACGTCTTTTTCATAGGGTTATAAATGCGATTCCCATAAATTTTTTAATTAAAATTTTGGGGTTGTTAATATTTTGAGTTTTTTTTTTTTTTTTAAATATTTATATAGGTTTTTAAATGGTTTTTTTTTTTAAATTTTTTTGTGTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGatat                                                                                                                                                                              <  1:306127/250‑1 (MQ=255)
                                                                        atatAGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                                                                                                                           >  1:683351/1‑250 (MQ=255)
                                                                           tgggggggTATAGATAGACGTCATTTTAATAGGTTTATAAATGCGACTACCATGAATTTTTTAATTGAAAGTTTTGGTTTGTTAATATTTTGAGTTTTTTTTTTTTTTTAAAATTTTTATATAGTTTTTTTAATGAATTTTATTTTAACATTTTTTTGTTTTAAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCttt                                                                                                                                                                        <  2:392997/248‑1 (MQ=255)
                                                                            gggggggTATAAATAGACCTTTTTTTCATTGTGTTAAAAATGCGAATACTATGAATTTTTTATTTAAAAGTATTTGTTTTCTTATATTTTGATGTTTTTTTTTCTTTTTAAATATTTATATGGGTCTGTTAATGGATTTTTTTTTTACATTTTTTTGTTTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCtttt                                                                                                                                                                       <  2:765618/249‑1 (MQ=255)
                                                                             gggggtAATAGATAGACGTCTTTTTAATAGGTTTATAAATGCGACTACCATAAATTTTTTAATTAAAATTTTGGGTTTGCTAATAATTTGAGCTTTTTTATTCTTTTAAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTAAAATTTTTTGTGTTTAGCCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTa                                                                                                                                                                      <  1:6359/250‑1 (MQ=255)
                                                                              gggggTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAATTATTTAAAAAATATAAATCAAACATGATAGCTTTTTATCGCtcttt                                                                                                                                                                      >  2:769477/1‑248 (MQ=255)
                                                                                gggTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAAATTAACAATCAAACCTTAAAGCTTGATATCGCTTTTAAAg                                                                                                                                                                   >  2:288452/1‑250 (MQ=255)
                                                                                 gttaagaatagaCTTAATTTTCATAGTTTTATAAATGCGACTACAATGATTTTTTTAATTAAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTTTTTTTTTAAAATATTTATATAGGTATTTAAATGTATTTTTTTTTTACAATTTTTTGTGTTTAGGCAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGt                                                                                                                                                                  <  2:249637/243‑1 (MQ=255)
                                                                                   tataGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTTAATATTTATAACTTAGCAATCAACCATTAACGCTTCATATCCCCTTTTAAgaca                                                                                                                                                                >  2:629542/1‑247 (MQ=255)
                                                                                                    tttCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAAAATAATTATAAATTAACAAACAACAAATAACGCTAGATAACGCTTTTAAAATTCCGTTTTTTATATTCTg                                                                                                                                                >  2:347585/1‑250 (MQ=255)
                                                                                                      tCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGta                                                                                                                                              >  1:223468/1‑250 (MQ=255)
                                                                                                              ttATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAAATTAGCAATCAACCATTAAAGCTTGTTATCGCTTTTTAAATAGCGGTTTTCATATCCTGTATAAAGCt                                                                                                                                       >  2:65797/1‑249 (MQ=255)
                                                                                                               tataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTTATATCCTGTATACAGCTGa                                                                                                                                     >  2:272247/1‑250 (MQ=255)
                                                                                                                ataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGAAATCAAACATTTACCCTTGATATTGCTTTTTAAGTCGCGCTTTTCTTATTCTCATTACAGCTgta                                                                                                                                    >  2:660398/1‑248 (MQ=255)
                                                                                                                   aaTGCAATTCCCATGATTTTTTTAATTAAATGTTTTGGTTTGTTGATATTTTGGGTTTTTCTTTTCTTTTAAAATATTTATATAGTTCTTTTAATGTTTTTTTTTTTTACATTTTTTTGTTTTAAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcgcg                                                                                                                                 <  2:191420/250‑1 (MQ=255)
                                                                                                                    aTGCGCCTACCATGAATTTTTTAATTAAAAGTATTGGGTTGTTGATAATTTGAGTTTTTCTTTTCTTTTTAAATATCTATATAGTTTTGTTAATGTTTTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCgg                                                                                                                                <  1:273369/250‑1 (MQ=255)
                                                                                                                               aTGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTATAACTTAGCAATAAACCATTAACGCTTGATATCTATTTTAAATACTCGTTTTTCTTTTCCTTTTTACAACCTCCGCGGACGGGGCATCt                                                                                                                     >  2:614582/1‑250 (MQ=255)
                                                                                                                                  agttttttAATGAAATTTATTGGGTTTTTAATATTTTGGTCTTTTTTTTTTTTTTAAAATATTTATATAGGTTTTTTATTGTATTTTATTTTTAAATTTTTTTGTTTTTAAGAATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCa                                                                                                                  <  2:737244/247‑1 (MQ=255)
                                                                                                                                  aaGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAAATTAGCAATCAAACATTAACGCTTGATATTGCTTTTTAAGTCGCGGTTTTTATATCCTGTATACAGCTGATGCGGACGGGCAATCTTct                                                                                                                  >  1:70498/1‑249 (MQ=255)
                                                                                                                                  aaGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTTTATACAGCTGACGCGGACTGCCCATCTTCa                                                                                                                  >  1:764435/1‑250 (MQ=255)
                                                                                                                                             tGAAGGTTTTTGGTTTTTGATAATTTGAGCTGTTCTTTTCTTATTATATATCTGTATTTGTTTTTTAATTTATTTTTTTTTTTCATTTTTTTGTGTTTAAGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACt                                                                                                        <  2:606515/249‑1 (MQ=255)
                                                                                                                                             tGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTTCTATCCTGTATACAGCTGACGCGGACGGGCAATCTTCTTACCGTCCCtt                                                                                                       >  2:601567/1‑249 (MQ=255)
                                                                                                                                                    atgggtttgttaaTAATTTGAGCTGTTCTATTTTTTTAAAATTTCTATATAGTTCTGTTAATGGATTTTATTTTTACATTTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCAtt                                                                                                 <  1:578542/237‑1 (MQ=255)
                                                                                                                                                                 tAATTTTAGCTTTTTTATTTTTTTTAAATATCTATATAGGTTTTTAAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGAAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGc                                                                                   <  1:724765/250‑1 (MQ=255)
                                                                                                                                                                   aTTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATTATCTGTATACAGCTTACCCGGACGGGCAATTTTCATACCGTCCCTGTGCATTTTCCTCCACTGGGCgg                                                                                 >  2:209314/1‑249 (MQ=255)
                                                                                                                                                                    tttGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGGTATCGCTTTTTAAGTCGCGTTTTTTATATCCTGTATACAGCTTAAGCGGACCGGCAATCTTCATACCGTCCCTGTGCATTTCGCTCCCGTGGGCGAt                                                                                >  2:423520/1‑250 (MQ=255)
                                                                                                                                                                               cTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTATAACTTAGCAAACAACCATTAAAGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCCTGATACAACTTACGCGGACGGGCCAACTTTATACCCTCACCGTTCATTTCGCTCCAGTGGGCGATCCCGCCAACgg                                                                     >  2:459637/1‑250 (MQ=255)
                                                                                                                                                                                   tCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCCACGGTACg                                                                 >  1:394452/1‑250 (MQ=255)
                                                                                                                                                                                      ttttAAATTTTTTTATAGGTTTGTTAATGTTTTTTTTTTTTACATTTTTTTGTTTTTAGGCAAAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTg                                                               <  1:95448/249‑1 (MQ=255)
                                                                                                                                                                                         tAAATATTTATATATTTCTTTTAATTTTTTTTTTTTTTATATTTTTTTGTGTTTAGACATATAAAAACCCACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCa                                                            <  2:749215/249‑1 (MQ=255)
                                                                                                                                                                                         tAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGAAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGGCGGGCAAACTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAACCAACGGTACGTGCCAt                                                           >  1:380497/1‑250 (MQ=255)
                                                                                                                                                                                           aaTATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAATGGCGGGTTAAAAAAATTAAAACTTAACAATCAACCATCAACGCCTGATATTGCTTTTTTAGACGTGTTTTTTAAATACTGTATACAGCTGTAACTGCCCGGGCATCTTTTTTCCGTCCCTATGCTTTTTGATCCCTTGGGCGAACCATCCTTATGTACGTGCCCttt                                                         >  2:295843/1‑249 (MQ=255)
                                                                                                                                                                                           aaTATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTATCACTTATAAATCAACCATTAATGCTTGTTTTAGCTTTTTAATAAGTGTTTTTTATATCCTGTTTAAAGCTGACGCGCACACGCAATGTTCTTATCTTCTCAGGGACATTCGGCCCAGTAGGCGCTCCCACCAACGGTActagccctcg                                                         >  2:845540/1‑241 (MQ=255)
                                                                                                                                                                                               tCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGTTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATTTCCTGTATACCGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGCGGGCGATCCAGCCAACGGTACGTGCCATTGCGa                                                      >  1:766787/1‑249 (MQ=255)
                                                                                                                                                                                                   tataTGTCTATTAATGTTTTTTTTTTTAATATTTTTTTTTTTTTAGAAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATg                                                 <  2:52415/250‑1 (MQ=255)
                                                                                                                                                                                                    atggttttGTTAATGGTTTTTATTTTAAAAATTTTTTGTTTTAAGAAAAAAAAAAACAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGa                                                <  1:668568/243‑1 (MQ=255)
                                                                                                                                                                                                                        tattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAACAATCAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGGAGGGCAATCTTCATATCGTCACTGTGCATTTTCCTCCAGTGGGGGATCCAACCAACGGTTCGTGCCATTGCGAAAATGAAGGTGAACATGGAAGACgg                             >  2:419024/1‑249 (MQ=255)
                                                                                                                                                                                                                         attttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAAATTAGCAATCAAACATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTTATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGTATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAAAATGGAAGACGGa                            >  2:724797/1‑249 (MQ=255)
                                                                                                                                                                                                                         attttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTAAAAATTAGCAATCAAACATTTACGCTTTATATCACTTTTAAAATCGCGTTTTTTATATCCTGTTTAAAGCTGACTCGGGCCGGCAAACTTCATACCCTCACTTTGCATTTCGCCCCAGTGGGCGTATCATGCAACGGGACGTTCCATTGCGTAAATGACGGTGAACACGGAAggagga                            >  2:552006/1‑244 (MQ=255)
                                                                                                                                                                                                                            tttACATTTTTGTGTTTTACGACATAAAAAAAAAACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATa                         <  1:743770/249‑1 (MQ=255)
                                                                                                                                                                                                                                tatttttttGTTTTTAGACATATAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCa                     <  2:349078/246‑1 (MQ=255)
                                                                                                                                                                                                                                aatttttttttttttAGGAATAAAAAAATTAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCa                     <  1:788297/246‑1 (MQ=255)
                                                                                                                                                                                                                                     ttttGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGGGAACATGGAAGACGGAATACCCATCGCtt               >  1:270877/1‑250 (MQ=255)
                                                                                                                                                                                                                                     ttttGTGTATAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCt                <  1:795593/249‑1 (MQ=255)
                                                                                                                                                                                                                                        tgtgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTAAAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGATGACGCGGACGGGCAATCTTCATTCCGTCCCTGTGCCTTTTTCTCCAGTGGGCGATCCCGCCAACGGTTCGTGCCCTTTCGAAAATGACTGTGAACATGGAAGACGGGATACCCACCGCttt              >  2:822018/1‑248 (MQ=255)
                                                                                                                                                                                                                                                 gaataaAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATAc   <  1:657931/248‑1 (MQ=255)
                                                                                                                                                                                                                                                   atatAAAAATCAACCCGCCATATGAACGGCGGGGTTAAATATTTAAAACTTAGCAATCAACCCTTTACGCCTGTTATTGCTTTTAAAGTCGCGGTTTTCCTCTCCTTTATATAGCTGACGCTGCTGGGCAATCCTCATATCCTCACTGTGCATTTTTCTCCCGTGGGCGCCCCCGCCACCGGTCCTCGCCATTTCTAAAATGACGGTGGACATGGGAGACGGAATACCCCTCTCTTTCAGGGCGATAcc  >  2:38561/1‑249 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
AAAAATATTGTTAGTTTTATTAAATGCAAAACTAAATTATTGGTATCATGAATTTGTTGTATGATGAATAAAATATAGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGATGATACC  >  NZ_CP009273/748336‑748826

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