Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	23,032	C→T	F214F (TTC→TTT)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	23,032	0	C	T	98.3%	188.2 / NA	59	F214F (TTC→TTT)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (37/21);  minor base G (0/1);  total (37/22)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCAC  >  NC_000913/22788‑23274                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   aGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGcg                                                                                                                                                                                                                                              >  1:88010/1‑250 (MQ=255)   gTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGcggt                                                                                                                                                                                                                                            <  1:305248/251‑4 (MQ=255)     tGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTc                                                                                                                                                                                                                                           <  1:46954/250‑1 (MQ=255)       aCGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGAt                                                                                                                                                                                                                                        <  1:34155/251‑1 (MQ=255)                 cAAAGACTTTATCCGGCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCAAACATCATCCGCGCGCAGGGCAAACTAATCGGCACCGGTCACCTGCCCAAAGGCGCGAAGCCGGTACACTGGTGCGTTGCCTGCCGTTCTACGCTGGCGGAAGCGGACGTTGAGTAATAAGACAAAACTTCTCCGTCCATCGACGTAGCTTTTCAGGCGGACGATCAGGaagcc                                                                                                                                                                                                                               >  2:352058/1‑246 (MQ=255)                     gccTTTAGCCGTCTGGGCGTGCGGGGCGACTGGTCGCACCCGTACCTGACCATGGACGTCAAAACTGAAGCCAACATCAGCCGCGCGCTGGGCAAAATCAGCGGCAACGGTCACCTGCACAAAGGGGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTGCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTGTCAGGCGGTCGATCAGGATGCACTGaa                                                                                                                                                                                                                          <  2:196439/249‑1 (MQ=255)                      aCTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGACCAGGATGCACTGaa                                                                                                                                                                                                                          >  1:264489/1‑250 (MQ=255)                              cGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAAtt                                                                                                                                                                                                                 >  1:209806/1‑251 (MQ=255)                                      cGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTa                                                                                                                                                                                                          >  1:346915/1‑250 (MQ=255)                                                    ggtccaaccctgaccataCCATGGTCTTCACACCTGAAGCCAACATCCTCCGCGCGCTGGGCAAAATCATCGGCAACGGGCACCTGCACAAATGCGCGAAGCCGGTGCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAATCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCACCGACGTTGCGTTTCAGGCGGTCGAGCAGGATGCACTGAAAGAAAAATTTGCCGTAAGCAACGTTAACg                                                                                                                                                                                             <  2:202698/232‑1 (MQ=255)                                                       cGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCa                                                                                                                                                                                        <  1:309558/251‑1 (MQ=255)                                                        gCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCa                                                                                                                                                                                        >  1:79114/1‑250 (MQ=255)                                                                 ccTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCGGGGCAAAATCATCGGCAACGGGCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTgg                                                                                                                                                                              <  2:88010/251‑1 (MQ=255)                                                                     aCCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGACCCTCGACGTTGCTTTTCAGACGGTCGATCAGGATGCACTGAAAGACAAATTTGCCGTAAGCAACGTTAACGGCCACATCTCGCTGGtgat                                                                                                                                                                          >  2:31070/1‑248 (MQ=255)                                                                         tGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTCCGACAAAACTTCTCCGTCCATCGACGCTGCTTTTCAGGCGGTCGATCAGGATGCCCTGAAAGCACAATTTGCCGTAATCACCGTTAACGGCCAAATCTCGCTGGTAAttagg                                                                                                                                                                      >  2:154179/1‑247 (MQ=255)                                                                            aCTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGCAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTTGacc                                                                                                                                                                   >  1:199055/1‑251 (MQ=255)                                                                              ttCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGACGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCCCTGCAAACAAAATTTGCAGTAAGCAACATTAACGGCCCATTCTCGcgggggatctgggcccc                                                                                                                                                                 >  2:256491/1‑235 (MQ=255)                                                                                        aaGCCAACAACATCCGCGCGCTGGGCAAAATCAGCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTgg                                                                                                                                                       <  2:484539/251‑1 (MQ=255)                                                                                         aGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTgg                                                                                                                                                       >  1:5974/1‑250 (MQ=255)                                                                                          gCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCAGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAAAGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGAc                                                                                                                                                     >  1:12686/1‑251 (MQ=255)                                                                                             atcaccatcCGCGCGCTGGGCAAAATCCTCGGAAACGGTAACCTGCAAAAAGGCGAGAAGCCGGTAAATGGGTGCGTTGACTGCCGTGCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGActct                                                                                                                                                  <  2:422941/246‑1 (MQ=255)                                                                                             aaCATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGCCAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGActct                                                                                                                                                  >  1:298042/1‑251 (MQ=255)                                                                                                atcatcCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTcagcc                                                                                                                                               >  1:242823/1‑251 (MQ=255)                                                                                                        cgcgcTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTctgcctgcgaccc                                                                                                                                       >  1:103800/1‑246 (MQ=255)                                                                                                             tGGGCAAAATCATCGGCAACGGTCTCCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCgcgc                                                                                                                                   >  1:123181/1‑250 (MQ=255)                                                                                                                    aaTCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTa                                                                                                                           <  1:239649/251‑1 (MQ=255)                                                                                                                     atcatcGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTAt                                                                                                                          <  1:227793/251‑1 (MQ=255)                                                                                                                        atcGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGc                                                                                                                       >  1:143641/1‑251 (MQ=255)                                                                                                                          cGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCAc                                                                                                                     >  1:296228/1‑251 (MQ=255)                                                                                                                            gCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCAcct                                                                                                                   <  1:394655/251‑2 (MQ=255)                                                                                                                                 ggTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGAttt                                                                                                              >  1:117997/1‑251 (MQ=255)                                                                                                                                 ggTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGGTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGAttt                                                                                                              >  2:98103/1‑251 (MQ=255)                                                                                                                                            caAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTctgcc                                                                                                                                               <  2:12921/207‑1 (MQ=255)                                                                                                                                            caAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTctgcc                                                                                                                                               >  1:12921/1‑207 (MQ=255)                                                                                                                                                 gcgcGAAGACGGTACACTGGTGCGTTGACTGCCGTCCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAACACAAAATTTGCCGTAAGCAACGTTAAAGGCCCAGTCTCGCTGGTGATCTGTCCCCCTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCAATTTCCCCTGCTTTCGACTATGCGCTGGTg                                                                                              >  2:32458/1‑251 (MQ=255)                                                                                                                                                    cgAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGCCTATGCGCTGGTGc                                                                                             >  1:89218/1‑249 (MQ=255)                                                                                                                                                        gccggtacACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTAGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGc                                                                                             >  1:233736/8‑245 (MQ=255)                                                                                                                                                           ggtacACTGGTGCTTTGACTGCCGTTCTGCGCTTGCGGAAGCGGGCGTTGCGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGCTCAGGCTTCACTGAAAGCAAAATTTACCGTAAGCAACGTTAACGGCCCCCTCTCGCTGTTGATCTATACCCCTACTCCGTTGACTCTGCCTCCCAACCGCTCACTCTCTTTTTCCCCTTTTTTCGTCTATGCGCTGGTGCAggtccacc                                                                                    >  2:35812/5‑244 (MQ=255)                                                                                                                                                                              gtcggTTCTGCGCTTGCGGAAGCGGAAGTTGAGTATTACGACAAAACCTCTCCGTCCATCGACGTGGAGTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTCGCCGTAAGCAACGTGAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAAGCTCTACTGCACCTGATTTCGTCTAGGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTgga                                                                 <  2:166884/247‑2 (MQ=255)                                                                                                                                                                                            ggCGTAAGCGGAAGTTGTGTATTACGACAAAACTTCTCCGTCCAGCGACGTTGCTTTTCAGGCGGTCGATCATGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGGTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCTTGTACTCTGCCTGCGAACCGCGACATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCATATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGtt                                                    <  2:407693/250‑1 (MQ=255)                                                                                                                                                                                               ggaagcgtaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCCCGCTGGTGATCTGGACCACTACGCCGTGGAGTCTGACTGCGAAACGGCCAATCTCTATTGCAACTGATTTCGACTACGCGCTGGTGCAGGTCGACGGTCAGGCCGTGATTTTGGCGAAAgctagagttgac                                                 >  2:153141/1‑239 (MQ=255)                                                                                                                                                                                                   gcggaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGt                                            >  1:269122/1‑251 (MQ=255)                                                                                                                                                                                                         gTTGTGTATTACGACAAAACTTCCCCGTCCAACGACGTTGCTTTTCAGGCGGTAGATCAGGATGCACTGAAAGCAAAATTCGCCGTAAGCAACGTTAACGGACCAATCTCGCAGGGGATCTGGACCACTACGCCGTGGACTCTGCCTGCGACCCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGAGCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTTGTCGAAAGCGTAATGCa                                      <  2:106332/251‑1 (MQ=255)                                                                                                                                                                                                          ttGAGTATTACGCCAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGCCGGTCCGGCCGTGATTCTTGCGAAAAATCTGTTTGAAAGCGTAAtccgg                                     >  2:355536/1‑247 (MQ=255)                                                                                                                                                                                                           tGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGCCGGTCAGGCCGTGTCTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGc                                    >  2:206350/1‑251 (MQ=255)                                                                                                                                                                                                                   ccGACAAAACTTCGCCGTCCAGCGTCGTTGCTTTGCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTTACGGCCCAATCGCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGc                            <  2:221645/250‑1 (MQ=255)                                                                                                                                                                                                                   aCGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCACCGGTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCgg                             <  2:273930/250‑1 (MQ=255)                                                                                                                                                                                                                       cAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTgg                        >  2:151927/1‑250 (MQ=255)                                                                                                                                                                                                                        aaaaCTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAc                       >  1:189799/1‑251 (MQ=255)                                                                                                                                                                                                                        aaaaCTGCACCGTCCAGCGACGTTGCTTTTCAGGCGGTCTATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGGTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAc                       <  2:23536/251‑1 (MQ=255)                                                                                                                                                                                                                         aaaCTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGGTCTGGTTGAAAGCGTAATGCAGCGTGTCGGCGTGAcc                      >  2:127087/1‑251 (MQ=255)                                                                                                                                                                                                                         aaaCTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAcc                      <  1:6549/251‑1 (MQ=255)                                                                                                                                                                                                                            cTTCTCCGTACATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCACAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACAGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGCTCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCggt                   >  2:436398/1‑249 (MQ=255)                                                                                                                                                                                                                                tcCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGCTTac                >  1:260096/1‑250 (MQ=255)                                                                                                                                                                                                                                  cGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGCTGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGAAGCGTATCGGCGTGACCGATTACAcc             >  2:150008/1‑251 (MQ=255)                                                                                                                                                                                                                                          gACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtctcag     <  2:346915/251‑3 (MQ=255)                                                                                                                                                                                                                                            cGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtctcagcc   >  1:5520/1‑247 (MQ=255)                                                                                                                                                                                                                                            cGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCa   <  1:366276/251‑1 (MQ=255)                                                                                                                                                                                                                                             gTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCCGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtctcaaccc  >  2:266037/1‑246 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   AGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCAC  >  NC_000913/22788‑23274

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