Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,634,676	0	A	G	57.8%	‑5.4 / 44.4	45	D70D (GAT→GAC)	tfaQ	putative tail fiber assembly protein TfaQ
Reads supporting (aligned to +/- strand):  ref base A (12/7);  new base G (11/15);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTCTGAATGCCATT  >  NC_000913/1634432‑1634886                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                    cTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTc                                                                                                                                                                                                                >  2:55608/1‑250 (MQ=31)  tGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCa                                                                                                                                                                                                               >  2:3470/1‑250 (MQ=31)   gcagcaTCCTGAAGCGGCGCAATGTGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAAc                                                                                                                                                                                                             <  2:194074/251‑1 (MQ=11)   acagcaTCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCtt                                                                                                                                                                                                                 <  1:312880/246‑1 (MQ=11)    cagcagCCTGGAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTGTCTTCCGCCTCCCGGATCCCGAACAGTTTTTCTGCTTCCGTATCCTTCACCCATTCTGTGCCCTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGTCCGAGTTCAGAAATAAATAGCGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATTATCTTCAACg                                                                                                                                                                                                            <  2:363742/251‑1 (MQ=11)             gAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCGTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGcc                                                                                                                                                                                                     <  1:194458/249‑1 (MQ=11)                       aaTATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCtcat                                                                                                                                                                                         <  1:377911/251‑1 (MQ=11)                         gatGCTCGCTGGCTACCTGCATCAGGCTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGGTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAAAGAGATGCCACGATGCCtcatca                                                                                                                                                                                       <  2:225059/250‑1 (MQ=11)                            gCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTg                                                                                                                                                                                    <  1:136754/251‑1 (MQ=11)                                       aCCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAg                                                                                                                                                                         <  1:376778/250‑1 (MQ=31)                                        ccTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACACCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGcc                                                                                                                                                                       >  2:64250/1‑251 (MQ=11)                                              tCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATGCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGAAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGACGCCGGAAGCCACGTGATAAAACGTTTTACCCCGATGATCTTCTACGAGATGCCACGCTGACTCATCACGGTTGA‑AAACAGCCACGa                                                                                                                                                                   >  2:276555/1‑249 (MQ=11)                                                     gttTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAGCGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAagccagcc                                                                                                                                                           <  2:255161/250‑1 (MQ=14)                                                      tttgttGTTTCCTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCTAGTTCAGAAATAAAGAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCGCGAGGACTCGGTACTTTTGA‑AAACATCCATGAAGCCAGCCg                                                                                                                                                          <  2:151796/246‑1 (MQ=11)                                                       tttttGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCgg                                                                                                                                                         >  1:93356/1‑250 (MQ=11)                                                                      ccTCCCGGATCCGGAACAGTTTTTCTGCTCCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCa                                                                                                                                         <  1:185702/251‑1 (MQ=11)                                                                             gATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCgg                                                                                                                                  <  1:3470/251‑1 (MQ=31)                                                                              aTCCGGCACAGATTTTCTGCTTCAGTGTCCTTCACCCAGGATGTGCCGTTCCTCTTCTGAAACTCCCCTTCCGGCGATACCCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGCAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACTATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGt                                                                                                                                 >  2:233292/1‑251 (MQ=11)                                                                                    aaccgTTTTTCTTCTTCCTTATCGTCCACCCGGGCTGTGCCTTTCCACTTCTGATATTCCCCCCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGAACCGAGTTCAAAATTAAATACCCCGTCGCCGAAAGCCACGTCATAGAC‑GGTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGttt                                                                                                                          <  2:365508/247‑1 (MQ=11)                                                                                      cAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTg                                                                                                                         <  1:205293/251‑1 (MQ=11)                                                                                                        aTCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATg                                                                                                         >  2:108261/1‑249 (MQ=11)                                                                                                                    ggCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATg                                                                                           >  2:73111/1‑251 (MQ=11)                                                                                                                       tgtgCCGTTCCACTTCTCAAACTCCCCTTCCGGCGATAACCCGGAAACATTTTCCGGTAACGGACCGAGTTCAGAAATCAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCAGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGTTGCAATATCGGTCCTGCTTGCTGGCAGACCTGTATGAGGCGGAATATATGCAt                                                                                        >  1:351066/1‑251 (MQ=11)                                                                                                                                ccACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCa                                                                                <  2:71250/250‑1 (MQ=31)                                                                                                                                     tctgatattCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGTCGGAATATATTCGTCGCCTTCACCAAt                                                                           <  1:82552/242‑1 (MQ=11)                                                                                                                                     tCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAATCCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAAt                                                                           >  1:373253/1‑250 (MQ=14)                                                                                                                                          aaCTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGCAATAAATAACGCGTCGACGGAAGCCACGTCCTAAACCGTTTTAACCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGCGCAATATCGGTACTGTTTGCTGGCAGAACTGTATTCGGCGGAATATATGCCTCACCTTCACCAATAAATTc                                                                     >  2:146645/1‑251 (MQ=11)                                                                                                                                             tacccTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAATCCACGTCATAGACGGTTTTACACCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTGTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCAt                                                                   <  2:76618/248‑1 (MQ=11)                                                                                                                                              ccccTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCGCTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGATGCGGAATAGATGCATCACCTGCACCAATAAATTCAt                                                                   <  2:233735/249‑1 (MQ=9)                                                                                                                                                cctcccGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGt                                                               >  2:70946/5‑251 (MQ=11)                                                                                                                                                cctcccGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGt                                                               >  2:201437/5‑251 (MQ=11)                                                                                                                                                 ctcccGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGGATGAGGGGGAATATATGCGTCACCTGCACCAATAAATTCATTAGtt                                                              >  2:294115/4‑251 (MQ=11)                                                                                                                                                      ggCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGc                                                         >  2:64722/1‑251 (MQ=31)                                                                                                                                                       gCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGTGACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGcc                                                        <  2:240020/251‑1 (MQ=11)                                                                                                                                                                aGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGAtt                                                >  2:196474/1‑250 (MQ=31)                                                                                                                                                                 ggTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATtat                                              <  2:117849/251‑1 (MQ=11)                                                                                                                                                                      aaTTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTCTGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAAtt                                         <  2:364233/250‑1 (MQ=11)                                                                                                                                                                              ggTAACGGACCGAGTTCAGAAATAAATAACGAGTCGCCGGAAGCAACGACATAGACGGTATTACCCCGATGGTCTTCCACGAGCTGCCCAGCAGCCTCATCACTGTTGACAAAC‑GCCACAAAGCCCGCAGGAAAATCTGGCGGTGCAATATAGGTACTGTTTGCAGGCAGACCGGTCTGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGt                                 >  1:221240/1‑251 (MQ=2)                                                                                                                                                                                tAACGGAACGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCCTTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTcc                               >  2:217669/1‑251 (MQ=11)                                                                                                                                                                                   cGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCAACGTCATAAACAGTTTTACCCCGATGATCTTCAACGAGATGCAACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTCCTGTTTGCTGGAAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCCTTAGTTCCAGCCAGCAGATTATAAATTTTTATGGTCCGTg                            >  1:59684/1‑251 (MQ=11)                                                                                                                                                                                          aGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTc                     <  2:320875/251‑1 (MQ=11)                                                                                                                                                                                                gAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCAt               >  1:22930/1‑251 (MQ=11)                                                                                                                                                                                                 aaataaatAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCAtt              >  1:79928/1‑251 (MQ=11)                                                                                                                                                                                                  aataaataACGCGTCGCCGGAAGCCACGTCATACACCGTTTTACACCGATGATCTTCAACGAGATGCCACGATAACTCATCCCTGTTGA‑AAACAGCCACGAAGCAAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTc             >  1:210167/1‑251 (MQ=9)                                                                                                                                                                                                              gTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTATTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTCTGAATGCCAtt  >  2:256949/1‑250 (MQ=14)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                    CTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTCTGAATGCCATT  >  NC_000913/1634432‑1634886

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