Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,347,198	T→C	E43E (GAA→GAG)	fumB ←	fumarase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,347,198	0	T	C	95.9%	169.7 / ‑6.0	49	E43E (GAA→GAG)	fumB	fumarase B
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (24/23);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATC  >  NC_000913/4346953‑4347434                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             ccgccgGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAA‑GAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACt                                                                                                                                                                                                                                        >  1:276374/1‑251 (MQ=255)           acacGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATg                                                                                                                                                                                                                               >  2:157719/1‑251 (MQ=255)            cacGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGGGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATgg                                                                                                                                                                                                                              >  2:366017/1‑251 (MQ=255)                 gcTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCg                                                                                                                                                                                                                         <  1:90075/251‑1 (MQ=255)                 gcTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCg                                                                                                                                                                                                                         >  1:378416/1‑251 (MQ=255)                  cTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCGCCACTTTCAGGATGGTTTCGc                                                                                                                                                                                                                        >  2:242945/1‑251 (MQ=255)                              aCCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAAGAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCg                                                                                                                                                                                                             <  2:422530/250‑1 (MQ=255)                               ccGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGc                                                                                                                                                                                                           <  2:211393/251‑1 (MQ=255)                                 gacgaTGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCaa                                                                                                                                                                                                         <  1:161439/251‑1 (MQ=255)                                  acgaTGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAAc                                                                                                                                                                                                        >  2:148445/1‑251 (MQ=255)                                  acgaTGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAAc                                                                                                                                                                                                        >  2:250379/1‑251 (MQ=255)                                       gatCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTaa                                                                                                                                                                                                   <  1:315784/251‑1 (MQ=255)                                                                    gcagcaCGCCTTTGGCGTCGATTTCGGCGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTTGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTTTTGTGCCGTGCGGAGCATAAAAGAAGCGTCGTGAGAGGCTTGCTGCGCCAGCAGGGTCAGGTCTTCTGGCTCCACTTTCATGATGGTTTCGGCTTCGAAGTCGGTAACGCTAACGTAAGCGGATGTGAGTatatttgattcg                                                                                                                                                                      <  2:66101/251‑12 (MQ=255)                                                                        caCGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCgggat                                                                                                                                                                  >  2:319697/1‑248 (MQ=255)                                                                               ttGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCtt                                                                                                                                                            <  1:392710/250‑1 (MQ=255)                                                                                tGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCttt                                                                                                                                                           <  1:250379/250‑1 (MQ=255)                                                                                                gTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGTAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTg                                                                                                                                          <  2:303026/251‑1 (MQ=255)                                                                                                  ttCTTAAGAATTGCAGCGCCACGTGCTTGTCGGTTTCGCTGGCTTCTGGATCTTGAAGAAGAGCCGCAACCTGTTTCTGGTGGGCCGGTCGGAGCATAAAAGATGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCGCTTTCAGGATGTTTTGGCCGGCGAAGTCGGCAACGCTAACGTAAACGGAAGTGAGTAGATAGTATTCGGTATTGTCTGTCCCCATGGTGAATGGTGc                                                                                                                                         <  2:217261/250‑1 (MQ=255)                                                                                                  ttCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGc                                                                                                                                         <  2:298773/250‑1 (MQ=255)                                                                                                          aaTTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGATAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGt                                                                                                                                    >  1:440195/1‑247 (MQ=255)                                                                                                           aTTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATAGTGAAGAATAGACGCAACCTGTTTCTGGCGTGCCAGCCGGCGCATAAAAGAAGCGTGGTGAACGGCTTGCTGCGCCAGCAGGGACAGGGCTTCTGGCGCCACTTTCAGGATGGTTAAGACGTCGAAGTCGGCAACGCTAACGTAATCGGACGTGAGTAGATGGAATTCGGCATTGTATTTACGCATCGGGAAAGGTGCCAGGTANAt                                                                                                                                >  2:421984/1‑250 (MQ=255)                                                                                                             tGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATaa                                                                                                                              <  1:157719/250‑1 (MQ=255)                                                                                                              gCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGACGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATaa                                                                                                                              >  2:35563/1‑249 (MQ=255)                                                                                                                        gTACTTGTCGTTTTCGCTGGCGTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATTGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGtttgtttg                                                                                                                  <  1:436582/251‑1 (MQ=255)                                                                                                                         aaCTTGTCGTTTTCGCTGGCTTCTGGGTCGTGAAGAATAGCCCCAACCTGATTCTGGTGTGCCGGGCGGAGCATAAAAGAAGAGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCCACTGACTGCACTTTCAGGCTGGATTCGCCGTCGAAGTCGGCAACGCTAACGAAATCGGTAGTGAGTAGAGAGTAGTCGGTCTTGTCCTTCCCCATCGGGAACGGCACCTGGGCGATAAAGGGtttgtttg                                                                                                                  >  2:105917/2‑250 (MQ=255)                                                                                                                          aCTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGa                                                                                                                 <  2:7433/250‑1 (MQ=255)                                                                                                                             tGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTGGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTGTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCCGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACAt                                                                                                              <  1:318827/250‑1 (MQ=255)                                                                                                                              gTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACAt                                                                                                              >  1:323414/1‑249 (MQ=255)                                                                                                                                cGTTGTCGCTGGCTTCTGGATCGTGGAGAATAGCCGCAACCGGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCt                                                                                                          <  1:153756/251‑1 (MQ=255)                                                                                                                                 gTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGGGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCtt                                                                                                         >  2:49508/1‑251 (MQ=255)                                                                                                                                       gCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAg                                                                                                                                                                                                                <  1:42722/142‑1 (MQ=255)                                                                                                                                       gCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAg                                                                                                                                                                                                                >  2:42722/1‑142 (MQ=255)                                                                                                                                           gCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGCGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTa                                                                                               <  1:242945/251‑1 (MQ=255)                                                                                                                                                 ggATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTg                                                                                         <  2:378416/251‑1 (MQ=255)                                                                                                                                                          agaaTAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGc                                                                                                           >  1:352601/1‑224 (MQ=255)                                                                                                                                                          agaaTAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGc                                                                                                           <  2:352601/224‑1 (MQ=255)                                                                                                                                                            aaTAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTAtt                                                                              <  1:228196/251‑1 (MQ=255)                                                                                                                                                                 ccGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGGAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGaa                                                                          >  1:240919/1‑250 (MQ=255)                                                                                                                                                                       ccTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCg                                                                   >  1:134590/1‑251 (MQ=255)                                                                                                                                                                               tGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGAGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCNCGt                                                           >  2:421873/1‑251 (MQ=255)                                                                                                                                                                                         ggCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGggc                                                 >  1:41857/1‑251 (MQ=255)                                                                                                                                                                                          gCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGggcg                                                <  2:134590/251‑1 (MQ=255)                                                                                                                                                                                                   aaaaGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCaa                                        >  1:223317/1‑250 (MQ=255)                                                                                                                                                                                                    aaaGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCCCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCaa                                        >  2:381140/1‑249 (MQ=255)                                                                                                                                                                                                          gCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCGCCACTTTCAGGAGGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCTCGTAGTCACGGGCGGCACGCAAAGTGCAtt                                <  2:301267/251‑1 (MQ=255)                                                                                                                                                                                                                    aGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGGGCATTTATAAGAAcc                      >  2:259471/1‑251 (MQ=255)                                                                                                                                                                                                                               gcCAGCAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATcgcg           <  1:29282/251‑1 (MQ=255)                                                                                                                                                                                                                               gcCAGAAGGGTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAAAGCTAACGTAATCGGAAGTGCGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGACTGGTAGCTAACGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTCTAAGAACCCGTACATcgcg           >  2:221460/1‑251 (MQ=255)                                                                                                                                                                                                                                        gTCAGGGCTTCTGGCTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATc  <  1:300958/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             CCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATC  >  NC_000913/4346953‑4347434

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