Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,161,1480GA25.0% 128.9 / 15.3 60E28K (GAA→AAA) fabRDNA‑binding transcriptional repressor FabR
Reads supporting (aligned to +/- strand):  ref base G (39/6);  new base A (0/15);  total (39/21)
Fisher's exact test for biased strand distribution p-value = 1.02e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACC‑GGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGG  >  NC_000913/4160906‑4161397
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         
gatattaagcagaatttcccTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACC‑GGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAcc                                                                                                                                                                                                                                                   <  1:39905/234‑1 (MQ=255)
 cTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACC‑GGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAccc                                                                                                                                                                                                                                                  >  2:116519/1‑251 (MQ=255)
     ttCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACC‑GGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgcc                                                                                                                                                                                                                                               <  1:182584/250‑1 (MQ=255)
       cTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACC‑GGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGtt                                                                                                                                                                                                                                            >  2:388614/1‑251 (MQ=255)
                    tGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACC‑GGCCAAAGAATTGCAGCAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGc                                                                                                                                                                                                                               >  1:340334/1‑251 (MQ=255)
                              ttGGATTTTTATCGACGATACTCTCCGTTTAAGCGGAAGGTTTCCGCTGTACGTAAAAGAACC‑GGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCAAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCAGCCGTTCGCTGGTGGAAGCCGCaatt                                                                                                                                                                                                                        >  2:317245/1‑245 (MQ=255)
                                    ttttATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACC‑GGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAtt                                                                                                                                                                                                                >  1:266960/1‑250 (MQ=255)
                                        acggaggttACTCTCCTTTTAAGCGGCGGTTTTCCGCTGT‑CGTAAAAAAACCCGACCAAAGAATGGAAGTAATTATTTTTTATTGCGTTCCCTTTCTTTCACAAACCTGGAC‑CAACCCGGTATTTTCTTTCTCTGG‑AATAGGCCCAGC‑AGTACTTTTGGCAAGGATTCAGACATGGTGATGGGCGTAAGACCACAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGt                                                                                                                                                                                                            <  2:305370/242‑1 (MQ=255)
                                           gacgaTACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACC‑GGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTg                                                                                                                                                                                                        >  1:573033/1‑251 (MQ=255)
                                                  ctctCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACC‑GGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGc                                                                                                                                                                                                   >  2:498622/1‑249 (MQ=255)
                                                             aGCGGCAGGTTTCCGCGGTACGAAAAAGACCC‑GGCCAAAGAATTGCAGAAAATATTTTTTATTGCTTTACCGTTCTTTCACAATCCTGGAG‑CAACCCAGTATGTTCTTTCTCTGG‑TAAAGTGCAAGC‑AGTACTTTGGGAAAGGATTCAGACATCGTGATGGGGGAAAGAGCGCAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGcc                                                                                                                                                                                      <  1:61971/251‑1 (MQ=255)
                                                                   ggttttCCTCTGTCCGTAAAAGACCC‑GGCAAAAGAATTGCAGTAAATATTTTTTATTGCTTTCCCGTTCATTCACAATCCTGAAG‑CAACCCAGTATTTTCATTC‑CTGGGTATAGTGCCAGC‑ATTATTTTTGGCAAGGATTCAGACATCGTGATGGGCTTAAGAGCGCAACAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGtt                                                                                                                                                                                <  1:474163/248‑1 (MQ=255)
                                                                    ggTTTCCGCTGTACGTAAAAGAACC‑GGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGttt                                                                                                                                                                               >  1:548199/1‑251 (MQ=255)
                                                                        tCCGCTGTACGTAAAAGAACC‑GCCAAAAGAATTGAAGAAAATATGTTTTATTGCTTTACCTTTCATTCAAAATACTGGAG‑CAATCCAGTATGTTCATTCTCGGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGAAAGAGCGCAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCg                                                                                                                                                                            <  1:97652/250‑1 (MQ=255)
                                                                             tGTACGTAAAAGAACC‑GGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGt                                                                                                                                                                      >  1:561126/1‑251 (MQ=255)
                                                                              gTACGTAAAAGAGCC‑GTCCAAATAATTGCAGTAAATATGTTTGATTGCTTTCCCGTTCATTCACAATGCTGTAG‑CAATCCAGTATGTTCTTTCCCTGG‑TATAGTGC‑GGCGAGTTCTTTTGGTAGGAATTCAGACACCGGGATGGGCGTAAGAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGt                                                                                                                                                                      <  1:471376/250‑1 (MQ=255)
                                                                                       aGAACC‑GGCAAAAAAATGGCATTAAATAGGTTTTATTGCGTTACCGTTCATTCACAATACTGGAG‑CAACCCAGTATGTTCATTCTCGGG‑TATATTGCCAGC‑AGTTCTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGGGCACAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGaa                                                                                                                                                            <  1:446950/251‑1 (MQ=255)
                                                                                              ggCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGc                                                                                                                                                      >  1:215125/1‑251 (MQ=255)
                                                                                               cccaaaaaatttaaatAAATTTGTTTTATTTCGTTACCTTTCATTCAAAAGA‑TTGAGACAATCCAGTATTTTCATTCTTGGG‑AAAAGGGCCAGC‑ATTTCTTTGGGCAGGGATCCAGACATCGTGAGGGGGGAAAGGGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCa                                                                                                                                                     <  2:573033/236‑1 (MQ=255)
                                                                                                ccaaaaaatgtcAGTAAATATTTTTTATTGCTTTACCTTGCATTCACAATACTGGAG‑CAACCCAGTATGTTCATTCTCTGG‑TAGAGGCCCAGC‑TGTATTTTTGGAAGGGATTCAGACATCGTGAGGGGGGAAAGGGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCAt                                                                                                                                                    <  2:4054/240‑1 (MQ=255)
                                                                                                  aaaGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTg                                                                                                                                                  >  1:135218/1‑251 (MQ=255)
                                                                                                    aGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCt                                                                                                                                                >  2:369518/1‑251 (MQ=255)
                                                                                                     aaaTTGAAGAAAATTTGTTTTATTGGGTTCCCGTTCATTAACAATCCGGGAG‑AAATCCAGTATTTTCTTTCTCGGG‑TATGGTGCCAGC‑ATTACTTTTGGAAAGAATTAAGACATCGTGATGGGCGAAAGAGCGCAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTc                                                                                                                                               <  1:273590/250‑1 (MQ=255)
                                                                                                       aTTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTccc                                                                                                                                             >  1:54785/1‑251 (MQ=255)
                                                                                                             gTAAATTTTTTTTATTGTGTTCCCGTTCATTCACAATACTGGAG‑CAACCCAGTATTTTCATTCTCTGT‑TATGGGGCCAGC‑ATTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACctct                                                                                                                                       <  1:508066/251‑1 (MQ=255)
                                                                                                              aaaaTATGTTTTTTTGCGTTACCTTTAATTCACAATACGGGAA‑CAATCAAGTATGTTCATTCTCGGG‑TATAGTGCCAGA‑AGTACTTTTGGCAAGGATACAGAGATCGTGAGGGGCGTAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCtt                                                                                                                                      <  1:116519/250‑1 (MQ=255)
                                                                                                                          tttgcttttccttacattcaCAATGCGGGAG‑CAATCCAGTTTTTTCATTCTCGGG‑TATAGGGCCAGC‑AATACTTTTGGCAAGATTTCGAACTTTGTGATGGGCATAAGGGCCCAAAAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCAt                                                                                                                           <  2:548199/236‑1 (MQ=255)
                                                                                                                                ttACCGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCcgcg                                                                                                                    >  2:102479/1‑251 (MQ=255)
                                                                                                                                   ccGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGAGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACg                                                                                                                 >  1:334041/1‑251 (MQ=255)
                                                                                                                                    cGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGcc                                                                                                                                                                                                                    >  1:498971/1‑151 (MQ=255)
                                                                                                                                    cGTTCATTCACAATACTGGAG‑AAATCAAGTATTTTCATTCTCGGT‑TATAGTGCCAGC‑AGTCCTTTTGGCAAGGATCCAGACACCGTGAGGGGCGAAAGACCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGcc                                                                                                                                                                                                                    <  2:498971/151‑1 (MQ=255)
                                                                                                                                     gTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGt                                                                                                                >  2:335766/1‑250 (MQ=255)
                                                                                                                                          ttcaCAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACg                                                                                                           >  2:124901/1‑250 (MQ=255)
                                                                                                                                                 tACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGCTGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGAGGGCATTGGTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGt                                                                                                   >  2:448168/1‑251 (MQ=255)
                                                                                                                                                           aaTCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATg                                                                                          >  1:416705/1‑251 (MQ=255)
                                                                                                                                                           aaTCCAGTATGTTCATTCTCTGG‑CATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCAACTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGAccctg                                                                                          >  2:463149/1‑248 (MQ=255)
                                                                                                                                                             tCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGt                                                                                        >  1:81742/1‑251 (MQ=255)
                                                                                                                                                              ccAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGtt                                                                                       >  1:2966/1‑251 (MQ=255)
                                                                                                                                                                   aTGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtgatg                                                                                   >  1:437359/1‑250 (MQ=255)
                                                                                                                                                                            ctctGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGttt                                                                         >  1:136823/1‑251 (MQ=255)
                                                                                                                                                                              ctGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTa                                                                        >  2:447321/1‑250 (MQ=255)
                                                                                                                                                                                        tGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGcc                                                              >  2:437723/1‑251 (MQ=255)
                                                                                                                                                                                        ggcaagc‑agTATTTTTGGCAAGGATACAGACATCGGGGGGGGCGTAAGAGCCCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGcc                                                              <  2:35752/247‑1 (MQ=255)
                                                                                                                                                                                        gcccAGT‑AGTATTTGTGGCAAGATTTCAGACTTGGTGTTGGGCGTAAGAGCGCAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGcc                                                              <  2:416705/249‑1 (MQ=255)
                                                                                                                                                                                         cccAGC‑ATTACTTTTGGCAAGGATTCAAACATCGTGATGGGCGTAAGAGCACAACAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGcc                                                              <  1:124901/249‑1 (MQ=255)
                                                                                                                                                                                           cagc‑agTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGACGAGAGCGGTTTAATGCTACGCCAAc                                                           >  1:267210/1‑251 (MQ=255)
                                                                                                                                                                                                agTACTTTTGGCAAGGATTCAGACATCGTGAGGGGCGTAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCAt                                                       <  1:369518/251‑1 (MQ=255)
                                                                                                                                                                                                     ttttGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGGTTAATGCTACGCCAACTCATGCGcc                                                  >  1:318886/1‑251 (MQ=255)
                                                                                                                                                                                                      tttGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGCGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGGACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGcc                                                  >  1:379498/1‑250 (MQ=255)
                                                                                                                                                                                                            aaGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACACGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGAGGTTTAATGCTACGCCAACTCATGCGCCAGGcgcg                                           >  2:200679/1‑251 (MQ=255)
                                                                                                                                                                                                                  tCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCCTGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCg                                     >  2:567924/1‑251 (MQ=255)
                                                                                                                                                                                                                   cAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCt                                                                                                 >  1:405820/1‑190 (MQ=255)
                                                                                                                                                                                                                                    gTAAGAGCGCAACACAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCAAGACTGCGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGAATTTCCGCGACGGAGACGAACTGGGTCTGACTATGGTTGAGGAGAGCGGTGTCATGCTACGCCAACGCATGCGCCAGGCTCGTCAGAGATTCCCCAAAGGCgggcg                   >  2:487547/1‑249 (MQ=255)
                                                                                                                                                                                                                                      aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAgtg                 >  1:81536/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ggcgcAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGa              <  2:108708/250‑1 (MQ=255)
                                                                                                                                                                                                                                            gcAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATcc           >  1:452474/1‑251 (MQ=255)
                                                                                                                                                                                                                                             caacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGAGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGCGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGTGGGAGTGTGATCCg          >  2:5912/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    aaGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATgc                                                     <  1:12920/201‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATgc                                                     >  2:12920/1‑201 (MQ=255)
                                                                                                                                                                                                                                                     aGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCgg  >  1:88656/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         
TCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACC‑GGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAG‑CAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGC‑AGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGG  >  NC_000913/4160906‑4161397

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: ‑