Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,995,3920GA43.8% 49.5 / 25.0 48intergenic (‑122/+199)fliA/fliCRNA polymerase sigma factor FliA/H48 family flagellin FliC
Reads supporting (aligned to +/- strand):  ref base G (21/6);  new base A (0/21);  total (21/27)
Fisher's exact test for biased strand distribution p-value = 1.66e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCG‑TTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCT  >  NZ_CP009273/1995146‑1995636
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
tCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCG‑TTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGaa                                                                                                                                                                                                                                                    >  1:386696/1‑249 (MQ=255)
       ccACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCG‑TTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGAAAAAAAcc                                                                                                                                                                                                                                             <  1:125566/249‑1 (MQ=255)
               gcgGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCG‑TTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTg                                                                                                                                                                                                                                    >  1:555663/1‑250 (MQ=255)
                  ggCAGTCAAACCTGCGGGCGAAATGTTTCGGGACGCCCCAGGGGGACATACCGCTGCCACAGCGAGTGTTTATCCATTACACCTTCGGGGTTATAGAGTGAATTCACGATAAACAGCCCTGCT‑TTATATGAGTTATCGGCATGATTACCGTTTTCTGCAGGGTTTTAAATCGGACGATTAGTGGGTGAAAGGAGGGTTTTTTTGGGGGTTACAGGTAAATCCC‑AGGC‑AAAAAAAACCCCCGCCGGTGGCg                                                                                                                                                                                                                                 <  2:140572/250‑1 (MQ=255)
                                   gcgcAATGCTTCGTGACGCCCCAGCGGAACATAACGCTGCCCCAGCGAGTTTTTATCCATTACCCCTTAAGCGGTATAGAGTAAATTCACGATAAACAGCCCTGCG‑TTATATGAGTTATCGGAATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGGTTATTGGGGGGTTACAGGAAAATTCC‑AGGC‑AAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGa                                                                                                                                                                                                                <  1:425880/250‑1 (MQ=255)
                                     gcAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCG‑TTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGaa                                                                                                                                                                                                                                                    >  1:316127/1‑212 (MQ=255)
                                      cAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCG‑TTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGaaa                                                                                                                                                                                                                                                   <  2:316127/212‑1 (MQ=255)
                                                 gACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCG‑TTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTa                                                                                                                                                                                                  >  2:11150/1‑250 (MQ=255)
                                                   cGCCCCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACCCCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCT‑TTATATGATTTATCGGAATGATTATCCGTTTCTGCAGGGTTTTAAATCGGACGATTAGTGGGTAAAATGAGGGTTTATTTGGGGTTTACAGGTAAATTCC‑AGGC‑AAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATg                                                                                                                                                                                                <  1:608055/250‑1 (MQ=255)
                                                       ccAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCG‑TTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGc                                                                                                                                                                                            >  1:235754/1‑250 (MQ=255)
                                                         aGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCG‑TTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                                                                                                                          >  1:200725/1‑250 (MQ=255)
                                                               aCATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCG‑TTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTCGCCGGATgcg                                                                                                                                                                                    >  2:90218/1‑250 (MQ=255)
                                                                          ccAAAGCGAGTTTTTATCCATTACCCCTTCAGCGGTATAGAGTGAATTCACGAAAAACAGCCTTGCT‑TTATATGATTTATCGGCATGATTATCCGTTTCGGCAGGGTTTTTAATCGGACGATTAGTGGGTAAAAGGAGGGGTTTTTGGGGGGTTACAGGAAATTTCA‑AGGC‑AAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                                                                                                                           <  1:552162/248‑1 (MQ=255)
                                                                            acaGGGGGTTTTTATCCTTTACCCTTGCGGGGGTATAGGGTGAATTCACTAAAAAC‑GCCCTGGGTTTATATGATTTATGGGCATGTTTTTCCTTTTCTGGGGGGTTTTAAATCGGACGATTAGTGGGGAAAAGGGGGGTTTTTTTGGGGGTTACAGGTAAATCCC‑AGGA‑AAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                                                                                                                       <  2:561158/250‑1 (MQ=255)
                                                                            acaGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCG‑TTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                                                                                                                        >  2:214166/1‑249 (MQ=255)
                                                                                 gAGTGTTTATCCATTCCACCTTCAGCGGTAGAGAGTTATTTCACGATAAACACCCCTGCT‑TTATATGATTTATCGGCATGATTTTCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTGGGGGGTTACAGGTAAATCCC‑AGGC‑AAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                                                  <  1:427168/250‑1 (MQ=255)
                                                                                        tATCCATTCCCCCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCT‑TTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTAAAATGAGGGTTTATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                                                                                            <  1:344786/249‑1 (MQ=255)
                                                                                         aTCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCG‑TTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                                                                                                           >  2:258079/1‑249 (MQ=255)
                                                                                           ccATTACCCCTCCAGGGGTATAGAGTGAATTCACAAAAACCGGCCCTGTT‑TTATAGGAGTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAAGCGGACGATTAGTGGGTGAAAGGAGGGGTTTTTTGGGGGTTACAGGTAAATTCC‑AGGA‑AAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTCATCCGGCCTACaaaa                                                                                                                                                        <  2:72421/250‑1 (MQ=255)
                                                                                              ttccacTTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCT‑TTATATGATTTATCGGAATGTTTATCCTTTTCTGCAGGGTTTTTATTCGGACGATTAGTGGGTAAAATGAGGGTTTATTTGGGGTTTACAGTTAAATCCC‑AGCA‑AAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                                                                                                                     <  1:286693/247‑1 (MQ=255)
                                                                                                            gTATAGAGTAAATTCACGATAAACCGCCCTTGT‑TTATATGAGTTATCGGCATGTTTATCGTTTTCTGGAGGGTTTTTAATGGGACGTTTGGGGGGTAAAAGGGGGGTTTTTTTGGGGGTTAAAGGTAATTCCC‑AGCA‑AAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                                                                                                                       <  2:372427/250‑1 (MQ=37)
                                                                                                                          cACGATAAACAGCCCTGCG‑TTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGGGCCAATTCAATAAATTGCAATTCAACt                                                                                                                         >  2:416691/1‑250 (MQ=255)
                                                                                                                           agtaTAAACGGCCTGTCG‑TTATATGGTTTATGGCAATGATTGTCCGTTTGTGCGGGGTTTTAAAGGGGAAGTTGGGGGGGGGAAAGGAGGGGTTTTTTGGGGGTAAAAGAAAAATCCC‑AGGA‑AAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                                                                                                                        <  2:357835/247‑1 (MQ=25)
                                                                                                                           aCGATAAACAGCCCTGCG‑TTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTAAAATGAGGGGTTTTTGGGGGGTTACAGGAAAATCCC‑AGGC‑AAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACt                                                                                                                         <  1:534242/249‑1 (MQ=255)
                                                                                                                             aaaaaccAGCCCTGGT‑TTTTATGGGTTATCGGCGTGTTTACCCGTTTCTGCAGGTTTTTTAATGAGAGGATTAGGGGGTAAAAGGAGGGTTTATTTGGGGTTTACAGGTAATATCC‑AGGA‑AGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                                                                                                                      <  2:263105/244‑1 (MQ=37)
                                                                                                                              aaaaaCAGCCATGCT‑TTATATGATTTATCGGCATATTTATCCGTTTCTGCGGGTTTTTTATTCGGAGTATTAGTGGGAAAAATGGGGGGTTTTTTGGGGGTCAAAGGTAAATCCC‑AGGC‑AAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                                                                                                                      <  2:286530/247‑1 (MQ=37)
                                                                                                                                       ccTGCG‑TTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                                                                                                             >  1:529053/1‑249 (MQ=255)
                                                                                                                                                    gattttgcGGCAGGATTATCGTTTTCTGAGGGGTTTTTAATGGGACGTTTAGTGGGGAAAATGAGGGTTTATTTGGGGGTTACAGGTAAATCCA‑AGGA‑AAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt                                                                                                <  2:419583/243‑1 (MQ=37)
                                                                                                                                                             gtaTGTTTCTCCTTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTGAAGTGGGGGTTTTTTTGGGGGTTACAGGAAAATCCC‑AGGC‑AAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAAtt                                                                                       <  2:444441/248‑1 (MQ=37)
                                                                                                                                                              aaTGATTATCCGTTTCTGGAGGGTTTTTAGTCGGACGATGAGTGGGTAAAAGGGGGGTTTTTTTGGGGGTTA‑GGGGAAATTCCAAGGA‑AAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAAttt                                                                                      <  2:509908/249‑1 (MQ=37)
                                                                                                                                                               aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTg                                                                                     >  1:169524/1‑250 (MQ=255)
                                                                                                                                                                tGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTgg                                                                                    >  1:503809/1‑250 (MQ=255)
                                                                                                                                                                                cAGGGTTTTTAATGGGACGATTAGTGGGTAAAATGAGGGGTTATTTGGGGGTTACAGGAAAATTCC‑AGGC‑AAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGtc                                                                      <  1:576007/248‑1 (MQ=255)
                                                                                                                                                                                    tttttaAATCGGACGATTAGTGGGTGAAAGGGGGTTTTTTTTGGGGGTTAAAGGTAAATCCC‑AGGA‑AAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGtt                                                                <  2:604078/249‑1 (MQ=37)
                                                                                                                                                                                                aCGATTGGGGGGTGAAAGGGGGGTTTTTTTGGGGGTTACAGGAAAATCCC‑GGGA‑AAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACa                                                    <  2:555663/250‑1 (MQ=37)
                                                                                                                                                                                                       gggggCGAAATGAGGGTTTTTTGGGGGGTTAAAGATAAATTCC‑AGCC‑AGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCTGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTa                                             <  1:373523/248‑1 (MQ=38)
                                                                                                                                                                                                        tGGGTGAAGTGGGGGGTTTTTGGGGAGTTA‑AGGGAAAATCC‑CGGCAAGAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGAGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTaa                                            <  1:65122/250‑1 (MQ=38)
                                                                                                                                                                                                           gaaaaaggggggtTTATTTGGGGGTAACAGTAAAATCCA‑AGGC‑AAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAACCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAAccc                                         <  2:593605/241‑1 (MQ=38)
                                                                                                                                                                                                                      gttgttttGGGGGTTACAGTAAAATGCC‑GGGA‑AAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAgaga                              <  2:187746/245‑1 (MQ=38)
                                                                                                                                                                                                                        ttATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAgaga                              >  2:573200/1‑248 (MQ=255)
                                                                                                                                                                                                                         tATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACa                            <  1:129412/249‑1 (MQ=255)
                                                                                                                                                                                                                          aTTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAg                           >  2:561756/1‑249 (MQ=255)
                                                                                                                                                                                                                               ggggTTACAGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACctg                     >  2:235450/1‑250 (MQ=255)
                                                                                                                                                                                                                                  gTTACAGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACctgctg                  >  2:322061/1‑250 (MQ=255)
                                                                                                                                                                                                                                   ttACAGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAACGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCtgctgc                 >  1:190067/1‑250 (MQ=255)
                                                                                                                                                                                                                                       aGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTa             >  1:203515/1‑250 (MQ=255)
                                                                                                                                                                                                                                        ggTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTAc            >  2:402745/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   aGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCt  >  1:374093/1‑250 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
TCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCG‑TTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGC‑AGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCT  >  NZ_CP009273/1995146‑1995636

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