Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,092,703	A→G	F41L (TTC→CTC)	frvA ←	putative PTS enzyme IIA component FrvA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,092,703	0	A	G	97.7%	131.9 / NA	44	F41L (TTC→CTC)	frvA	putative PTS enzyme IIA component FrvA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (18/25);  minor base T (0/1);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAATAATTTTGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGATTTGCGGAGTAAT  >  NC_000913/4092456‑4092948                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       ggaagtattgtGCGACACATGGGGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCAGCGCTTGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGTCGCGTGTACCGTTCGGCACGGCGACACCGTAACGAAATCCAGTAGAGTGCATTTTTTCGCGCAGCGGCAGGGTTGGCAGGAGCTg                                                                                                                                                                                                                                                     <  2:210677/241‑1 (MQ=255)   aataatTTTGCGACACATTGTGCCGATGATTTTGCCCTGATCCTCTTCGCCTCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACTGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTGCGCGCAGCAGCAGGGTTTGCAGGAGCTGgt                                                                                                                                                                                                                                                   <  1:376871/251‑1 (MQ=255)   aataatTTTGCGACACAGTCTGCCGATGATTGTGACCTTATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGgt                                                                                                                                                                                                                                                   <  2:227637/251‑1 (MQ=255)     aaatTTTTCGACACCGGGTGCCGATGATTTTGACCTGATCCCCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGGCTTCGCCATCGCTGGCTGTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTGGGTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACTCCGGAACCAAATCCAGTAGAGTGCATTTTGTCGCGCATCAGCAGGGTTTGCAGGAGCTGgtgt                                                                                                                                                                                                                                                 <  2:274027/250‑1 (MQ=255)                 acaGTGTGCCGATGAGTTTGACCTGATCCTCTTCGCCGCTTTGCGTCACGCCGAGGCAGAGCCAGCAATTGACGGCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATa                                                                                                                                                                                                                                     <  2:99379/251‑1 (MQ=255)                   aGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAACTCCAGTAGAGTGCATTTTTTCGCGCAGCAGCCGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATaa                                                                                                                                                                                                                                    >  2:108012/1‑250 (MQ=255)                      gtgCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAAcc                                                                                                                                                                                                                                <  1:391515/251‑1 (MQ=255)                           gatgatTTTGCCCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGtt‑tt                                                                                                                                                                                                                          <  1:321090/251‑1 (MQ=255)                              gatTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGt                                                                                                                                                                                                                        >  2:293478/1‑250 (MQ=255)                              gatTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTa                                                                                                                                                                                                                       <  1:129086/251‑1 (MQ=255)                                         gATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTc                                                                                                                                                                                                            <  2:72439/251‑1 (MQ=255)                                          aTCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCGTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCg                                                                                                                                                                                                           <  1:410770/251‑1 (MQ=255)                                               ctTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAAc                                                                                                                                                                                                      >  1:274910/1‑251 (MQ=255)                                                          tGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATGTTCTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGaa                                                                                                                                                                                            <  1:361364/250‑1 (MQ=255)                                                                aCGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAAt                                                                                                                                                                                      >  1:131815/1‑250 (MQ=255)                                                                      aGGCAAATCCAGCAATTAACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGACTAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGCCACAGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAATCTGTGATAAAACCGTT‑TTGTATCGCTATTATCTCCAACTGCTTCAGAACGTAATAAGCGcc                                                                                                                                                                               >  2:341464/1‑251 (MQ=255)                                                                        gCAGATCCAGCAATTGACGTCTTCGCCATCGCTAGCTTTCCAGTCAATCGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAAGCGCTTTTACCGAGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCAt                                                                                                                                                                             >  1:196374/1‑251 (MQ=255)                                                                           gATCCAGCAATTGACGTCTGCGCCATCTCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACTGCGACACCGGAACCAAAGCCAGTAGAGTGCATTGTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGc                                                                                                                                                                          <  2:344706/251‑1 (MQ=255)                                                                              ccAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCt                                                                                                                                                                       <  2:312418/251‑1 (MQ=255)                                                                                    aatgatgTCTTCGCCATCGCTCGCCTTACAGACAATAACATGTGTTTTGAGGGAGAATAATACGCACGGATGTTTCACGCCGGCCCTTTTACGGTGCGGCACGAAAACACCGGAACCAAATCCAGTAGAGTGCATTTTTTAGCGCAGAAGCAGGATTTGCAGGAGCTGGTGTGAGTCTGTGATAAAACATTT‑TTGTAACGCTATTGGCGCCAACTGCTTCAGAACGGAATATACTCCATTGCCGTGAatat                                                                                                                                                                 >  2:428430/7‑251 (MQ=255)                                                                                                      gCTGGCTTTCCAGTTAATAGCCTGCGCTTTGCGGGCTAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATGTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTg                                                                                                                                               <  2:219818/251‑1 (MQ=255)                                                                                                       cTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGc                                                                                                                                              >  1:135800/1‑251 (MQ=255)                                                                                                               ccAGTCAATAGCCTGCGCTATGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCTCCAACTGCGTCAGAACGGAATAAGCGGCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAagag                                                                                                                                      <  2:186379/251‑1 (MQ=255)                                                                                                                      aTAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCAt                                                                                                                               >  2:192768/1‑251 (MQ=255)                                                                                                                             gcgtgttGCGGGCGAATAATACGAACGGTGGTTTAACGCAGGCGCTTTTACCGGGCGGCACGGCGACATCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGGTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCGGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGcctcggctt                                                                                                                          <  2:131815/244‑8 (MQ=255)                                                                                                                             gcgcTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATAGTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATc                                                                                                                          <  1:215513/249‑1 (MQ=255)                                                                                                                                    gCGGGCGAATAATTCTAACGGTTGTTTAACGCAGTCGGTTGTACCGTGCGTCACGGCGACATCGGAACCAAAGCCAGTTGAGTGTATTGGTTCGCGCAGCAGCAGGGTTTGCAGGAGCGGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAGCTGCTTCGGAACGGAAGAAGCTCTATTGCTCTGAATATGAAGGTCAAGATAGCTTGAAGTAGGGTCTGGCAGTAATCGGACTCCa                                                                                                                  <  2:425903/250‑1 (MQ=255)                                                                                                                                            ataatcCGCACGGTTGTTTAACGCAGGCGATTTTACCGTGCGGCACGGAGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGCGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAatat                                                                                                          >  2:430556/1‑250 (MQ=255)                                                                                                                                                                   cAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTGCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGa                                                                                  <  1:390556/251‑1 (MQ=255)                                                                                                                                                                    aGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAAtggtgg                                                                                   >  1:84125/1‑249 (MQ=255)                                                                                                                                                                      gcgcTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGTAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCTATTAATCGGACTCCTTAAAATATATTTTAATCAATTGTAATGGTGGAGCg                                                                               <  2:344625/251‑1 (MQ=255)                                                                                                                                                                        gcTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCCTTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGACATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGAt                                                                             >  2:319165/1‑251 (MQ=255)                                                                                                                                                                                 gTGCGGCACGGCGACCCCGGAACAAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGCACGGAATAAGCGCCATTGCCCTGACTATTCAGGTCACTACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTACTGGTGGAGCGATTAAGtgatg                                                                    >  2:390945/1‑251 (MQ=255)                                                                                                                                                                                                acCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGAAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCg                                                     <  1:41982/251‑1 (MQ=255)                                                                                                                                                                                                aactgAACCAAATCCAGTATAGTGCATTTTTTTGCGCAGCAGCAGGGTTTGCAGGATGTTGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCATCAGAACAGAATAAGCGCCATTGCCCTTAATATTCAGGTCAATACAGCTTGCAGTAATAGCGGTCGGTAATCGGATTTCATAAAATTGATTGTAATCGACTGGAATGGTGGAGCGATTAAGTGATGAGAGGCGCTATGTCg                                                     <  2:271040/247‑1 (MQ=255)                                                                                                                                                                                                   ggAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGCGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGt                                                   >  2:20200/1‑250 (MQ=255)                                                                                                                                                                                                     aaCCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTaa                                                 >  2:276434/1‑250 (MQ=255)                                                                                                                                                                                                                  agagTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAAttt                                    <  2:376207/250‑1 (MQ=255)                                                                                                                                                                                                                    agTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATACAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACACTTTga                                  >  2:417052/1‑250 (MQ=255)                                                                                                                                                                                                                        cATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGAc                               >  1:213829/1‑249 (MQ=255)                                                                                                                                                                                                                             tttCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTg                         <  1:380451/250‑1 (MQ=255)                                                                                                                                                                                                                                         gcagGATTAGCAGGAGCTGGTGTGAGGCGTTGATAAAACCGTTGTTG‑AACGCTATTGTAGACAACTGCTTCAGATCGTAATAAGCGCCATTCTCCTGAATATTCAGGTAAATACAGCGTGCAGTAAGAGCTGCCATTAATCGGAATCCATAAAATATATTGTAATCAATTGTAATGGTGGAGCGATTAAGTGATTAGAGGCGATAGGTCTCTATACAATGACAACTTGTGCCAGAGTTAGACGTCAGAttt            <  2:140485/251‑1 (MQ=255)                                                                                                                                                                                                                                           agGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGGT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGCTTTGc          >  2:17254/1‑251 (MQ=255)                                                                                                                                                                                                                                                   cAGGAGCTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGATTTGCGGAGTAAt  <  1:293478/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       TGAATAATTTTGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTT‑TTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGATTTGCGGAGTAAT  >  NC_000913/4092456‑4092948

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