Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,439	C→T	V37V (GTG→GTA)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,439	0	C	T	98.3%	171.4 / NA	59	V37V (GTG→GTA)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/38);  minor base A (1/0);  total (21/38)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCGGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑A‑AACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG‑CATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCA  >  NC_000913/3816201‑3816685                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                          cGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCg                                                                                                                                                                                                                                               <  1:403903/251‑1 (MQ=255)    aGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGcc                                                                                                                                                                                                                                             >  1:306125/1‑250 (MQ=255)        cAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTAAACAATACTTTGTTATCGCCAAAtt                                                                                                                                                                                                                                        >  2:240380/1‑251 (MQ=255)            gTAATGGTGAAATCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTATTCAATAGAGGCGGTACACAATACTGTGGGATCGCCAAATTctac                                                                                                                                                                                                                                    >  1:88851/1‑248 (MQ=255)              aaTGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGAcc                                                                                                                                                                                                                                   <  1:200923/250‑1 (MQ=255)                  gTGACCTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTAAACAATACTTTGGTATCGCCAAAttttcccagccc                                                                                                                                                                                                                              >  2:35267/1‑241 (MQ=255)                      aCTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGTCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTGTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAg                                                                                                                                                                                                                          <  2:136764/251‑1 (MQ=255)                                    ttCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGtt                                                                                                                                                                                                             <  1:222902/250‑1 (MQ=255)                                      cAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATGTCCATTGTGCGGCCACCCGGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGtttt                                                                                                                                                                                                           <  2:139371/250‑1 (MQ=255)                                             aCTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTTCCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTCCTGCATGTTTTGTATAg                                                                                                                                                                                                     <  2:325160/249‑1 (MQ=255)                                             aCTGCCGAGCGAAGTGCACGGGCGACCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACAGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATTCTTTGGTATCGCCAAATTCGACCATCACCGAGCCTTCTGCATGTTTTGTATAg                                                                                                                                                                                                     <  2:7044/249‑1 (MQ=255)                                               tGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACCCAATACTTTGGTATCGCCAAATTCGACCCACCCCGCGCCTTCTGCCTGTTTTTTATAGttcc                                                                                                                                                                                                 >  2:30754/1‑249 (MQ=255)                                                gCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTGCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGTTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTGGTATAGTTAc                                                                                                                                                                                                 <  2:164877/250‑1 (MQ=255)                                                       gAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTAACCAATACTTTGTTATCTCCAAAATCGACCCGCACCCAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAgg                                                                                                                                                                                         >  2:80745/1‑251 (MQ=255)                                                             cACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCTCCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACg                                                                                                                                                                                   <  2:243691/251‑1 (MQ=255)                                                              aCGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTTGACCCGCACCGAGCCCTCTGCATGTTTTGGATTGTTTCGAGTAAGGGTAACgg                                                                                                                                                                                  >  2:120068/1‑251 (MQ=255)                                                               cGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACggg                                                                                                                                                                                 >  1:299716/1‑251 (MQ=255)                                                                          aCGCTGGATTTCCATTGTGCGTCCACCCTGCTTCCCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGat                                                                                                                                                                      <  2:289358/251‑1 (MQ=255)                                                                          aCGCTGGATTTCCATTGTGCGTCCACCCTGCTGACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGat                                                                                                                                                                      <  1:402223/251‑1 (MQ=255)                                                                                  tttCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCATGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGGAACGGGACGCACCTGATTATTGCt                                                                                                                                                               >  1:182864/1‑250 (MQ=255)                                                                                   ttCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTGTGTATAGTTACGAGTCAGGGTAACGGGACGCGCCTGATTAGTGCTa                                                                                                                                                              <  2:190173/250‑1 (MQ=255)                                                                                   ttCCATTGTGCGGCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACGCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCATCACCGAGCCTTCTTCATGTTTTGTATAGTTACGAGTCAGGGTAACGGTACGCATCTGATGAGTGCTa                                                                                                                                                              <  2:399401/250‑1 (MQ=255)                                                                                    tCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAggg                                                                                                                                                                                        <  1:133278/223‑1 (MQ=255)                                                                                     ccATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACgg                                                                                                                                                           <  1:90046/251‑1 (MQ=255)                                                                                         tgtgCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCct                                                                                                                                                        <  1:49218/250‑1 (MQ=255)                                                                                         tgtgCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCct                                                                                                                                                        <  1:301654/250‑1 (MQ=255)                                                                                         tgtgCGTCCACCCGGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCGACGGCct                                                                                                                                                        <  2:42037/250‑1 (MQ=255)                                                                                              gTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGGATCGCCAAATTCGGCCAGCACCGGGCCTTCTGCATGTTTTGTATAGTTACGAGTAAGGGTAAAGGGAAGCACCTGATTATTGCTACGGCCTGCTggg                                                                                                                                                  >  1:129834/1‑250 (MQ=255)                                                                                                ccACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGGATCGCCAAATTTGACCAGCACCCAGCCTTCTGCCTGTTTTGTATAGTTACGGGTCCGGGTAAAGGGCCGCACCTGATTATTGCTACGGCctgcctggcg                                                                                                                                                >  1:17819/1‑245 (MQ=255)                                                                                                    ccTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCata                                                                                                                                            <  2:176700/251‑1 (MQ=255)                                                                                                    ccTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCata                                                                                                                                            <  1:75994/251‑1 (MQ=255)                                                                                                       gCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCTCCTGATTATTGCTACGGCCTGCTGGACGCATAtt                                                                                                                                          <  1:372273/250‑1 (MQ=255)                                                                                                       gCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCatat                                                                                                                                           <  1:35267/249‑1 (MQ=255)                                                                                                           aCCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGaaa                                                                                                                                      <  1:134325/250‑1 (MQ=255)                                                                                                                   ctgcTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACATTGGTATCGACAAATTCGACCAGCACCGAGCCTTCTGGATGTTTTGTATAGTTACGAGTCAGGGTAACGGGGCGGACCTGATTATTGCTACGGCCTGCTGGACGCAAATTGAAATCTCCGGCt                                                                                                                             >  1:20466/3‑251 (MQ=255)                                                                                                                       ttCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑a                                                                                                                         >  1:331027/1‑250 (MQ=255)                                                                                                                          aCGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTCCCCAAAACTTTGGTATTGCCCAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGGGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGATGCATATTGAAATCTCCGGCTTG‑A‑AAc                                                                                                                     >  1:148260/1‑250 (MQ=255)                                                                                                                                      ggTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATGATTGCTACGGCCTGCTGGACGCAGATTGAAATCTCCGGCTTG‑A‑AACAAATGTGGCTgcg                                                                                                        <  2:29974/251‑1 (MQ=255)                                                                                                                                            ggTAGAACGCGGCAGCATACCGGACTCTGCGGTGATCCAGCCATGGCCCTGTCCTTTCAGGAAGCGAGGCTCGCCTTCTTCAATAGAGGCTGTACACACTACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTCCGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGAATATTGAAATCTCCGGCTTG‑A‑AACAAATGTGGCGGCGCATtat                                                                                                  <  2:235132/251‑1 (MQ=255)                                                                                                                                                gAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑A‑AACAAATGTGGCTGCGCATTATACgg                                                                                              <  2:51288/251‑1 (MQ=255)                                                                                                                                                gAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCCAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCCCCTGATTATTGCTACGGCCTGCTGGACGCATATTGGAATCTCCGGCTTG‑A‑AACAAATGTGGCTGCGCATTCTACgg                                                                                              >  1:184155/1‑251 (MQ=255)                                                                                                                                                 aaCGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGGGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCCCCTGATTATTGCTACGGCCTGCTGGGCGCATATTGAAATCTCCGGCTTGTA‑AACaat                                                                                                                  >  2:72180/1‑230 (MQ=255)                                                                                                                                                 aaCGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑A‑AACAAAt                                                                                                                 <  1:72180/231‑1 (MQ=255)                                                                                                                                                 aaCGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑A‑AACAAATGTGGCTGCGCATTATACGGa                                                                                             <  1:278577/251‑1 (MQ=255)                                                                                                                                                            ataccGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGATGCATATTGAAATCTCCGGCTTG‑A‑AACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGt                                                                                  <  2:148260/248‑1 (MQ=255)                                                                                                                                                              accGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑A‑AACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTa                                                                                <  2:331027/250‑1 (MQ=255)                                                                                                                                                                 gTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACCATAATTTTGTATCGCCAAATTCGACCAGCACCGGGCCTTCTGCATGTTTTGTCTAGTTACGAGTCAGGGTAAAGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑A‑AACAAATTTGGCTGCGCATTATACGGACTTCCGGCGGTTAttt                                                                             >  2:416108/1‑250 (MQ=255)                                                                                                                                                                   aCTCGGCGGTGATCCAGCCCTGGCCCTGCCCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTGCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCGACGTCCTGCTGGACGCATATTGAAATCGCCGGCTGG‑T‑AACAAAGGTGGCTTTGCATTATACGGACTTCCGGCGGGTATTCCt                                                                           <  2:396543/251‑1 (MQ=255)                                                                                                                                                                         cGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACACTACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCCTGTTTTGTATAGTTACGAGTCAGGGGAACGGGACGCACCTGATTATTTCTACGGCCTGCTGGACGCATATTGCAATCTCCGGCTTG‑A‑AACAAATGTGGCTGCGCATTATACGGACTTCccgcgggtttttctattcct                                                                     >  2:142526/1‑237 (MQ=255)                                                                                                                                                                           gTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCATGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑A‑AACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc                                                                   <  1:69555/251‑1 (MQ=255)                                                                                                                                                                              aTCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGCGCCTTTTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGa                                                                                                                                        >  1:359940/1‑181 (MQ=255)                                                                                                                                                                              aTCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGa                                                                                                                                        <  2:359940/181‑1 (MQ=255)                                                                                                                                                                                    ccctggccctgACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑A‑AACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG‑CATCg                                                         <  1:207745/251‑1 (MQ=255)                                                                                                                                                                                                              aggCACGCCTTCTTCAATAGAGGCGGTACCCAATACTTTGGTATCGCCAAATTCGACCAGAACCGGGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑ACAACAA‑TGTGGCTGCGCATTTTACGGACTTCCGGCGGGTATTCCTATCCTGACAAGGACA‑CGATGGCTATAATCCTTTCACCcctcct                               >  2:358764/2‑251 (MQ=255)                                                                                                                                                                                                               ggCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑A‑AACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG‑CATCGATGGCTATAATCCTTCCACCTCTCCtt                              <  1:352252/251‑1 (MQ=255)                                                                                                                                                                                                                           tcAATAGAGGCGGTACACAATTCTTTGGTATCGCCAAATTCGGCCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑A‑AACAAATGTGGCTGCGCATTATACGGACTGCCGGCGGTTATTCCTATCCAGACAAGG‑CATCGATGGCTATAATCCTTCCACCTCTCCTTTTGTAAACAg                    <  2:261153/249‑1 (MQ=255)                                                                                                                                                                                                                                 gagGCGGTACACAATGCGTTGGTAGCGCCAAAGTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑A‑AACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTAGCCTGACAAGG‑CATCGATGGCTATAATCCGTCCACCTCTCCTTGTATAAACAGGAACGTCt            <  2:359794/251‑1 (MQ=255)                                                                                                                                                                                                                                      ggTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑A‑AACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG‑CATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGAt       <  1:416108/251‑1 (MQ=255)                                                                                                                                                                                                                                         acacAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑A‑AACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG‑CATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCg    >  1:35240/1‑251 (MQ=255)                                                                                                                                                                                                                                           acaataCTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCGTGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑A‑AACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG‑CATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCa  <  1:193784/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                          CGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCGGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑A‑AACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG‑CATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCA  >  NC_000913/3816201‑3816685

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