Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,441,065	Δ1 bp	coding (611/621 nt)	rpsD ←	30S ribosomal subunit protein S4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,441,065	0	T	.	97.7%	195.6 / ‑3.0	45	coding (611/621 nt)	rpsD	30S ribosomal subunit protein S4
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (21/22);  total (22/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 = 9.96e-01

CATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGATAAGAAGCGATGTTAACAACACGACCGTTTACCAT  >  NC_000913/3440831‑3441306                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  cATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCa                                                                                                                                                                                                                                    <  1:147179/249‑1 (MQ=255) cATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCa                                                                                                                                                                                                                                    <  1:244135/249‑1 (MQ=255)      aatactccATCAATCTCAACCTCGGTCTCCTCGCAACCCGGCAGCGATGAGAGGAGAATACGTCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTGACCTTGGCGTGCGGCTAACTCACTTGCTCGATATCAACCAGTCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCGTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTTGACGATCTGGTGTTc                                                                                                                                                                                                                             <  2:415951/245‑1 (MQ=255)       gtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCg                                                                                                                                                                                                                            <  1:429404/251‑1 (MQ=255)         acacCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGtt                                                                                                                                                                                                                          >  1:102489/1‑251 (MQ=255)                              gTCACCGAGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCCACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATc                                                                                                                                                                                                        >  2:307039/1‑248 (MQ=255)                                    gcgcAACCCGGCATAGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTc                                                                                                                                                                                               >  1:107568/1‑251 (MQ=255)                                        aaCCCGGGATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCTAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGGCGAACTCACTTGCTCTATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGAGCAGGTGTTCGTGAATGGCCGCAGACAGATCATAACGCTCCGGc                                                                                                                                                                                           <  2:306499/251‑1 (MQ=255)                                             ggCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACGCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCGCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACg                                                                                                                                                                                      <  2:44806/251‑1 (MQ=255)                                                aTCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCtt                                                                                                                                                                                   >  1:127094/1‑251 (MQ=255)                                                       agagCAGCATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGGCCGCAGACAGCTCAGAACGCTCCGGCTTACGCTTAAAcat                                                                                                                                                                             >  2:444751/1‑248 (MQ=255)                                                        gagCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTAc                                                                                                                                                                           >  1:384189/1‑251 (MQ=255)                                                         aacAGAATACGGCGCAGTGCGTTCCCCCGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGGGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTACGCTTTACTTGG‑GGAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGAAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTAcc                                                                                                                                                                          >  2:130813/3‑251 (MQ=255)                                                           cAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCtt                                                                                                                                                                        >  1:72689/1‑251 (MQ=255)                                                             gAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAACGCCACGCTCTAAATGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTAACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTACGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTcc                                                                                                                                                                      >  1:413156/1‑251 (MQ=255)                                                                        aGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGCCAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGc                                                                                                                                                           >  2:412199/1‑251 (MQ=255)                                                                            cGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCGTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCa                                                                                                                                                       <  2:170521/251‑1 (MQ=255)                                                                                  ccAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCTCTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTcc                                                                                                                                                 <  1:252708/251‑1 (MQ=255)                                                                                   cAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTcc                                                                                                                                                 >  2:11878/1‑250 (MQ=255)                                                                                    agagTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTcc                                                                                                                                                 >  1:469347/1‑249 (MQ=255)                                                                                              cAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGc                                                                                                                                     <  1:290589/251‑1 (MQ=255)                                                                                                aaGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCtt                                                                                                                                   >  1:2052/1‑251 (MQ=255)                                                                                                 aGCCACGCTCTAAAGGCTCAAGGGTCACCTTTGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCttt                                                                                                                                  <  2:229906/251‑1 (MQ=255)                                                                                                       gCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgc                                                                                                                            >  1:6701/1‑251 (MQ=255)                                                                                                         tctAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgctg                                                                                                                          >  1:440827/1‑251 (MQ=255)                                                                                                            aaaGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgctgct                                                                                                                        <  1:272168/250‑1 (MQ=255)                                                                                                                          cACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGA‑TAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTTCTCAGCCAGCTCCagag                                                                                                         <  2:361294/251‑1 (MQ=255)                                                                                                                            ccTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCg                                                                                                       >  1:303056/1‑251 (MQ=255)                                                                                                                                  cgtgcgtCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTGTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTc                                                                                                 <  1:451196/251‑1 (MQ=255)                                                                                                                                          gAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCAcgcg                                                                                            >  2:412751/1‑248 (MQ=255)                                                                                                                                            aCTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACt                                                                                       <  1:56825/251‑1 (MQ=255)                                                                                                                                                   tGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTGGGTACTAAGCTTTACGTGGAG‑AAAGCTCGACGATCAAGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGGTTACGCTTAAACGTTCCTTCCAGCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTGTTCACGCGGCGCAGCCAGCTCCAGAGCGGCTTGCACGCGAGACTGcttct                                                                                 <  1:292649/250‑1 (MQ=255)                                                                                                                                                   tGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCGGCATTGTGTCCTCTCGTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTGCGCTTAAACGTGCCTTCCATCGTGCCAGCATCAGCTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTGCAGAGCGGCTTTCGCGCGGGACTGcttc                                                                                  <  2:148292/249‑1 (MQ=255)                                                                                                                                                       cGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCGTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTGTTCGCGCTGCTCAGCCATCTCCAGAGCGTCTTTCACGCGATACTGCTTCTTCg                                                                              <  2:463874/249‑1 (MQ=255)                                                                                                                                                                   ggCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATg                                                                >  2:433892/1‑251 (MQ=255)                                                                                                                                                                      gcgGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCtt                                                             >  2:234376/1‑251 (MQ=255)                                                                                                                                                                       cgGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCtt                                                             >  1:21097/1‑250 (MQ=255)                                                                                                                                                                        ggTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTAc                                                           >  1:172022/1‑251 (MQ=255)                                                                                                                                                                           tttAGAATCTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTTGGACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCTGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACTCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACaa                                                         <  1:454589/250‑1 (MQ=255)                                                                                                                                                                              aGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTc                                                     <  1:153188/251‑1 (MQ=255)                                                                                                                                                                                 aaCTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCAtt                                                  <  1:248566/251‑1 (MQ=255)                                                                                                                                                                                     atGTCACAGAACCCTGCATTGAGTCCTCTCTATGGGACTAAGCTTTACTTGGAG‑AAAGCGCGACGATCAGGTGTTCGTGAATGTCCGCAGACAGATCAGAACGCTTCGTCGTACGCTTAAACGTCCCTTCGGTCTTGCCAGCAGCGACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGTCTTTCACGCGATACTGCTGCTTCGCTTTCTCGCGAATGCTTACAACGTCATTCgt                                               <  2:286384/249‑2 (MQ=255)                                                                                                                                                                                                         gtgtCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGATAAGAAGCGAt                          <  1:4724/251‑1 (MQ=255)                                                                                                                                                                                                                    ttGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGGCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGGCATTCGGACTAACCTGATAAGAAGCGATGTTAacaaca                <  2:60446/250‑1 (MQ=255)                                                                                                                                                                                                                                 ttACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGCTAAGAAGCGATGTTAACAACACGACCGTTTACCAt  >  1:280033/1‑251 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  CATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGATAAGAAGCGATGTTAACAACACGACCGTTTACCAT  >  NC_000913/3440831‑3441306

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