Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9110GA28.3% 73.7 / 14.3 46intergenic (‑122/+1)fliA/fliXRNA polymerase sigma factor FliA/small RNA FliX
Reads supporting (aligned to +/- strand):  ref base G (29/4);  new base A (0/13);  total (29/17)
Fisher's exact test for biased strand distribution p-value = 2.34e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTT‑ATCCATTAC‑ACCTTCAGCGGTATAGAGTGAATTCACGATAAA‑CAGCCC‑TGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTT  >  NC_000913/2001679‑2002157
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        
cgcgGGCAGGCTAACCGGCAGGCGCAATGCTTCGTGACGCTCCAGCGGGACAAATCTGTGCCAAAGCGAGTGTTG‑ATCCATATC‑ACCTTCAGCGGTAGAAATTGAATTCACG‑TAAAGCAGCCC‑TGGGTTATATGAGTTATCGGCATGATTATCTGTGT‑CTGCAGGGTTTTTAA‑T‑CGGGGGATTAGTGGGGGAAATGAGGGGTTTTTTGGGGT‑TC‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCgg                                                                                                                                                                                                                                         <  2:487506/249‑1 (MQ=255)
  cgGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTT‑ATCCATTAC‑ACCTTCAGCGGTATAGAGTGAATTCACGATAAA‑CAGCCC‑TGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGc                                                                                                                                                                                                                                     >  1:534997/1‑251 (MQ=255)
          gaaaCCGGCAGGCGCATTGCTGCGTGACGCCCCAGCGGGACATAACGCTGCCACAGGGAGTGTTT‑ACCCTTTCC‑CCCTTCAGGGGTATAGAGTAAATTCACAAAAAA‑CGGCCC‑GGCTTTATATGAGTTATCGGCATGATTATCCGTTT‑CGGCAGGGTTTTTAA‑T‑CGGACGATTGGGGGGTAAAATG‑GGGGTTTTTTGGGGGGTT‑AAAGGTAAA‑TCCCAGGAA‑AAAAAAACCCCC‑GCCGGTGGCGGGGAAGc                                                                                                                                                                                                                             <  1:704994/249‑1 (MQ=255)
          cGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTT‑ATCCATTAC‑ACCTTCAGCGGTATAGAGTGAATTCACGATAAA‑CAGCCC‑TGCGTTATATGAGTTATCGGCATGATTATCCGATT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGTGGGTGAACTGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAAACCC‑GCCGGTGGCGGGGaaac                                                                                                                                                                                                                             >  2:574578/1‑249 (MQ=255)
              ccTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTT‑ATCCATTAC‑ACCTTCAGCGGTATAGAGTGAATTCACGATAAA‑CAGCCC‑TGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGGGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGt                                                                                                                                                                                                                         >  2:334829/1‑251 (MQ=255)
                 gCAGGCGAAATGCTTCGTGACGCACCAGCGGGAAATAACGCTGCCACAGCGAGTGTTT‑ACCCATTAC‑CCCTTCAGCGGTATAGAGTGAATTCACGAAAAA‑CAGCCC‑TGCGTTATATGAGTTACCGGCTTGATTACCCGTTT‑CTGCGGGTTTTTTAA‑T‑CGGACGATTAGTGGGTGAAATGGGGGGTTATTGGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGc                                                                                                                                                                                                                      <  1:302772/251‑1 (MQ=255)
                          atgcgtcggacctcgccgggggaaCATAACGCTGCAACAGCGGGTGTTT‑ATCCTTGTC‑CCCATCAGGGGTATAGAGGAAATCCACGATAAA‑CAGCCC‑TGCGTTATAGGAGTGATCGGCAGGATTATCGGTTG‑CT‑AGGGGTTTTTAAAT‑CGGACGATGAGTGGGTAAAAGGAGGGGTTTTTGGGGGG‑TT‑A‑AAGCAAA‑TTTCAAGCAAGAAAAATCCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAAtg                                                                                                                                                                                                            <  1:366104/228‑2 (MQ=255)
                           ttttcgttGACGCCCCAGCGGGACAAACCTCTGCCTCAGGGTGTTTTT‑ATCC‑TTACCACCTCCAGCGTTAAATGGTGAATTCAGTAAAAA‑CAGTCC‑TGCGTTATATGTGTTATCGCAATTATTATTCGTTT‑CTGC‑GGGTTTTTAAAT‑CGGACAATTAGGGGGTGAAATGGGGGTTTTTTTGGGGG‑TT‑ACAGGAAAA‑TCCCGGGAA‑AAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAAtt                                                                                                                                                                                                            <  1:432184/243‑1 (MQ=255)
                                  tGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTT‑ATCCATTAC‑ACCTTCAGCGGTATAGAGTGAATTCACGATAAA‑AAGCCC‑TGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGTGGGTGAACTGCGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGGGGCGGGGGAGCACGTTGCTGCACAAAttcactttt                                                                                                                                                                                                     >  2:614951/1‑244 (MQ=255)
                                  ggcccccccAGCGGGACACACCGATGCCAAAGGGAGTGTTA‑ATACATTCC‑CCCTTCAGGGGTAAAGAGTAAATCCACGATAAA‑CGGCCCGTG‑GTTATATGAGTTATCGGAATGATT‑TCCGTTTTCTGCAGGTTTTTAAA‑T‑GGGAGGATTAGTGGGTAAAATGGGGGGTTATTGGGGGG‑TT‑ACAGGTAAA‑TTCCAGGAA‑AAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTa                                                                                                                                                                                                    <  1:407608/244‑1 (MQ=255)
                                     cGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTT‑ATCCATTAC‑ACCTTCAGCGGTATAGAGTGAATTCACGATAAA‑CAGCCC‑TGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGt                                                                                                                                                                                                 >  1:225804/1‑251 (MQ=255)
                                               ggACATAACGCTGCCCAAGCG‑GTGTTTTATCCTTTAC‑ACCTTCAGCGTTAGAGAGTGAATCCATAATAAA‑CAGCCC‑GGCGTTATATGAGTTTTCGTAATGATTATCCGTTT‑CTGGAGGTTTTTAAA‑T‑GGGAATATTAGTGGGTGAAATGAGGGTTTTTTGGGGGG‑TT‑ACGGTTAAA‑TTCCAGGAA‑AAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATgc                                                                                                                                                                                       <  1:558893/251‑1 (MQ=255)
                                                        gctgctACAGCGAGTGTTT‑ATCCATTAC‑ACCTTCAGCGGTATAGAGTGAATTCACGATAAA‑CTGCCC‑TGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGTGGGTGAAATGAGGGGTTGTTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACACC‑GCCGGCGGCGGGGAAGAACGTTGCTG‑CCAAATTGCGCTTTATGTTCCCGGATGCGGCGTaacc                                                                                                                                                                              >  2:637476/1‑249 (MQ=255)
                                                         cTGCCACAGCGAGTGTTT‑ATCCATTAC‑ACCTTCAGCGGTATAGAGTGAATTCACGATAAA‑CAGCCC‑TGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                                                                                                                             >  1:109906/1‑251 (MQ=255)
                                                                     gtgtTT‑ATCCATTAC‑ACCTTCAGCGGTATAGAGTGAATTCACGATAAA‑CAGCCC‑TGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                                                 >  1:824821/1‑251 (MQ=255)
                                                                                              cgggataGAGTGAATTCACGATAAA‑CAGACC‑TGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGGGGGTGAAATGTGGGGTTATTTGGGGG‑TT‑AACGGTAAA‑TTCCAGGCA‑GCAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTTCTG‑ACACATTGCGCTTTATGTTGCCGGATGCGGCGTAACCGCCTTATCCGGCCTACAACACTGTGCCAATTCAAt                                                                                                                                           >  2:497238/5‑250 (MQ=255)
                                                                                               ggTATAGAGTGAATTCACGATAAA‑CAGCCC‑TGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACCAAAATGTGCAAATTCAATa                                                                                                                                          >  1:614759/1‑250 (MQ=255)
                                                                                                gTATAGAGTGAATTCACGATAAA‑CAGCCC‑TGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCCGGGTTTTTAT‑T‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGAA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAACACGTTTCTG‑ACAAATAGAGATTTCTGTTACCGGGTGCGACGTAAACGATTTATCCGGCATACAAAAATGTGCAAATTCAATaa                                                                                                                                         >  2:926984/1‑250 (MQ=255)
                                                                                                gTATAGAGTGAATTCACGATAAA‑CAGCCC‑TGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGTGGGTGACATGAGGGGTTATTTGGGGG‑TT‑ACAGGAAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGGAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                                                                                                                         >  1:928435/1‑250 (MQ=255)
                                                                                                  ataGAGTGAATTCACGATAAA‑CAGCCC‑TGCGTTATATGCGTTATCGGCAGGATTATACGTTT‑CGGCAGGGTTTTTAA‑T‑CGGACGATTAGAGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGAGGCGGGGAATCACGTTTCTG‑TCAAATTGCGCTTTATTGTGCGGGATTCGGCGAAAACGCCATATCCGGCCTACAAAAATAGGCAAATTAAATAAAtt                                                                                                                                      >  2:924067/1‑251 (MQ=255)
                                                                                                        tGAATTCACGATAAA‑CAGCCC‑TGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGTGGGTGAACTGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGATGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATTGTGCCAATTCAATAAATTGCAAtt                                                                                                                                >  2:179475/1‑251 (MQ=255)
                                                                                                        tGAATTCACGATAAA‑CAGCCC‑TGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                                                                                                                >  2:61999/1‑251 (MQ=255)
                                                                                                                ataaaaa‑CAGCCT‑TGTTTTATATGAGTTATCGGTATTTTTTCCCTTTT‑CGGCAGGGTTTTTAA‑G‑CGGACGATTAGGGGGGAAAATGGGGGGTTATTGGGGGG‑TT‑AGAGGAAAA‑TCCCAGGAA‑AAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                                                                                                                        <  1:50922/247‑1 (MQ=255)
                                                                                                                                      tatGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc                                                                                                    >  1:387579/1‑251 (MQ=255)
                                                                                                                                            ttATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGAGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGAGGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAAAAATTGCAATTCAACTTGTAGGCCTAATAAGCGCAGCGCATCAgg                                                                                              >  2:150626/1‑251 (MQ=255)
                                                                                                                                                cGGCATGATTATCCGTTT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGGAGCCCGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAAt                                                                                          >  1:399841/1‑251 (MQ=255)
                                                                                                                                                   cATGATTATCCGTTT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCA‑GCCGGTGGCGGGGAAGCAAGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTg                                                                                       >  2:860887/1‑251 (MQ=255)
                                                                                                                                                      gATTCTCCGTTT‑CTGCAGGTTTTTTAA‑T‑CGGAAGATGAGGGGGTGAAATGGGGGTTTATTTGGGGT‑TT‑AAAGGTAAA‑TTCCAGGAA‑AAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGGTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGc                                                                                     <  1:702857/250‑1 (MQ=37)
                                                                                                                                                            cctttt‑CTGCTGGTTTTTAAA‑T‑CGGATGTTTGGGGGGTGAAATGAGGGTTTATTTGGGGT‑TT‑ACAGGTAAA‑TTCCAGGAA‑AAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGcc                                                                               <  1:890502/247‑1 (MQ=37)
                                                                                                                                                                 t‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAgg                                                                         >  2:527587/1‑250 (MQ=255)
                                                                                                                                                                    tGCAGGGTTTTTAA‑T‑CGGACGAGTAGTGGGTGAAATGAGGGGGTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACACC‑GCCGGTGGCGGGGAAGCCCGTTGCTG‑ACAAATTACGCTTTATGGTGCCGGAGGCGGCGTCAACGCCTTATCCGGCCTACAAAAATGGGCAAATTAAATCAATTGCAATTCAAATTTGAGGCCTGATAAGCGCAGCGCCACAGGCCATTTGGCGttccattcggact                                                                       >  2:951184/1‑240 (MQ=255)
                                                                                                                                                                        gggTTTTTTA‑TGCGGAAGATTGGGGGGAGAAATGAGGGGTTTTTGGGGGG‑TT‑ACAGTAAAA‑TCCCAGACA‑AAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAg                                                                    <  1:877015/251‑1 (MQ=37)
                                                                                                                                                                          gTTTTTAA‑T‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGtt                                                                  >  1:889132/1‑250 (MQ=255)
                                                                                                                                                                               tttagtggtactataaggtgggaaaaTGAGGGTTAATTGGGGGG‑TA‑ATAGGAAAA‑CTCCGGACA‑AAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCGGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAAt                                                               <  1:604542/227‑1 (MQ=25)
                                                                                                                                                                                 a‑T‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACACAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGt                                                          >  2:738749/1‑251 (MQ=255)
                                                                                                                                                                                      ggACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTAc                                                       >  1:632782/1‑251 (MQ=255)
                                                                                                                                                                                             gaGTGGGAGAAATGGGGGGTTATT‑GGGGG‑TT‑AACGGTAAAATTCCAGGAA‑AAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGAtt                                                <  1:39575/250‑1 (MQ=37)
                                                                                                                                                                                                    gtgaaTGGGGGGGTATTGGGGGG‑TTCAAAGGTAAA‑TTC‑AGGCA‑AAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCt                                          <  1:574578/248‑1 (MQ=37)
                                                                                                                                                                                                     gAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCACTTTGGCGTTGCCGTCAGTCTCCGTTAATCCGGTTACAACGATTAACCCTgc                                        >  2:400579/1‑251 (MQ=255)
                                                                                                                                                                                                         tGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACCAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTgccgc                                     >  2:259001/1‑247 (MQ=255)
                                                                                                                                                                                                           gggggTTATTTGGGTG‑TA‑ACGGGAAAA‑TTCCAGGAA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCaga                                  <  1:620958/250‑1 (MQ=255)
                                                                                                                                                                                                           aGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTgcagc                                     >  1:527130/1‑248 (MQ=255)
                                                                                                                                                                                                                    ttGGGGG‑TT‑ACAGGGAAA‑TTACAGGAA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCAAGTTGCTG‑ACAAATTGAGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCACTAAATTGCAATTCAACTTGTAGGCCTAATAAGCGCAGAGCATCAGGCAATTTGGCGTTGCCGTCAGTATCAGTTAATCAGGTTACAAAGCTTAACCCTgcagcatagccagcacc                         >  2:329432/1‑240 (MQ=37)
                                                                                                                                                                                                                                cAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCCATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGGGACAGAACCTGCTGCGGt                >  2:354602/1‑250 (MQ=255)
                                                                                                                                                                                                                                  atgaaatatCCAGGAA‑AAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTAc              <  1:628277/243‑1 (MQ=255)
                                                                                                                                                                                                                                              ggCA‑AAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCttt  <  1:61999/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        
CGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTT‑ATCCATTAC‑ACCTTCAGCGGTATAGAGTGAATTCACGATAAA‑CAGCCC‑TGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTTAA‑T‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAA‑TTCCAGGCA‑GAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTG‑ACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTT  >  NC_000913/2001679‑2002157

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