Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,482,901	C→T	V313M (GTG→ATG)	valS ←	valine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,482,901	0	C	T	97.3%	227.1 / NA	73	V313M (GTG→ATG)	valS	valine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (38/33);  minor base G (0/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCTTTGTAACGCGGATCTTCCGGGTTAACGG  >  NC_000913/4482657‑4483137                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             aCCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTGCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGgc                                                                                                                                                                                                                                        <  1:270756/251‑1 (MQ=255)             gCCAGGCCATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGat                                                                                                                                                                                                                            >  2:297367/1‑251 (MQ=255)                aGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGCTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGCACATCTGGGCGCTTTCACGGATATc                                                                                                                                                                                                                         >  2:57018/1‑251 (MQ=255)                   aCATCGGCACGCTCGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTGCGTCCAGCAGGCCAAGCGCGGCATCTGCGGCAACGACTGCGTTACGTGCAGCAAATCGCTCCATTTTCTGGAACTCTGCATGGAGTTGGGTGGGAGAAACTTCAGATGCGTTACCTTTGGTATCAAACATCTGGgctgttgtagggaggtcgcc                                                                                                                                                                                                                      <  2:382513/251‑19 (MQ=255)                    cATCGGCACGAACGTACGACTGGTCGGGCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGTTTTGCGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTTGTATCGAACATCGGGGCGCTTTCACGGATATCGCCg                                                                                                                                                                                                                     <  2:358038/251‑1 (MQ=255)                     aTCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGt                                                                                                                                                                                                                    >  1:1076/1‑251 (MQ=255)                               ccgcactccTGGCCGGTCAGCAGCGGTTCGATAACACCGCCGCCACGGCCACCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCCGGAATAAACGGCAGATTCGTTACCTTGGGTATCGAACATCTGTGCGCTTGCACGGATATCGCCCTCAAAGGTCAg                                                                                                                                                                                                          <  2:51530/243‑1 (MQ=255)                                 gTACCACTGGTAGGTCAGCATTGGTTCGATACCTCCGCCGCCACGGTCGCCGTACGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCACAGGTCAgg                                                                                                                                                                                                         >  2:370410/1‑250 (MQ=255)                                   aCCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATg                                                                                                                                                                                                      >  2:389929/1‑251 (MQ=255)                                   aCCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGACAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATg                                                                                                                                                                                                      >  1:180040/1‑251 (MQ=255)                                    ccACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGt                                                                                                                                                                                                     <  1:304983/251‑1 (MQ=255)                                            cggtcAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCATCGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATAGCTCCGTCAAAGGTCAGGATGTTGatcat                                                                                                                                                                                              <  1:296874/250‑1 (MQ=255)                                              gtcAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCgg                                                                                                                                                                                           >  2:232110/1‑251 (MQ=255)                                                 aGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTGCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGGTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAg                                                                                                                                                                                        <  1:159312/251‑1 (MQ=255)                                                  gccttGGTTCGATAACTCCGCCGCCACGGTCGCCGTAAGGTACAGTCAGGTCGTGCGGTTTAATTTCTTCCAGCATGCCAAGCGCGTCAACTGCGGCAACGACTGATTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCAGGGCGCTTTCACGGATATCGCCGTCGAAGGTCAGGATGTTGATCCTCGGCAgg                                                                                                                                                                                       >  2:135201/4‑251 (MQ=255)                                                   cattgccacaaTAACTACGCCCCCACGCTCGCCGCTAGGCACGGTCAGGTCGTACCGTTTAATTTCTTCCACCAGGCCAAGCGCGTCAACTCCCCCAACGACTGCTTTACGTGCAGCAAGACGCTCCAGCTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTCTCGACCATCTTGGCGCTTTCACGGATATCGCCGTCGAAGGTCAGGATGTTGATCAGCGGCAgcg                                                                                                                                                                                      >  2:163890/11‑249 (MQ=255)                                                    aTTGGTTCGATAACTACGCCGCAACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTATTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTAAGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGc                                                                                                                                                                                     <  1:35563/251‑1 (MQ=255)                                                        gTTCGATAACTACGCCGCCACGGTCGCAGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTg                                                                                                                                                                                  >  2:398149/1‑250 (MQ=255)                                                          gcGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACg                                                                                                                                                                               <  2:181399/250‑1 (MQ=255)                                                           cGATAACTTCGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCATCAACTGCGGCAACGACTTCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACGCTGCAGGGATTTCGCTGGAATAAACGTCATATGCGGTACCTTTGGTAGCGAACTTCTGGTCGCTTTCACGGGTATCGGCGTCAAAGGTCGGGATGTTGATCATCGGCTGTGCGTGACGt                                                                                                                                                                              <  2:260182/251‑1 (MQ=255)                                                                                   cGTAAGGAACGGTCAGGTCGTTCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGCCTGCTTTACGTGCAGCAAAACGCTCCAGTTTCGGGAACTCTGCAGGGATTTCGCTTGAATAAACGTCAGATTCGTTACCTGTGGTATCGAACATCTGGGCTCTTTCACGGAGATCGCCGTCAAAGGTCAGGATGTTTTTCATCGGCATGGCTTGACGTTTACCCACTTCATAGTCGTTAAAg                                                                                                                                                      <  2:381620/251‑1 (MQ=255)                                                                                           aCGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTgcgc                                                                                                                                              <  1:162125/251‑1 (MQ=255)                                                                                                 aGGTCGTGCGGGTTAATCTCTTCCGGCAGGCCAAGCGCGTCAGCTGCGGCAACAACTGCTTTACGTGCAGCAAACCGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCGGGAATGAACGTCAGAGTCGTTACCTTTGGTATCGAACATCTGGGCGCTGTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTAACGTTTACCCGCTTCATAGTCGTTAAAGTCGTGCGCCGGAGt                                                                                                                                        <  2:214434/251‑1 (MQ=255)                                                                                                  ggTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTg                                                                                                                                       <  1:398149/251‑1 (MQ=255)                                                                                                  ggTCGTGCGGTTTAATTTCTCCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGTCGGGACGTTTACCCACTTCATAGTCGTTAAGGTCGTGCGCCGGAGTg                                                                                                                                       <  1:28754/251‑1 (MQ=255)                                                                                                     cGTGCAGTTTAATTTCTTCAACCAGGCCAAGCGCGTCAACTGCGGAAACGACAGCTTTACGTGCAGCACAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGCTTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGTATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGCCGATTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGGGTGAtt                                                                                                                                    >  2:205999/1‑251 (MQ=255)                                                                                                      gTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGAttt                                                                                                                                   <  2:112054/251‑1 (MQ=255)                                                                                                       tGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGCTTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGAtttt                                                                                                                                  >  2:89491/1‑251 (MQ=255)                                                                                                        gCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTc                                                                                                                                 >  1:323336/1‑251 (MQ=255)                                                                                                             ttAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGc                                                                                                                             <  2:333327/250‑1 (MQ=255)                                                                                                               aaTTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCAACTTCCTAGTCGTTAAAGCCGTGCGCCGGAGTGATTTTCACGCAGc                                                                                                                          >  2:85339/1‑251 (MQ=255)                                                                                                                aTTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTGTGGTATCGAACAGCTGGGCGCTTTCACGGATATCGCCGGCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCTCTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGcc                                                                                                                         <  1:262308/251‑1 (MQ=255)                                                                                                                     ttcCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGc                                                                                                                    >  1:301501/1‑251 (MQ=255)                                                                                                                            agGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTc                                                                                                             >  1:294007/1‑251 (MQ=255)                                                                                                                               ccAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCAt                                                                                                          <  1:91717/251‑1 (MQ=255)                                                                                                                                    cgcgTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCgg                                                                                                     <  1:294142/251‑1 (MQ=255)                                                                                                                                    cgcgTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGGCAGATTCGTTACTTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTccttggcgg                                                                                                     <  1:294131/251‑8 (MQ=255)                                                                                                                                       gTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCgg                                                                                                  >  2:311390/1‑250 (MQ=255)                                                                                                                                         cAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCCTTTCGGCgtgt                                                                                                >  2:289047/1‑251 (MQ=255)                                                                                                                                              gCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATGTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTtcgt                                                                                            <  1:39434/250‑1 (MQ=255)                                                                                                                                                 gCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGcc                                                                                        >  1:367646/1‑251 (MQ=255)                                                                                                                                                  cAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTCAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGcc                                                                                        >  2:97991/1‑250 (MQ=255)                                                                                                                                                      gACTGATTTACGTGCAGCAAAAAGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGGCGTGCGCCGGAGTGCTTTTCCCGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACgc                                                                                   >  2:172355/1‑250 (MQ=255)                                                                                                                                                             ttACGTGCAGCAAAACGCTCCAATTGCTTGAACTCTGCAGGGATTTCGCTTGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTGTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTGTACCCTCGTCATAGTCATTAAAGTCGGGCGCCTGAGTGATTATCACGCAGCCGGTGTCTTTTTCCATGTCGGCGTGTACGTCGCCAACATTCGGaa                                                                             <  2:244635/250‑1 (MQ=255)                                                                                                                                                               aCGTTCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTTCCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCTTTAAAGTCGTGCGCCGGAGTGATGTTCACGCAGCCGTTGCCTTGTTCCATGTCGGCGTGTTCGTCGCCAACGATCTGAATac                                                                          <  2:367737/251‑1 (MQ=255)                                                                                                                                                               aCGTGCAGCAAAACGCTCCATTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATGCGTTACCTTTGGTATCGATCATCGGGGCGCTTTCACGGAGATCGCTGTCAAAGGTCAGGATGTTGATCAGCGGCAGGGCGTGACGTTGACCCACTTCATAGTCGTTAAGGTCGTGCGCCGGAGTTATTTTCACTCAGCCGTTGCCTTTTGCCATGGCGGCTTGTTTGGCGCCAACGATCGGAATa                                                                           <  2:27994/250‑1 (MQ=255)                                                                                                                                                               aCGTGCAGCAAAACGCTACAGTTTCTGGACCTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCCCCAACGATCGGAATac                                                                          >  2:367331/1‑251 (MQ=255)                                                                                                                                                                 gggcagcaAAACTCTCCTGTTTCTGGAACTCTGCATGGATGTCGCTGGAATAAACGTCAGATTCGTTACCTTGTGTAGCGTACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGTATGGGGATCATCGGCCGGGCGTGACGTTTACCCGCTTTATAGTCTGTAAAGTCGTGCGCCGGAGTCATTTTCACACAGATGGTGCCTATTTCCATGTCGGCGTGTTCGTCTCCACCGTTCGTAATacta                                                                        <  2:261710/249‑3 (MQ=255)                                                                                                                                                                 gTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCCGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGGCGCCCACGATCGGAATacg                                                                         >  2:125113/1‑250 (MQ=255)                                                                                                                                                                   gcagcaAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCa                                                                                                                            >  1:294196/1‑197 (MQ=255)                                                                                                                                                                   gcagcaAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCAGAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCa                                                                                                                            <  2:294196/197‑1 (MQ=255)                                                                                                                                                                        aaaaCGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACCTCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCCTCGGCAGGGCGTGACGCTTCCCCACTTCCTAGTCGTTACAGTCGTGCGCCGGAGTGTTTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTa                                                                  >  1:400546/1‑250 (MQ=255)                                                                                                                                                                        aaaaCGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTa                                                                  >  2:151099/1‑250 (MQ=255)                                                                                                                                                                                       caggatctctGCAGGGAATTCACTGGGATAAGCGTCAGATTAGTTACCTTTGGCATCGACCATCTGGGCGCTTTCACGAGTATCGCCGACAAAGGTCAGGATGTTGATCATCGCCAGGGCGTGACGTTTACCCAATTCCCAGTCGTTCAAGACGTGCGCCGGAGAGATTTTCACTCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGGCGCCAACGATCGGACTAAGACGTTTTACCAGCgctagtacg                                                   >  2:367684/7‑242 (MQ=255)                                                                                                                                                                                          gAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTGGGGATCGAACATCTGGGCGCTGTCACGGATAGCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACat                                               <  2:1076/251‑1 (MQ=255)                                                                                                                                                                                           aaCTCTGCAGGGATTTCGCTGGAAGAAACGTCAGATGCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACata                                              <  2:324345/251‑1 (MQ=255)                                                                                                                                                                                                  cAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGcc                                        >  1:91314/1‑250 (MQ=255)                                                                                                                                                                                                       aTTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATc                                    >  2:167051/1‑249 (MQ=255)                                                                                                                                                                                                           cGCTGGAATCAACGTCAGATTCGTTACCTTTGGTATCGACCATCTGGGCGCTTTCACGGATCTCGCCGTCATAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGATAAAGTCGTGCGCCGGAGTGATTTTCCCGCAGCCGGTCCCTTTTTCCATGTTGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATCACATATTTGCCAATCAGaact                              >  2:138867/1‑248 (MQ=255)                                                                                                                                                                                                             cTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCttt                            >  2:368710/1‑251 (MQ=255)                                                                                                                                                                                                              tGGAATAAACGTCAGATTCGTTCCCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCAACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGCTCGGCATACGAAGGTTAACCAGCGGCCGAATGACATATTTGCCAATCAGATCTTTg                           >  1:93614/1‑251 (MQ=255)                                                                                                                                                                                                              tGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCTTTg                           >  1:270776/1‑251 (MQ=255)                                                                                                                                                                                                                gAATAAAAGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGAGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCAATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGAGGAAGAATGACATATTTGCCAATCAGCTCTTTGt                          >  1:148167/1‑250 (MQ=255)                                                                                                                                                                                                                  aTAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCTTTGTAAc                       <  1:389929/251‑1 (MQ=255)                                                                                                                                                                                                                   tAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCTTTGTAAcg                      >  1:101299/1‑251 (MQ=255)                                                                                                                                                                                                                        gTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCTTTGTAACGCGGAt                 <  2:294007/251‑1 (MQ=255)                                                                                                                                                                                                                             aTTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAACGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGCGTGCTTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCTTTGTAACGCGGGTCTTcc            >  2:38092/1‑251 (MQ=255)                                                                                                                                                                                                                              gtcggtacgtttGGTATAGAACAGCTGGGCGCTTTCACGGATAGCGCCGTCAAGGGTGAGGGTGTGGATCATCGGGGGGGCGGGACGTTTACCGACTTCGTAGGCGTTAAAGTCGTGGGGCGGAGTGATTTGCACGCAGCCGGTTGCTTGTTGCATGTCGGCGGGTTCGGCGCGATCGATCTGAAGACGACGGTTAACCAGTGGCAGAATGACATATTTGCCAATCAGATCTTTGTAACGCGGATCTTCCg           <  1:268603/242‑1 (MQ=255)                                                                                                                                                                                                                                 gTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCa                                                        <  2:132027/203‑1 (MQ=255)                                                                                                                                                                                                                                 gTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCa                                                        >  1:132027/1‑203 (MQ=255)                                                                                                                                                                                                                                   tACCTTTGGTATCGAACATCTGGGCGCTGTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCTGCAGAATGACATATTTGCCAATCAGATCTTTGTAACGCGGATCTTCCGGGtt       <  2:186394/250‑1 (MQ=255)                                                                                                                                                                                                                                    aCCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCTGCGGGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCTTTGTAACGCGGATCTTCCGGGTTa      <  1:350160/250‑1 (MQ=255)                                                                                                                                                                                                                                       tttGGTCTCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCATGGCGTGACGTTTACCCAATTCATAGTCGTTAAAGTCGTGCGCCGGAGGGATTTTCACGCCGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCCGAATACCACGGTTAACCAGCGGCAGAATGACATCTTTGCCAATCAGATCTTTGTAACGCGTATCTTCCGAGTTAACgg  >  2:98899/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             ACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCTTTGTAACGCGGATCTTCCGGGTTAACGG  >  NC_000913/4482657‑4483137

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