Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,185,777	C→T	S143F (TCC→TTC)	rpoC →	RNA polymerase subunit beta'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,185,777	0	C	T	98.5%	231.5 / NA	68	S143F (TCC→TTC)	rpoC	RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (42/25);  minor base G (0/1);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTG  >  NC_000913/4185536‑4186021                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     gggCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGACTTCTATGTGGt                                                                                                                                                                                                                                             >  1:37808/1‑251 (MQ=255)                 aCGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGGGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATg                                                                                                                                                                                                                             >  2:80731/1‑251 (MQ=255)                 aCGAGTGCCTGTGCGGTAAGTACAAGCGCCGGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGTCGGTAt                                                                                                                                                                                                                              <  2:284838/250‑1 (MQ=255)                          tgtgCGGTAAGTACAAGCGTCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTg                                                                                                                                                                                                                    >  2:266814/1‑251 (MQ=255)                          tgtgCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGACGTGAGCGTATGGGCCACATCGAACTGGATTCCCCGACTGCGCACATCTGGATCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTg                                                                                                                                                                                                                    >  1:182450/1‑251 (MQ=255)                                     tACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTcagc                                                                                                                                                                                                          >  1:22810/1‑250 (MQ=255)                                        aaGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTcagc                                                                                                                                                                                                          >  1:48067/1‑247 (MQ=255)                                                 aaaCACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAACGTACGCCGTGAGCGTATGGGCCACATAGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGACGGCGGTATGCCCAACCTGGAACGTCAGCAGATCctgactg                                                                                                                                                                                              >  2:454053/1‑250 (MQ=255)                                                  aaCACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCctgactga                                                                                                                                                                                             >  1:275886/1‑250 (MQ=255)                                                        gTGGCGTCATCTGTGAGAAGTGAGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGc                                                                                                                                                                                        >  1:406017/1‑249 (MQ=255)                                                               cATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGCTATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGACGAGCAGTATCTgg                                                                                                                                                                               >  1:22722/1‑251 (MQ=255)                                                                    gtgAGAAGTGCGGCGTTGAAGTGACCCAGACTCAAGCACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTACCGGCCCGTATCTGTCCGCTGCTCGCTATCCCGCTCCGCGATATCGCACGCGTACTGTACTTGTACTTCTCTGTGGTTATCGCAGGCGGTATGACCAACCTGGAACGTCAGCAGATCctccctgccgcgccgtctctgcccccg                                                                                                                                                                          >  2:247235/1‑226 (MQ=255)                                                                    gtgAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGAc                                                                                                                                                                             <  2:173985/248‑1 (MQ=255)                                                                       agaAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTg                                                                                                                                                                       >  2:236432/1‑251 (MQ=255)                                                                              cggcgTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGt                                                                                                                                                                >  1:334723/1‑251 (MQ=255)                                                                               ggcgTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGCATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGCAGAGtt                                                                                                                                                               >  2:160462/1‑251 (MQ=255)                                                                                 cgTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCg                                                                                                                                                             >  1:152751/1‑251 (MQ=255)                                                                                    tGAAGTGACCCAGTCTAAAGTACGCCGTGGGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGGTATCGAAGGCGGGATGACCAACCTGGAACGGCAGCAGATCCTGACTGAAGAGCAGTAGCTGGACGCGCTGGAAGAGTTCGGt                                                                                                                                                           <  2:58462/250‑1 (MQ=255)                                                                                     gAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGCATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTg                                                                                                                                                          >  1:444374/1‑250 (MQ=255)                                                                                      aaGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAAGCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGAc                                                                                                                                                        <  2:31827/251‑1 (MQ=255)                                                                                          gACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGGAGAGTTCGGGGACGAAt                                                                                                                                                    >  2:351658/1‑251 (MQ=255)                                                                                           aCCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTAACGAAGGCGGGATGACCAACCTGGAACGTCAGCAGATACTGACTTAAGAGCGGTAGCTGGACGCGCTGGAAGAGTTCGGTGACGAAtt                                                                                                                                                   >  2:254391/1‑251 (MQ=255)                                                                                               aGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGAc                                                                                                                                               <  2:346940/251‑1 (MQ=255)                                                                                                gACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGAcg                                                                                                                                              >  1:365110/1‑251 (MQ=255)                                                                                                    aaaGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCgcagcga                                                                                                                                           >  2:373986/1‑244 (MQ=255)                                                                                                          cGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATggg                                                                                                                                    >  1:133138/1‑251 (MQ=255)                                                                                                            ccGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGCCGAATTCGACGCGAAGATggggg                                                                                                                                  >  2:181790/1‑251 (MQ=255)                                                                                                              gTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCg                                                                                                                                >  1:255322/1‑251 (MQ=255)                                                                                                                    gTATGGGCCACATCGAACTGGCTGCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTGTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCa                                                                                                                          <  1:58562/251‑1 (MQ=255)                                                                                                                        tggCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCACGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGTCGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATc                                                                                                                       <  2:37808/249‑1 (MQ=255)                                                                                                                           ccACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAgg                                                                                                                   >  2:340006/1‑251 (MQ=255)                                                                                                                             acaTCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATcc                                                                                                                      <  1:433402/246‑1 (MQ=255)                                                                                                                             acaTCGAACTGGCTTCCACGACTGCGCACATCTGGTTCCTGACATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGAGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATACTGACTGAAGAGCCGTATCTGGACGCGCTGGAAGAGTTCGGTGACGCATTCGACGCGACGATGGGGGCGGAGGCAATcc                                                                                                                      >  2:433402/1‑246 (MQ=255)                                                                                                                                tCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCGGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTTGAACGTCAGCATATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGTGGGCGGAAGCAATCCAGGctct                                                                                                               <  1:107907/250‑1 (MQ=255)                                                                                                                                          ttCCCCGACTGCGCACCTCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGc                                                                                                     >  2:215873/1‑250 (MQ=255)                                                                                                                                           tCCCCGACTGCGCACATCTGGTTCCTGAAAGCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCa                                                                                                    <  1:637/250‑1 (MQ=255)                                                                                                                                            ccccGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCAt                                                                                                   <  1:130096/250‑1 (MQ=255)                                                                                                                                             cccGCCTGCGCACATCGGGTTCGTGAAATCGCCGCCGTCCCGTATCGGTCTGCTTCTCGATATGCCTCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACTAATTCGCCGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTCAAGAGCATgg                                                                                                 <  2:76129/251‑1 (MQ=255)                                                                                                                                                gACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGCACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGAAGAATTCGACGCGAAGATGGGGTAGGAAGCAATCTAGGTTCTGCTGATGAGCATGGATc                                                                                              >  1:116937/1‑251 (MQ=255)                                                                                                                                                gACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATc                                                                                              >  1:451507/1‑251 (MQ=255)                                                                                                                                                      gcACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGc                                                                                        >  1:336469/1‑251 (MQ=255)                                                                                                                                                      gcACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCGGCTGCGCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGAGCCGGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATGCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGc                                                                                        <  2:406017/251‑1 (MQ=255)                                                                                                                                                        acaTCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCaa                                                                                      <  1:454053/251‑1 (MQ=255)                                                                                                                                                           tCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTCCTTTGAATTCTATGTGGTTATCGCAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAag                                                                                     >  1:283632/1‑249 (MQ=255)                                                                                                                                                                ttCCTGACATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGCTATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGATGAAGAGCATGGATCTGGAGCAAGAGTGCgg                                                                              >  2:62030/1‑250 (MQ=255)                                                                                                                                                                              gCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGa                                                                  >  1:49063/1‑249 (MQ=255)                                                                                                                                                                                cGTCCCGTATCGGTCTTCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTTAATTCTATGTGGTTATCGAAGTCGGTATGACCATCCTGGAACGTCAGCAGATCCTTACTGATGATCAGTATCTGGACGCGCTTGATGAGTTCGGTGACGAATTCGACGCGAAGATGGGTTCGGAAGCAATCCAGCCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGAagag                                                              <  1:53648/251‑1 (MQ=255)                                                                                                                                                                                 gTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGc                                                             >  1:42713/1‑251 (MQ=255)                                                                                                                                                                                   cccGTATCGGTCTGCTGATCGATATGCCGCTGCGCGATATCGACCGCGTACTGTACTTTGCATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGACCGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTg                                                           >  2:158335/1‑251 (MQ=255)                                                                                                                                                                                   cccGTATCGGGCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTTCTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTTGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCt                                                            <  2:18244/250‑1 (MQ=255)                                                                                                                                                                                        aTCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGa                                                      <  2:275743/251‑1 (MQ=255)                                                                                                                                                                                            gTCTGCTGCTCGATATGCCGCTGCGCGATCTCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAAcc                                                  >  1:131963/1‑251 (MQ=255)                                                                                                                                                                                                gctgctCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAATAGCAGTATCTGGACGCGCGGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGAGCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACt                                              <  1:62030/251‑1 (MQ=255)                                                                                                                                                                                                gctgctCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCa                                                 <  2:127373/248‑1 (MQ=255)                                                                                                                                                                                                 ctgctCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGACTTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGCTGGGGGCGGAAGCAATCCAGGCTATGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTc                                             >  1:272279/1‑251 (MQ=255)                                                                                                                                                                                                  tgctCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTcc                                            >  2:341203/1‑251 (MQ=255)                                                                                                                                                                                                   gctCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGCGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGCTGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGGGCAAGAGTGCGAACAGCTGCGGGAAGAGCTGAACGAAAACAACtacg                                           >  2:286795/1‑248 (MQ=255)                                                                                                                                                                                                          atGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCaa                                    <  1:339850/251‑1 (MQ=255)                                                                                                                                                                                                              cGCTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGt                                <  1:219835/251‑1 (MQ=255)                                                                                                                                                                                                                cTGCGCGATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTa                               >  1:22345/1‑250 (MQ=255)                                                                                                                                                                                                                     cgATATCGAGCGCGTACTGTCCTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGAAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGCCGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGACGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACCGCTGCGTGAAGAGCTGAAAGCAACAAACTCAGAACACAAGCGTAAAAAGc                         >  2:266006/1‑251 (MQ=255)                                                                                                                                                                                                                      gATATCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGAGCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCt                        <  1:53644/251‑1 (MQ=255)                                                                                                                                                                                                                         atCGAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGAc                     <  2:103631/251‑1 (MQ=255)                                                                                                                                                                                                                            gAACGCGTACTGTACTTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCa                   >  1:110434/1‑250 (MQ=255)                                                                                                                                                                                                                               cgcgTACTGTACTTTGAATGCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCg               <  2:22722/251‑1 (MQ=255)                                                                                                                                                                                                                                  gtactgtactTTGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGTCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCTAACAGCTGCGTGAAGAGCTTAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTAt            <  1:200480/251‑1 (MQ=255)                                                                                                                                                                                                                                   tacggtactTTGACTTCTATGAGGTTATCGCAGGCGGCATGTCCTACCTGGTACGTCAGCAGATACTGCCTGAAGAGCAGAATCTGGACGCGATGGAAGAGATCAGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTATGCTGAAGAGCAAGGAACTGGAGCAAGCGTGCGAACAGAGGCGTGAAGAGCTGATCGAAACCAACTCCGAACCCAAGCGTAAAGATCTGACCAAGCGTATc           >  2:88575/5‑251 (MQ=255)                                                                                                                                                                                                                                             tGAATTCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCATActgctg  <  1:102602/250‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     GGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACGCGCTGGAAGAGTTCGGTGACGAATTCGACGCGAAGATGGGGGCGGAAGCAATCCAGGCTCTGCTGAAGAGCATGGATCTGGAGCAAGAGTGCGAACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTG  >  NC_000913/4185536‑4186021

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