Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,911,779	0	C	T	26.7%	81.8 / 10.7	45	S266S (TCG→TCA)	htpX	protease HtpX
Reads supporting (aligned to +/- strand):  ref base C (8/25);  new base T (12/0);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 4.38e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCTAAATCAGTTTGGAGATAATGGTACACACGTCTCGCTGATGGCTGCGGCTATTCTGGCAGTGATTGCTGCCTGTGTCAGTGGTTTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACCGCAAAGTAGATCAGCGGGTTGCCGTTGCTCTCTTCACCTTCAT  >  NC_000913/1911532‑1912029                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                gCTAAATCAGTTTGGAGATAATGGTACACACGTCTCGCTGATGGCTGCGGCTATTCTGGCAGTGATTGCTGCCTGTGTCAGTGGTTTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGAc                                                                                                                                                                                                                                                              <  2:8250/251‑1 (MQ=255)              ggAGATAATGGTACACACGTCGCGCTGATGGCTGCGGCTATTCTGGCAGTGATTGCTGCCTGTGTCAGTGGTGTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTGTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGc                                                                                                                                                                                                                                                <  2:128055/251‑1 (MQ=255)                              acGTCTCGCTGATGGCTGCGGCTATTCTGGCAGTGATTGCTGCCTGTGTCAGTGGTTTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCAT‑GAACAAATCACTGAGCGGTTTTGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCt                                                                                                                                                                                                                               >  1:443727/1‑251 (MQ=255)                               cGTCTCGCTGATGGCTGCGGCTATTCTGGCAGTGATTGCTGCCTGTGTCAGTGGTTTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTg                                                                                                                                                                                                                              <  2:455556/251‑1 (MQ=255)                                            gCTGCGGCTATTCTGGCAGTGATTGCTGCCTGTGTCAGTGGTTTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCg                                                                                                                                                                                                                 <  2:404527/251‑1 (MQ=255)                                                      ttCTGGCAGTGATTGCTGCCTGTGTCAGTGGTTTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCAT‑GAACAAACCAATGAGGGATTTTTTCTTTACGGTTA‑TGGAGAGAGGCATC‑ATGCTGGGTGCTTTTTGTGG‑TTTATAGcg                                                                                                                                                                                                      >  2:296779/1‑250 (MQ=255)                                                            cagggATTGCTGCCTGTGTCAGTGGTTTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑tttt                                                                                                                                                                                               <  2:326471/247‑1 (MQ=255)                                                            cAGTGATTGCTGCCTGTGTCAGTGGTTTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCAT‑GAACAACACACTGAGCGGATTCGACTT‑ACCGTTAATGCAGAGGGCCATC‑ATGCTGGTTTCTTCTTGCGG‑GTCATAGCTGG‑tttt                                                                                                                                                                                               >  1:145999/1‑251 (MQ=255)                                                                gATTGCTGCCTGTGGCAGTGGTTTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCTCCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTGCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAg                                                                                                                                                                                            <  1:192711/250‑1 (MQ=255)                                                                   tgctgcCTGTGTCAGTGGTTTACGTAACACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGgcg                                                                                                                                                                                         <  2:259781/250‑1 (MQ=255)                                                                                      ttACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGAAT‑GGACAACTCCCTGCGGGGTTTTGACTT‑AACGTTTATTCAGAGAGCCATC‑ATGCTGGGTGCTTCTTGCGG‑TTCATATCTGG‑GTTTTTGGCGCTGCAGCGCGGGAATCAttt                                                                                                                                                                     >  1:222127/1‑251 (MQ=255)                                                                                        aCGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAAAAACTCCCTTAGCGATTTTGACTT‑ACCGTTAATGCAGAGAGACATC‑ATGCTGGTTTCTTCTTGCGG‑GTTATAGCTGGGTTTTTAG‑CGGTGCAGCGCGGCAAtacttttt                                                                                                                                                                   >  1:268683/1‑244 (MQ=255)                                                                                                         gATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCTT‑CACCAACTCCCTGCGCGATTTTTACTT‑TCCGTTAATTCCGCAAGCCCTC‑ATCCTGGGTTCTTCTTTCTG‑TTTATAGCTCG‑TTTTCCGGCGCTGCAGCCCGGCAATCCATTTCTCTCTaccaaccccttt                                                                                                                                                  >  2:431554/1‑246 (MQ=255)                                                                                                         gATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCAT‑GAACAACTCCCTGAGCGATTTTGACTT‑ACCTTTAATGCCGAGAGCCATC‑ATGCTGGTTGCTTCTTTCGG‑TTCATAGATGG‑TTTTCAGGCGCTGCAGCGCGGCCAGCATTTTTTCGCGACCAACCCGttt                                                                                                                                                  >  2:324184/1‑251 (MQ=255)                                                                                                             cAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGcc                                                                                                                                              <  1:85834/251‑1 (MQ=255)                                                                                                              aGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCAT‑GAAAAACTCCCTGAGCGATTTTGACTT‑ACCGTTAATGCAGAGAGCCCTC‑ATTCTGGTTGCTTCTTTCGG‑TTTATAACCGG‑TTTTCAAGCGCGGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCg                                                                                                                                             >  1:185400/1‑251 (MQ=255)                                                                                                                   aTAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCg                                                                                                                                        <  1:394578/251‑1 (MQ=255)                                                                                                                              gcgTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGACT‑GAACAACTCCCTCCGAGATTTTGCCTT‑ACCGTTAATGAAGAGAGCCATC‑ATGCTGGGTGGTTCTTTCGG‑GTCATAGCTGG‑GTTTCAGGAGCTGACGCGCGGCACTCATTTTCTCGCGcccaaccacatttaccgaaaccgcacaagcatg                                                                                                                             >  1:438916/1‑226 (MQ=255)                                                                                                                                  cGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGaa                                                                                                                          <  1:32033/250‑1 (MQ=255)                                                                                                                                    aTCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCAT‑GAAAAAACCACTGAGCGATTTTGACTT‑TCCGATAATGCAGAGAGCCACC‑CTGCTTGGTGCTTCTTTCGGGTCCAT‑GCTGG‑GTTTCCGGCGCTGCAGCGCGGCAATCACTTTCTCGCGACCCACCAGTTTTTCCGAACCCGCCTCAGCATCGaaatc                                                                                                                       >  1:101356/1‑248 (MQ=255)                                                                                                                                         gACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTTGGGCTT‑GAACAAACCACCTAGCGGTTTTTACTT‑ACCCTTAATGCCGAGAGCCATCCAT‑CTGGTTGGTTTTTTCGG‑GTTATTGCTTG‑TTTTCAAGGGCCGCAGCGCGGCAATCCTTTTTTCGCGCCCCACCCGTTTTGTCGCACCCGCCCCAACCTGGAATT‑CACGa                                                                                                                  >  2:109765/1‑250 (MQ=255)                                                                                                                                             cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGAcg                                                                                                             <  2:83933/251‑1 (MQ=255)                                                                                                                                               cTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCAT‑TAACAACTCCCTGAGCGGTTTTTACTT‑ACCGTTTATGGAGAGAACCATC‑ATGCTGGTTTCTTCTTTCGG‑TTCATTGCTGG‑GTTTTAGGCGCTGCCGCGGGGCAAACATTTTTTCGCGACCAACCCGTTTTGTCGAACCGGCATCAGCAAGGAATTTCACG‑TGAcgcg                                                                                                           >  1:340440/1‑251 (MQ=255)                                                                                                                                                tctgtatttttttcTTCAGATCACTTCCGGTATTCCCCCGTTCGCCTAGCTTCAATTCGTTGCTCCATCGGCTGGTGGGTCAT‑GAACAACTCCCTGAGCGTTTTCGACTT‑ACCGTTAATGCAGAGCGCCATC‑ATGCTGCTGGCTCCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCTGTAATCATTTTTTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATC‑CACGATGAcgcg                                                                                                           <  2:1577/237‑1 (MQ=255)                                                                                                                                                        aTTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATTGAA‑AACTCACTTAGCGATTTTGACTT‑ACCGTTAATAGAGAGAGCCATC‑ATGCTGGGTTCTTCTTTCGG‑GTCCTAGCTGG‑GTTTCCGGCGCTGCAGCGGGGCCATCATTTTCTCGCGACCAACCAGGTTTGCCGGAACGGCATCAGCATGGGATT‑CACGATGACGCgaggaacacc                                                                                                  >  1:349709/1‑244 (MQ=255)                                                                                                                                                           tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGCGGGGCAT‑GAACAAATCACTGAGAGCTTTCGACTT‑AACGTTAATGCAGAGAGCCCTC‑ATTCCGGGTGCTTCTTGCGG‑GTCATAGCTGG‑TTTTCAAGCGCTGCAGCGCGGGAATCATTTTCTCTCGACCCACCCGTTTTTTCGGACCGGCATTAACCTGGAATT‑ctataatacgcgggaaacacacgg                                                                                               >  1:475347/1‑228 (MQ=255)                                                                                                                                                           tACTTCATATTACTTCAGGTAGTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑CCCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCACTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGACGCGAGAACCACATg                                                                                                <  1:38074/250‑1 (MQ=255)                                                                                                                                                            aCTTCATATTACGTCAGTGATTCAACCGTGNGCTGAGCTTCAATTCGTGTAGCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCAGC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGACGCGAGAACCACATgg                                                                                               <  2:245304/250‑1 (MQ=255)                                                                                                                                                               tCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGACGCGAGAACCACATGGTGa                                                                                            <  2:443727/250‑1 (MQ=255)                                                                                                                                                                  tatTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑Ca                                                                                                                     <  1:384741/222‑1 (MQ=255)                                                                                                                                                                   atTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGACGCGAGAACCACATGGTGATAAt                                                                                        <  1:108488/250‑1 (MQ=255)                                                                                                                                                                                             agagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAACAACTCTCTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGTGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCCTGGAATT‑CACGATTACGCGAGAACCCCATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTc                                                             <  1:450013/251‑1 (MQ=255)                                                                                                                                                                                                   tCAATTCGTTTATCCAGCGGCGGGTGGGGCAT‑GAACCACTCACTGAGCGGTTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGGTTCTTCTTTCGG‑TTCATAGCTGG‑GTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGGCCAACCCGTTTTTCCGAACCGGCATCAGCATGGAATT‑TACGATGACGCGCGAACCAAAAGGTGATAATGCTCGCCAGAATACCAAAAAACAGTTCCAGaa                                                        >  1:418163/1‑250 (MQ=255)                                                                                                                                                                                                            ttATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCaa                                               <  1:64495/250‑1 (MQ=255)                                                                                                                                                                                                             tATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAAc                                              <  2:145999/250‑1 (MQ=255)                                                                                                                                                                                                                   gcggcgGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACCGCAAAg                                       <  1:52359/251‑1 (MQ=255)                                                                                                                                                                                                                      gcggttgtgtCAG‑GAACAACTCACCGAGCGATTTCGACTG‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTGCTTGCGG‑TTCATAACTGG‑TTGTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGTATCATCATGGAATT‑CACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATGCCAAACGCCAGTTCCAGATCCGTTGCAACCTCAAAGTAg                                    <  2:34330/243‑1 (MQ=255)                                                                                                                                                                                                                         ggtgggtCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACCGCAAAGt                                      >  1:86944/1‑246 (MQ=255)                                                                                                                                                                                                                          gtgggtCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACCGCAAAGTAGATCa                                <  1:296779/251‑1 (MQ=255)                                                                                                                                                                                                                              gtCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACCGCAAAGTAGATCAGCgg                            >  1:412137/1‑251 (MQ=255)                                                                                                                                                                                                                                 aT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACCGCAAAGTAGATCAGCGGGtt                         >  1:344894/1‑251 (MQ=255)                                                                                                                                                                                                                                 aT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACAGCAAAGTAGATCAGCGGGtt                         >  1:451283/1‑251 (MQ=255)                                                                                                                                                                                                                                     aacaacTCACTTAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATGGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTTCCGAACCGGCATCAGCATGTAATT‑CACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACCGCAAAGTAGATCAGCGGGTTGcc                      <  2:451014/251‑1 (MQ=255)                                                                                                                                                                                                                                                    gATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGACGCGAGAACCACATGGTGATAATGCTCGCCa                                                                                <  2:347287/178‑1 (MQ=255)                                                                                                                                                                                                                                                    gATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGACGCGAGAACCACATGGTGATAATGCTCGCCa                                                                                >  1:347287/1‑178 (MQ=255)                                                                                                                                                                                                                                                         cGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACCGCAAAGTAGATCAGCGGGTTGCCGTTGCTCTCTTCACCTtcat  <  1:388891/251‑1 (MQ=255)                                                                                                                                                                                                                                                         cGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTACAGAACCGTTGCAACCGCAAAGTAGATCAGCGGGTTGCCGTTGCTCTCTTCACCTtcat  >  2:358534/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                GCTAAATCAGTTTGGAGATAATGGTACACACGTCTCGCTGATGGCTGCGGCTATTCTGGCAGTGATTGCTGCCTGTGTCAGTGGTTTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAT‑GAACAACTCACTGAGCGATTTCGACTT‑ACCGTTAATGCAGAGAGCCATC‑ATGCTGGTTGCTTCTTGCGG‑TTCATAGCTGG‑TTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATT‑CACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACCGCAAAGTAGATCAGCGGGTTGCCGTTGCTCTCTTCACCTTCAT  >  NC_000913/1911532‑1912029

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