Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,634,569	0	A	T	73.5%	16.7 / 36.5	50	F106Y (TTT→TAT)	tfaQ	putative tail fiber assembly protein TfaQ
Reads supporting (aligned to +/- strand):  ref base A (7/6);  new base T (19/17);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.29e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCAACGAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATT  >  NC_000913/1634394‑1634818                                                                                                                                                                                |                                                                                                                                                                                                                                                          cagaatcGAGAGTTCTTCCTCCGTTTCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCATGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGGCGCCgg                                                                                                                                                                                <  2:408589/245‑1 (MQ=9)         agaTTTCTGCCTCCGTTGCATGTTCCAGATCTGCAGCACCCTGAAGCGGCGCAATATTCTCACTGGCTTCCAGCATCTGGTTGGTTTTTGTTTCTCCCGCCTCCCGGGTCCGGAACAGTTTTTCTGCTTCCGTAGCTGTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGGAACATTTTCCGGTAACTGACCGAGGTCAGAAATAAATCACGCGGCGCTGTAAGCCACg                                                                                                                                                                        <  2:182191/251‑1 (MQ=11)          gaTTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGt                                                                                                                                                                       >  1:14870/1‑251 (MQ=31)              tcttcCTTCGTTGCAATTTCCCGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAATTCAGAAATAAATAACGCGTCGCCGGAACCCACGTCAt                                                                                                                                                                    >  2:228665/1‑250 (MQ=9)                 tcctccGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATagac                                                                                                                                                                >  1:163872/1‑248 (MQ=11)                 tcctccGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATag                                                                                                                                                                  <  1:455807/249‑2 (MQ=11)                      cGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGAGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCAACGTCATAGACGGttt                                                                                                                                                           >  1:267264/1‑251 (MQ=11)                        ttGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCCTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGtttt                                                                                                                                                          >  1:113341/1‑250 (MQ=14)                          ccAATTTCCAGATCTTAAGCATCCTGAAGCTGCGTAATATGCTCTCTTGCTACCTGCATCAGGCTTTTTTTTTTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCGGCTGCCTTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCTTCGCCGGAAGCCACGTCATAGACGGTTTTAcc                                                                                                                                                       <  2:34579/250‑1 (MQ=11)                            aaTTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGACTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAAACAGGTAAAATTTTACGGTAACTGCCCAAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTAATagaagggtttacccc                                                                                                                                                     >  2:353212/1‑237 (MQ=11)                                 ccAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCgatg                                                                                                                                                 <  1:104148/250‑1 (MQ=11)                                         gcagcaTCCTGTAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCaa                                                                                                                                         <  2:42159/250‑1 (MQ=11)                                           agcaTCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACg                                                                                                                                       >  1:10155/1‑250 (MQ=11)                                                ccTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATCACCAGGTAAAATTTTCTGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGGAGCCACGTCATCGACGGTTTTACCCCGATGGTCTTCAACGAGATGc                                                                                                                                 >  1:216165/1‑251 (MQ=11)                                                        ggcgCAATACTCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGAGGCCACTATgc                                                                                                                         <  2:342578/251‑2 (MQ=11)                                                              atatGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCtcatc                                                                                                                   <  1:398896/251‑1 (MQ=11)                                                               tatGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGCTGCCtcatc                                                                                                                   >  1:357487/1‑250 (MQ=11)                                                                  gCTCACTGGCTACCTGAATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTg                                                                                                               >  1:118114/1‑251 (MQ=11)                                                                              ccTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGCTGCCTCATCACTGTTGAAAACAGcc                                                                                                   >  1:314809/1‑251 (MQ=11)                                                                              ccTGAATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAAAAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACAAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAACTAAATAACGCGTCGACGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTACACGAGATGCCACGATGCCtcagccctattgacaaaaacc                                                                                                   >  2:391119/1‑233 (MQ=11)                                                                                 gcctcaggcttttttttGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGTCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATTCCACGATGCCTCATCACTGGTGAAAACAGCCAc                                                                                                 <  1:85297/239‑1 (MQ=11)                                                                                 gcagaaggcttttttttGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTTTCCTTCACCCAGGCTGTGCCGTTCCCCTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTCCCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCa                                                                                                  <  2:111333/238‑1 (MQ=11)                                                                                  cATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATACTTCACCCAGGCTGTGCCGTTCAAATTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTCACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGAAAGCCACGTCATAAACGGTTTTACCCCGATGGTCTTCACCGAGATGCCACGATGCCTCATCACTGCTGaaaaaatccaca                                                                                                >  1:52774/1‑242 (MQ=11)                                                                                         ttgttTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagccagc                                                                                        >  1:445757/1‑251 (MQ=31)                                                                                              ttttGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGa                                                                                    >  1:88245/1‑250 (MQ=31)                                                                                              ttttGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCATTCACCCAGGATGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAACATTTTCCGGAAACGGACCGTGTTCAGAAATAAATAAAGCTACGCCGGAAGCCACGTCATAGCCGGTTTTACCCCGATGGTCTTCAACTAGATGCCAATATGCATCATCACTGTTGaaacccaccacacagacagactgtt                                                                                   >  2:127949/1‑229 (MQ=2)                                                                                               tttGTTTCTGCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATACCGCGTCGCCGGAAGCCACGTCATAGACGGTGTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGaa                                                                                   <  1:159915/250‑1 (MQ=11)                                                                                                  gTTTCTTCCGCCACCTGGATCCGGAACAGTTTTTCTGCGTCCGTATCCTTCACCCAGGCTGTGCCGTTCCCCTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAAGAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACTAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCCCAAAGCCAGCCGTAatat                                                                                <  1:377651/250‑1 (MQ=11)                                                                                                    ttcttcCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCt                                                                              <  1:5485/250‑1 (MQ=11)                                                                                                     tcttcCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGGGTTCAGAAATAAATAACGCGGCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTg                                                                             <  1:264778/250‑1 (MQ=11)                                                                                                      cttcCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAAAAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGATCTTCAACGAGATGCCACCATGCCTCATCCCTGTTGAAAACATCCACAACGCCAGCCGGAATACCTggg                                                                           >  2:372255/1‑250 (MQ=11)                                                                                                        tcCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCg                                                                          <  1:77178/250‑1 (MQ=11)                                                                                                         ccGCCTACCGGATCCGGAACAGTTTTTCAGCTGCCGTATCCTTAACCCCGGCTGTGCCGTTCCAATTCTGACACTCCCCTTCCGGCGATAACCAGGAAAACTTGTCCGGTAACGGACCGAGTTCAGAAATAAATAACACGCCGCCGGAAGCGACGTCATACACCGTTTTACCCCGAGAATCTTCAACGAGATGCACCGATGACCCATCACGGTTGAAAACAGCCATGAAGCCAGACGGAatagttgaagcg                                                                        >  2:163318/1‑242 (MQ=2)                                                                                                           gCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGTACCGCGTTCAGAACTAAATAACGCGTCGCCGGAAGCCGCGTCATAGACGGTTTTACCCCGATAGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACCGCCACAAAGCCCACCGGAATATCAGGCTGTGc                                                                      >  2:100673/1‑251 (MQ=11)                                                                                                           gCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGCGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGc                                                                      >  1:71651/1‑251 (MQ=11)                                                                                                                  ggATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTAACACGATGATCTTCAACGAGATGCCACGATGACTCATAACTGTTGCAAACAGACACGAAGCCAGCCGGAATATCTGGCGGTGAAATATCg                                                               >  1:79858/1‑251 (MQ=255)                                                                                                                      ccGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTAc                                                           <  1:100673/251‑1 (MQ=11)                                                                                                                        ggAACAGTTTTTCTTCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAAc                                                                                                                                        >  1:263563/1‑172 (MQ=11)                                                                                                                        ggAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCTTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTg                                                         <  2:113341/251‑1 (MQ=14)                                                                                                                              gTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGc                                                    >  1:337490/1‑250 (MQ=11)                                                                                                                                ttttCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACATCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGcagg                                                 >  1:319809/1‑248 (MQ=11)                                                                                                                                ttttCGGCTTCCGTATCCTTCACCCAGGCTGGTCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGATATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGcag                                                  <  1:292600/250‑3 (MQ=11)                                                                                                                                  ttCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGc                                                <  1:276851/250‑1 (MQ=11)                                                                                                                                    cTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCAGGAAGCCACGTCATAGACGATTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCATCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTATCGGGACAATAACGGTACTGTTTGcaagaagc                                             >  1:262843/1‑244 (MQ=11)                                                                                                                                                cctgctaccatgCTGTGCCGTTCCACTGCTGAAACTCCCCATCCGGCGATAACCAGGTAACATTTACCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGAAAGCCACGTCATAACCCGTTTTACCCCGATGACCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAACCAGCCACGAAGCCAGCAGTAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAgg                                  <  2:88149/239‑1 (MQ=11)                                                                                                                                                         aGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAatatat                        <  1:124558/251‑1 (MQ=31)                                                                                                                                                                 ccGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAACACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCAcc                >  1:388186/1‑251 (MQ=11)                                                                                                                                                                      cctcTTCTGATATTCCCCGCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATTGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCa            <  1:423664/247‑1 (MQ=11)                                                                                                                                                                      ccACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCAc           >  1:412354/1‑251 (MQ=31)                                                                                                                                                                           tCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAAt       <  1:71102/250‑1 (MQ=31)                                                                                                                                                                                aaCTCCACTTCCGGCGATAACCAGGTACCATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTAACCCGATGATCTTCAACGAGATGCAACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGAAATATAGGTACTGTTTGCGGGCAGACCAGTATGAGGCGGGATATATGCAACACCTTCACCAATAAAtt  >  1:76295/1‑250 (MQ=11)                                                                                                                                                                                |                                                                                                                                                                                                                                                          CAGCAACGAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATT  >  NC_000913/1634394‑1634818

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