Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,995,3990AC51.1% 52.5 / 74.6 45intergenic (‑129/+192)fliA/fliCRNA polymerase sigma factor FliA/H48 family flagellin FliC
Reads supporting (aligned to +/- strand):  ref base A (19/3);  new base C (0/23);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 6.32e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-08
Rejected as consensus: Frequency below/above cutoff threshold.

CCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACG  >  NZ_CP009273/1995153‑1995648
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
ccACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcc                                                                                                                                                                                                                                                         <  1:125566/249‑1 (MQ=255)
        gcgGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTg                                                                                                                                                                                                                                                >  1:555663/1‑250 (MQ=255)
           ggCAGTCAAACCTGCGGGCGAAATGTTTCGGGACGCCCCAGGGGGACATACCGCTGCCACAGCGAGTGTTTATCCATTACACCTTCGGGGTTATAGAGTGAATTCACGATAAACAGCCCTGCTTTATATGAGTTATCGGCATGATTACCGTTTTCTGCAGGGTTTTAAATCGGACGATTAGTGGGTGAAAGGAGGGTTTTTTTGGGGGTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCg                                                                                                                                                                                                                                             <  2:140572/250‑1 (MQ=255)
                            gcgcAATGCTTCGTGACGCCCCAGCGGAACATAACGCTGCCCCAGCGAGTTTTTATCCATTACCCCTTAAGCGGTATAGAGTAAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGAATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGGTTATTGGGGGGTTACAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGa                                                                                                                                                                                                                            <  1:425880/250‑1 (MQ=255)
                                          gACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTa                                                                                                                                                                                                              >  2:11150/1‑250 (MQ=255)
                                            cGCCCCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACCCCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCTTTATATGATTTATCGGAATGATTATCCGTTTCTGCAGGGTTTTAAATCGGACGATTAGTGGGTAAAATGAGGGTTTATTTGGGGTTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATg                                                                                                                                                                                                            <  1:608055/250‑1 (MQ=255)
                                                ccAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGc                                                                                                                                                                                                        >  1:235754/1‑250 (MQ=255)
                                                  aGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                                                                                                                                      >  1:200725/1‑250 (MQ=255)
                                                        aCATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTCGCCGGATgcg                                                                                                                                                                                                >  2:90218/1‑250 (MQ=255)
                                                                   ccAAAGCGAGTTTTTATCCATTACCCCTTCAGCGGTATAGAGTGAATTCACGAAAAACAGCCTTGCTTTATATGATTTATCGGCATGATTATCCGTTTCGGCAGGGTTTTTAATCGGACGATTAGTGGGTAAAAGGAGGGGTTTTTGGGGGGTTACAGGAAATTTCAAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                                                                                                                                       <  1:552162/248‑1 (MQ=255)
                                                                     acaGGGGGTTTTTATCCTTTACCCTTGCGGGGGTATAGGGTGAATTCACTAAAAACGCCCTGGGTTTATATGATTTATGGGCATGTTTTTCCTTTTCTGGGGGGTTTTAAATCGGACGATTAGTGGGGAAAAGGGGGGTTTTTTTGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                                                                                                                                   <  2:561158/250‑1 (MQ=255)
                                                                     acaGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                                                                                                                                    >  2:214166/1‑249 (MQ=255)
                                                                          gAGTGTTTATCCATTCCACCTTCAGCGGTAGAGAGTTATTTCACGATAAACACCCCTGCTTTATATGATTTATCGGCATGATTTTCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTGGGGGGTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                                                              <  1:427168/250‑1 (MQ=255)
                                                                                 tATCCATTCCCCCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCTTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTAAAATGAGGGTTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                                                                                                        <  1:344786/249‑1 (MQ=255)
                                                                                  aTCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                                                                                                                       >  2:258079/1‑249 (MQ=255)
                                                                                    ccATTACCCCTCCAGGGGTATAGAGTGAATTCACAAAAACCGGCCCTGTTTTATAGGAGTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAAGCGGACGATTAGTGGGTGAAAGGAGGGGTTTTTTGGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTCATCCGGCCTACaaaa                                                                                                                                                                    <  2:72421/250‑1 (MQ=255)
                                                                                       ttccacTTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCTTTATATGATTTATCGGAATGTTTATCCTTTTCTGCAGGGTTTTTATTCGGACGATTAGTGGGTAAAATGAGGGTTTATTTGGGGTTTACAGTTAAATCCCAGCAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                                                                                                                                 <  1:286693/247‑1 (MQ=255)
                                                                                                     gTATAGAGTAAATTCACGATAAACCGCCCTTGTTTATATGAGTTATCGGCATGTTTATCGTTTTCTGGAGGGTTTTTAATGGGACGTTTGGGGGGTAAAAGGGGGGTTTTTTTGGGGGTTAAAGGTAATTCCCAGCAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                                                                                                                                   <  2:372427/250‑1 (MQ=37)
                                                                                                                   cACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGGGCCAATTCAATAAATTGCAATTCAACt                                                                                                                                     >  2:416691/1‑250 (MQ=255)
                                                                                                                    agtaTAAACGGCCTGTCGTTATATGGTTTATGGCAATGATTGTCCGTTTGTGCGGGGTTTTAAAGGGGAAGTTGGGGGGGGGAAAGGAGGGGTTTTTTGGGGGTAAAAGAAAAATCCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                                                                                                                                    <  2:357835/247‑1 (MQ=25)
                                                                                                                    aCGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTAAAATGAGGGGTTTTTGGGGGGTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACt                                                                                                                                     <  1:534242/249‑1 (MQ=255)
                                                                                                                      aaaaaccAGCCCTGGTTTTTATGGGTTATCGGCGTGTTTACCCGTTTCTGCAGGTTTTTTAATGAGAGGATTAGGGGGTAAAAGGAGGGTTTATTTGGGGTTTACAGGTAATATCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                                                                                                                                  <  2:263105/244‑1 (MQ=37)
                                                                                                                       aaaaaCAGCCATGCTTTATATGATTTATCGGCATATTTATCCGTTTCTGCGGGTTTTTTATTCGGAGTATTAGTGGGAAAAATGGGGGGTTTTTTGGGGGTCAAAGGTAAATCCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                                                                                                                                  <  2:286530/247‑1 (MQ=37)
                                                                                                                                ccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                                                                                                                         >  1:529053/1‑249 (MQ=255)
                                                                                                                                             gattttgcGGCAGGATTATCGTTTTCTGAGGGGTTTTTAATGGGACGTTTAGTGGGGAAAATGAGGGTTTATTTGGGGGTTACAGGTAAATCCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt                                                                                                            <  2:419583/243‑1 (MQ=37)
                                                                                                                                                      gtaTGTTTCTCCTTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTGAAGTGGGGGTTTTTTTGGGGGTTACAGGAAAATCCCAGGCAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAAtt                                                                                                   <  2:444441/248‑1 (MQ=37)
                                                                                                                                                       aaTGATTATCCGTTTCTGGAGGGTTTTTAGTCGGACGATGAGTGGGTAAAAGGGGGGTTTTTTTGGGGGTTAGGGGAAATTCCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAAttt                                                                                                  <  2:509908/249‑1 (MQ=37)
                                                                                                                                                        aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTg                                                                                                 >  1:169524/1‑250 (MQ=255)
                                                                                                                                                         tGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTgg                                                                                                >  1:503809/1‑250 (MQ=255)
                                                                                                                                                                         cAGGGTTTTTAATGGGACGATTAGTGGGTAAAATGAGGGGTTATTTGGGGGTTACAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGtc                                                                                  <  1:576007/248‑1 (MQ=255)
                                                                                                                                                                             tttttaAATCGGACGATTAGTGGGTGAAAGGGGGTTTTTTTTGGGGGTTAAAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGtt                                                                            <  2:604078/249‑1 (MQ=37)
                                                                                                                                                                                         aCGATTGGGGGGTGAAAGGGGGGTTTTTTTGGGGGTTACAGGAAAATCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACa                                                                <  2:555663/250‑1 (MQ=37)
                                                                                                                                                                                                gggggCGAAATGAGGGTTTTTTGGGGGGTTAAAGATAAATTCCAGCCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCTGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTa                                                         <  1:373523/248‑1 (MQ=38)
                                                                                                                                                                                                 tGGGTGAAGTGGGGGGTTTTTGGGGAGTTAAGGGAAAATCCCGGCAAGAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGAGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTaa                                                        <  1:65122/250‑1 (MQ=38)
                                                                                                                                                                                                    gaaaaaggggggtTTATTTGGGGGTAACAGTAAAATCCAAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAACCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAAccc                                                     <  2:593605/241‑1 (MQ=38)
                                                                                                                                                                                                               gttgttttGGGGGTTACAGTAAAATGCCGGGAAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAgaga                                          <  2:187746/245‑1 (MQ=38)
                                                                                                                                                                                                                 ttATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAgaga                                          >  2:573200/1‑248 (MQ=255)
                                                                                                                                                                                                                  tATTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACa                                        <  1:129412/249‑1 (MQ=255)
                                                                                                                                                                                                                   aTTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAg                                       >  2:561756/1‑249 (MQ=255)
                                                                                                                                                                                                                        ggggTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACctg                                 >  2:235450/1‑250 (MQ=255)
                                                                                                                                                                                                                           gTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACctgctg                              >  2:322061/1‑250 (MQ=255)
                                                                                                                                                                                                                            ttACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAACGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCtgctgc                             >  1:190067/1‑250 (MQ=255)
                                                                                                                                                                                                                                aGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTa                         >  1:203515/1‑250 (MQ=255)
                                                                                                                                                                                                                                 ggTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTAc                        >  2:402745/1‑250 (MQ=255)
                                                                                                                                                                                                                                            aGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCt              >  1:374093/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      aaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAac     <  1:529197/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       aaaCCGCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGATTTATGTTGCCGGATGCGGCGGAAACGCCTTATCCGGCCCACAAAAATGTGCAAATTCAATAAATTGCCAGTCCACTTGTAGGCCTGGTAAACGCAGGGCATCAAGCAATTTGGCCTTGCCGTCAGCCTCAGTTAATCAGGTTACAACGATTAACCCGCCAGCAGACACAGAACCTGCTGCGGAACCTGGTTAGCTTTTGCCAaca    >  2:409703/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        aaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAacac   >  1:512455/1‑250 (MQ=255)
                                                                                                                                                                                                                                                         aCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACg  >  1:277977/1‑250 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
CCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACG  >  NZ_CP009273/1995153‑1995648

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