Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,474,290	T→C	K88R (AAA→AGA)	rpsL ←	30S ribosomal subunit protein S12

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,474,290	0	T	C	95.7%	140.1 / ‑1.4	47	K88R (AAA→AGA)	rpsL	30S ribosomal subunit protein S12
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (22/23);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.28e-01

AATTTTACGCTGACCAATGACGCGACGACGTGGCATGGAAATACTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTG  >  NC_000913/3474046‑3474526                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             aaTTTTACGCTGACCAATGACGCGACGACGTGGCATGGAAATACTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGGGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAacac                                                                                                                                                                                                                                        <  1:82982/251‑1 (MQ=255)  aTTTTACGCTGACCAATGACGCGACGACGTGGCATGGAAATACTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGCCATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACg                                                                                                                                                                                                                                       >  1:16494/1‑251 (MQ=255)          cTGACCAATGACGCGACGACGTGGCATGGAAATACTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGCCCGCCAc                                                                                                                                                                                                                               >  1:37816/1‑251 (MQ=255)              ccAATGACGCGACGACGTGGCATGGAAATACTCCTTTGTTAATTCAGTATTGTCCAAATCTGTCCGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCTTTATTTAACGGAGAACCATTAAGCCTTAGGCCTCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACggat                                                                                                                                                                                                                           <  2:276690/251‑1 (MQ=255)               cAATGACGCGATGACGTGGCATGGAAATACTCCTTTGTAAATTCAGGATTGTCAAAAACTCTAGGAGTTTAGTTTGACATTTAAGTTAAAACGGTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATATTTGGAACGAGCCTGCTTACGGTCTTTAACGCGGGAGCAGCCAAGCGCACCTCGTACGGGGTGGTAAGGAACACCGGGGAGGGCTCTAACACGACCGCCACggatc                                                                                                                                                                                                                          <  2:144281/251‑1 (MQ=255)                   gACGCGACGACGTGGCATGGAAATACTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACggatcagga                                                                                                                                                                                                                      <  1:206099/251‑1 (MQ=255)                      gcgACGACGTGGCATGGAAATACTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACggatcaggatc                                                                                                                                                                                                                    >  1:365762/1‑250 (MQ=255)                         acgacgTGGCATGGAAATACTCCGTTGTTACTTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACAGCCACGGATCAGGATCACgg                                                                                                                                                                                                                >  2:319537/1‑251 (MQ=255)                             cgTGGCATGGAAATACTCCTTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTg                                                                                                                                                                                                            <  1:157586/251‑1 (MQ=255)                                ggCATGGAAATCCTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGATTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCTTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTACCACGACCGCCACGGATCAGGATCACGGAGTGCTc                                                                                                                                                                                                         <  1:251032/251‑1 (MQ=255)                                  aatgtaaaTACTCCGTGGTTAATTCGGGTGTGTCAAAAGCTCTAGGAGGTTATGTTGACATTTAAGTTAAAACGTTTGGCCTTGTTTAGCGGAGATCCATGAAGCCTTAGGAGGCTTCACTCCATACTTGGAATGAGCCTGGTTACGGGCTTTAACGCCGGAGTATTCGAGCGCACCACTTACGGGGTGGTAACGATCACCGGGGAGTTCTCGAACATGACCGCCACGGATCAGGATCACGGAGTGCTCCt                                                                                                                                                                                                       <  1:234245/246‑1 (MQ=255)                                     gcaaaTACCCCTTTGTTAAGACAGGTTTGTCCAAAGCTTTCCGATTTAAGTTAGACGTTTAAGTTAAACCGTTTGGCCTTACTCAACGAAAAACAATTAACCTTTAGTCCGTTTCACGCCATATGTGAAACGCGTCATCTTACCTTCTTTAACGCCGGATCATTCAAGCCCCCCGCGTACTGTGGTTTAAGGAACTCCCGGGAGGTCTCTAACACGCCCGCCACGGATCAGGATCACGGAGTGCTCCTGc                                                                                                                                                                                                     <  2:328631/248‑1 (MQ=255)                                               cGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGAcc                                                                                                                                                                                          >  1:327685/1‑251 (MQ=255)                                                 ttgttAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCtt                                                                                                                                                                                        >  1:77649/1‑251 (MQ=255)                                                 ttgttAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCt                                                                                                                                                                                         <  2:321768/250‑1 (MQ=255)                                                  tgttAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCt                                                                                                                                                                                         >  1:321768/1‑249 (MQ=255)                                                   gtaaaTTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGTCCTTACTTAACGGAGAACCATTAATCCTTAGGACGCTTCACGCCATACTTGGAACTAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCTCCACGTACGGTGGGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTca                                                                                                                                                                                      <  2:181452/248‑1 (MQ=255)                                                      aaTTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTcacc                                                                                                                                                                                    >  2:62132/1‑250 (MQ=255)                                                        ttCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTcacca                                                                                                                                                                                   >  1:90563/1‑249 (MQ=255)                                                        tacAGGATTGCCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTAAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTcacca                                                                                                                                                                                   <  2:90563/247‑1 (MQ=255)                                                         tCAGGAGTGTCCAAACCTCTACAAGTTTAGTTTGACATTTAAGTTAAAACTTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTACCGCCGGAGCAGTCAAGCGCACCACGTACGGTGGGGTAACGAACACCCGGGAGGTCTCTAACACGTCCGCCACGTATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCaccacc                                                                                                                                                                                 <  2:117990/250‑1 (MQ=255)                                                                                       tttGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTc                                                                                                                                                  >  2:131970/1‑251 (MQ=255)                                                                                                               ggCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGc                                                                                                                          >  1:122777/1‑251 (MQ=255)                                                                                                                   ttacttaACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAAAGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAgcg                                                                                                                      >  1:108723/1‑251 (MQ=255)                                                                                                                             gaaaaCCATTAAGCCTTAGGACGCTTCACGCAATACTGGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTCCGGTGTGGTAACGAACCCCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCgg                                                                                                            <  2:240534/248‑1 (MQ=255)                                                                                                                                aaCCATTAAGCCGTAGGACCCTTCCCGCCATACTTGGAACGAGCCTGCTTACGGACTTTAACGCCGGATCAGTCAAGCGCACCACGTACGGTGTGTTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACTCGGCATACTTTACGCAGCGCGGAGTTCGGttt                                                                                                         <  1:313281/251‑1 (MQ=255)                                                                                                                                  ccATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCCCCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGttttt                                                                                                       <  1:201343/251‑1 (MQ=255)                                                                                                                                    aTTACGCCTTATGACGCCTCACGCCATACTTGGACCGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCCCCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTa                                                                                                     <  2:36446/251‑1 (MQ=255)                                                                                                                                         gCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGt                                                                                                >  1:163264/1‑251 (MQ=255)                                                                                                                                         gCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTACCGCAGCGCGGAGTTCGGTTTTTTAGGAGt                                                                                                >  2:181111/1‑251 (MQ=255)                                                                                                                                              aGGACGCTTCACGCTATACTGGGAACTTGCCTGCGTACGGTCTGTAACGCCGGAGCATTCAGTCGCACCCCGTACGGTGTGGTAACGAACACCGGGGAGGTCTCTAACACGACCGCCCCGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGgta                                                                                            <  2:223567/250‑1 (MQ=255)                                                                                                                                                gACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGCCTTTCACGCCGGCGCAGTCAAGCGCACCCCGTCCGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCCTACTTTACGCCGCGCGGCGTTCGGTTTTTTAGGAGTGgtagt                                                                                          >  2:32552/1‑250 (MQ=255)                                                                                                                                                   gCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGtatat                                                                                      >  1:133067/1‑251 (MQ=255)                                                                                                                                                   gCTTCACGCCATACTTGGAACGAGCCTGCTTACGGACTTTAACACCGGCCCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGAAGGACTCTAACACGACAGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTAGAAACCGTTAGTCAGACGCACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTACtatat                                                                                      >  1:177418/1‑251 (MQ=255)                                                                                                                                                     ttCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATac                                                                                    >  1:84782/1‑251 (MQ=255)                                                                                                                                                      tCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGtatata                                                                                     <  1:319537/249‑1 (MQ=255)                                                                                                                                                               aCTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCGAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTGTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACAt                                                                          <  1:176704/251‑1 (MQ=255)                                                                                                                                                                  tGGAACGAGCCTGCTTACGGTCTGTAACGCCGGAGCATGAAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACCCGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATAc                                                                        <  2:365762/250‑1 (MQ=255)                                                                                                                                                                    gAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGcc                                                                     <  2:77649/251‑1 (MQ=255)                                                                                                                                                                        gAGCCTGCTTCCGGTCTTTAAAGCCGGAGCATTCAAGCGCTCCACGTCCGGTGTGGTAACGCACCCCCGGGAGGACTCTAACACGATCGCCTCGGATCAGGGTCCCGGAGGCCACCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTAGACAACGTTCGTCAGACGACCCCGGCCGACTTTACTCAGCGCGGAGCTCGGGTTTTTAGGAGTGGCAGTATAAACACGAGTaacaacgccaagc                                                                 >  2:233924/1‑239 (MQ=255)                                                                                                                                                                             tGCTTACGGTCTTTAACGCCGGAGCAGGCAAGCGCACCACGTACGATGTGGTAAAGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCAACGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGAGGAGTTCGGTTTTTTCGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTg                                                            >  2:60671/1‑251 (MQ=255)                                                                                                                                                                              gCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGc                                                           >  1:221282/1‑251 (MQ=255)                                                                                                                                                                                   cGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGc                                                       >  1:365721/1‑250 (MQ=255)                                                                                                                                                                                     gTCTTTAAGCCCGGAGCAGTCAAGCGCACCACGTACGGTGGGTTACCGAACACCCGGGAGGTTTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCTGTTTTTTAGGAGTGGTAGTATTTACACGAGTACATACGCCACGTTTTTGCGGGCATg                                                    <  1:124653/251‑1 (MQ=255)                                                                                                                                                                                                                accacGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCg                         >  1:73651/1‑251 (MQ=255)                                                                                                                                                                                                                         gggggggTAACGAACCCCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGc                <  2:365721/246‑1 (MQ=255)                                                                                                                                                                                                                                       cacCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTg  <  1:32552/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             AATTTTACGCTGACCAATGACGCGACGACGTGGCATGGAAATACTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTG  >  NC_000913/3474046‑3474526

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