Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,148	G→A	C134C (TGC→TGT)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,148	0	G	A	98.3%	201.0 / NA	58	C134C (TGC→TGT)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/30);  minor base C (1/0);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCTACAATGGATTCGATTCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCGGCACGCCTTCTTCAATAGAGGCGGTACACAACA  >  NC_000913/3815860‑3816440                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                     tCGCTACAATGGATTCGATTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCNCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCAccc                                                                                                                                                                                                                                                                                             >  1:517944/1‑299 (MQ=255)    cTACAATGGATTCGATTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAATACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTaa                                                                                                                                                                                                                                                                                         >  1:195306/1‑300 (MQ=255)      aCAATGGATTCGATTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATc                                                                                                                                                                                                                                                                                       <  1:375628/300‑1 (MQ=255)       cAATGGATTCGATTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATc                                                                                                                                                                                                                                                                                       >  1:646295/1‑299 (MQ=255)        aaTGGATTCGATTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCACTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCGGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCAACGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGa                                                                                                                                                                                                                                                                                     >  2:618341/1‑300 (MQ=255)                             ggCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTATCGCATCTGCCCGCGCCCCACAGGCACCCGTACTCGCAGCGGTACGCGTGCCCCCATc                                                                                                                                                                                                                                                                >  1:439526/1‑300 (MQ=255)                             ggCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATc                                                                                                                                                                                                                                                                >  1:34525/1‑300 (MQ=255)                                       aacagggttagtagtTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCGGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGc                                                                                                                                                                                                                                                      <  2:298596/291‑1 (MQ=255)                                       aacaaGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGc                                                                                                                                                                                                                                                      >  1:41347/1‑300 (MQ=255)                                       aacaaGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGc                                                                                                                                                                                                                                                      <  2:565483/300‑1 (MQ=255)                                       aacaaGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTACACTGCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCCGAAGc                                                                                                                                                                                                                                                      <  2:202159/300‑1 (MQ=255)                                                   aGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCATACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTcc                                                                                                                                                                                                                                          <  2:43048/300‑1 (MQ=255)                                                   aGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTcc                                                                                                                                                                                                                                          <  1:491084/300‑1 (MQ=255)                                                    gCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCa                                                                                                                                                                                                                                         >  1:235006/1‑300 (MQ=255)                                                      tctTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAg                                                                                                                                                                                                                                        <  2:570931/299‑1 (MQ=255)                                                       ctTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCATTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGc                                                                                                                                                                                                                                       >  2:659603/1‑299 (MQ=255)                                                        ttCATGGGTGAACGGCTCGCCTTCTGCCGTCTCCTGGACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATTTCGTTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGt                                                                                                                                                                                                                                     <  2:28150/300‑1 (MQ=255)                                                             gggTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTTTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATgg                                                                                                                                                                                                                                <  2:33739/300‑1 (MQ=255)                                                             gggTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATgg                                                                                                                                                                                                                                <  1:604074/300‑1 (MQ=255)                                                               gTGAACGGCTCGCCTTCTGCCGTCCCTTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCACATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTcc                                                                                                                                                                                                                                          <  2:42931/288‑1 (MQ=255)                                                               gTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCTTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTg                                                                                                                                                                                                                              <  2:692999/300‑1 (MQ=255)                                                                        tCGCCTTCTGCCGTCCCCTGCACTTCAATTATCCGCCCGTCTTCAGTCATCACTACTTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGAGCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTTAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCAccc                                                                                                                                                                                                                     <  2:159975/300‑1 (MQ=255)                                                                        tCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCAccc                                                                                                                                                                                                                     <  1:558380/300‑1 (MQ=255)                                                                             ttCTGCCGTCCCCTGCACTTCACTGATGCGCCCGTCTTCGGTCATCACTACGTGCATGTCGGTCTCTTCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCTACAGAAACCGCTGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTTTAGCACATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCACAGTCCAGCGTAATGGTGAACTCACCCAgcg                                                                                                                                                                                                                 <  2:368708/299‑1 (MQ=255)                                                                             ttCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAgcg                                                                                                                                                                                                                 <  1:643763/299‑1 (MQ=255)                                                                                      ccccTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATc                                                                                                                                                                                                       >  2:434812/1‑300 (MQ=255)                                                                                                              gTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACggg                                                                                                                                                                                <  1:686845/299‑1 (MQ=255)                                                                                                                  tCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTTTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATc                                                                                                                                                                           <  2:94159/300‑1 (MQ=255)                                                                                                                    ggTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACATAGGCACCCGTAATCGAAGCGGTCCGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACCCAACCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCa                                                                                                                                                                          >  2:387945/1‑299 (MQ=255)                                                                                                                    ggTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGTGTGCCACCATCCGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCAAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGAGATCa                                                                                                                                                                          >  2:388516/1‑299 (MQ=255)                                                                                                                    ggTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAAATCGAAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCCTCTGCAAGCGCCACACAGGCACCCGCAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAGTGGTGAACTCACCCAGCGCTTTCACATCTACTGCCGCGCGACGTGCACGGGCGATc                                                                                                                                                                           >  2:387748/1‑298 (MQ=255)                                                                                                                     gTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGa                                                                                                                                                                        >  1:126732/1‑300 (MQ=255)                                                                                                                         tcaCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTAACCGTTTACAATACCGACAGAAACCGAGGCTACCATAAATTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTCGCGCATCTGCCAGCGCCACCCCGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGACGCACGTCGCAGTCCAGCGTAATGGTGACCTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGTACGGGCGATCAGACGCt                                                                                                                                                                    >  1:215012/1‑300 (MQ=255)                                                                                                                           aCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCt                                                                                                                                                                    >  1:117632/1‑298 (MQ=255)                                                                                                                           aCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTNTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCNCCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATAAGCCTGAAGCACGTCGAAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCt                                                                                                                                                                    >  2:519519/1‑298 (MQ=255)                                                                                                                                            ctctGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTAACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTCCTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTcc                                                                                                                                                 >  1:485755/1‑300 (MQ=255)                                                                                                                                               tGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTTTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCAc                                                                                                                                               <  2:605928/299‑1 (MQ=255)                                                                                                                                               tGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCAc                                                                                                                                               <  1:613007/299‑1 (MQ=255)                                                                                                                                                gctgctgagTCGTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCAccc                                                                                                                                             <  1:100757/294‑1 (MQ=255)                                                                                                                                                        gTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTAcc                                                                                                                                     >  2:328552/1‑300 (MQ=255)                                                                                                                                                        gTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCAACCTGCTTAcc                                                                                                                                     >  2:447753/1‑300 (MQ=255)                                                                                                                                                        gTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCATTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTAcc                                                                                                                                     >  2:448038/1‑300 (MQ=255)                                                                                                                                                           ttCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCtt                                                                                                                                   <  1:549901/299‑1 (MQ=255)                                                                                                                                                               aCGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTcgct                                                                                                                              >  1:82738/1‑299 (MQ=255)                                                                                                                                                               aCGTATTCCAGATCGCAAATTGCTTCAACGTTTACAATACCGACAGAAACCGCGGCTCCCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGTGCATCTGCCAGCGCCAAACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCCCCCTGCTTCCCTTTcgct                                                                                                                              >  1:57401/1‑299 (MQ=255)                                                                                                                                                                                              ttacaATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAAc                                                                                               <  2:441972/298‑1 (MQ=255)                                                                                                                                                                                                                     gCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGc                                                                         >  1:401979/1‑299 (MQ=255)                                                                                                                                                                                                                       tACCATCCCTTTTATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTTAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGc                                                                         <  2:197517/297‑1 (MQ=255)                                                                                                                                                                                                                                 ttCATTGGATTTGTTTTCAGCTTGCCTTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTTAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCACTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAgcc                                                             <  2:348675/299‑1 (MQ=255)                                                                                                                                                                                                                                 ttCATTGGATTGGTTTTCAGCTTGCCGTGTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGTGCGATCAGACGCTTGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGc                                                                         <  2:198735/287‑1 (MQ=255)                                                                                                                                                                                                                                      tggATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTTCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAgccctggcc                                                       <  2:646295/299‑1 (MQ=255)                                                                                                                                                                                                                                           tGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGAc                                                  <  1:628907/300‑1 (MQ=255)                                                                                                                                                                                                                                            ggTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGAcc                                                 >  1:638481/1‑300 (MQ=255)                                                                                                                                                                                                                                                              ttCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCAc                               <  2:704601/300‑1 (MQ=255)                                                                                                                                                                                                                                                                     aGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCcttct                        >  2:613797/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                      tgccAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTaca       >  2:637835/3‑300 (MQ=255)                                                                                                                                                                                                                                                                                          aGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACacaat   <  1:322214/300‑2 (MQ=255)                                                                                                                                                                                                                                                                                           gcgcCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGNACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACacaata  >  1:510947/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                     TCGCTACAATGGATTCGATTCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCGGCACGCCTTCTTCAATAGAGGCGGTACACAACA  >  NC_000913/3815860‑3816440

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑