Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	276,899	0	G	A	50.8%	‑5.2 / 81.5	62	G62E (GGG→GAG)	mmuM	homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (15/15);  new base A (17/14);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 8.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTTGTGCTATGGTGCTTATTTCCTTACTCAACCCCGAAACGCAAAACA‑GGAGCCAGAACATGTCGCAGAATAATCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGATGGCTCTGAATACCGTGGCGATTATCATTGTAG  >  NC_000913/276656‑277139                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 gTTGTGCTATGGTGCTTATTTCCTTACTCAACCCCGAAACGCAAAACA‑GGAGCCAGAACATGTCGCAGAATAATCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCa                                                                                                                                                                                                                                            >  1:124262/1‑250 (MQ=255)  ttGTGCTATGGTGCTTATTTCCTTACTCAACCCCGAAACGCAAAACA‑GGAGCCAGAACATGTCGCAGAATAATCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGACGGGCATCTTGATTACTACCGGGCGGAGgccca                                                                                                                                                                                                                                            >  2:73132/1‑246 (MQ=255)  ttGTGCTATGGTGCTTATTTCCTTACTCAACCCCGAAACGCAAAACA‑GGAGCCAGAACATGTCGCAGAATAATCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAAt                                                                                                                                                                                                                                          >  2:64084/1‑251 (MQ=255)   tgtgCTATGGTGCTTATTTCCTTACTCAACCCCGAAACGCAAAACA‑GGAGCCAGAACATGTCGCAGAATAATCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAAt                                                                                                                                                                                                                                          <  1:246998/250‑1 (MQ=255)    acgCTATGGTGATTATTTCCATACTCAACCCCTAAACACACAACA‑GGAGCCAGAACATGTAGCAGAATAATCCGTTACGCGATCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGAGCGAGAGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAATAGCATCTTGATTAATACCGGGCGGGGGCGCACTgc                                                                                                                                                                                                                                        >  2:177679/3‑251 (MQ=255)           ggTGCTTATTTCCTTACTCAACCCCGAAACG‑AAAACAAGGAGCCCGAACATGTCGCACAATAATCCGTTACGCGCTCTTCTTGCTAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACAGAGCTGGCAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATc                                                                                                                                                                                                                                  >  1:179040/1‑250 (MQ=255)              gCTTATTTCCTTACTCAACCCCGAAACGCAAAACA‑GGAGCAAGAACATGTCGCAGAATAATCCGTTAAGCGCTCTTCTTGATAAACAGGATATACTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGAGCAATGCGCGATCACTg                                                                                                                                                                                                                              >  1:23293/1‑251 (MQ=255)               cTTATTTCCTTACTCAACCCCGAAACGCAAAACA‑GGAGGCAGACCATGTCGAAGAATAATCCGTTACGCGCTCTTATTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGATTATCCGCGAAGTGAATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGc                                                                                                                                                                                                                             >  1:14197/1‑251 (MQ=255)                 tATTTCCTTACTCAACCCCGAAACGCAAAACA‑GGAGCCAGAACATGTCGCAGAATAATCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGTTGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGcc                                                                                                                                                                                                                            >  1:151717/1‑250 (MQ=255)                  aTTTCCTTACTCAACCCCGAAACGCAAAACA‑GGAGCCAGAACATGTCGCAGAATAATCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGTAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGcc                                                                                                                                                                                                                            >  2:79621/1‑249 (MQ=255)                    ttCCTTACTCAACACCGAAAAGCAAAACA‑GGAGCCAGAACATGTCGCAGAATAATCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGGAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCt                                                                                                                                                                                                                        >  2:286181/1‑251 (MQ=255)                         tACTCAACCCCGAAACGCAAAACA‑GGAGCCAGAACATGTCGCAGAATAATCCGTTACGCGCTCTTCTGGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAATCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAAGCCGGAGCTTATCCGCTAAGTGCACCTTGATTACTACCGGCCTTAGTCGCAATGCGCGAACACTGCCAGCTATc                                                                                                                                                                                                                     <  2:194814/249‑1 (MQ=255)                                   cGAAACGCAACACA‑GGAGCCAGAACATGTCGCAGAATACTCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACACGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCcgg                                                                                                                                                                                                         >  2:120175/1‑251 (MQ=255)                                      aaCGCAAAAAA‑GGAGCCAGAACATGTCGCAGAATAATCCGTAACGCGCTCTTCTTGCTAAACAGGATATCCAGCTGATGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGAAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGCGGCGCAATGCGCGCTCACTGCCAGCTATCAGGCGACGCcggcgg                                                                                                                                                                                                      >  2:32311/1‑251 (MQ=255)                                       aCGCAAAACA‑GGAGCCAGAACATGTCGCAGAATAATCCGCTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCGGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGTGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCctgcgtg                                                                                                                                                                                                     <  2:72496/251‑3 (MQ=255)                                           aaaaCA‑GGAGCCAGAACATGTCGCAGAATAATCCGTTACGCGCTCTTCTTGATAAACAGGATATCGTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTc                                                                                                                                                                                                 <  1:294872/251‑1 (MQ=255)                                                      ccAGAACATGTCGCAGAATAATCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCACGATCACTGCCAGCTATCAGGCGACGCCGGCGTGCTGCGCCGCGCGAg                                                                                                                                                                                       <  1:219475/251‑1 (MQ=255)                                                       cAGAACATGTCGCAGAATAATCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAgg                                                                                                                                                                                      <  1:161842/251‑1 (MQ=255)                                                         gAACATGTCGCAGAATAATCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTc                                                                                                                                                                                    >  1:25002/1‑251 (MQ=255)                                                           aCATGTCGCAGAATAATCAGTTACGCGCTCTTCTTGATAAACAGGATCTCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGACAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCCGCTATCAGGCGACGCCGGCGGGCTTCGCCCCCCGAGGGCTg                                                                                                                                                                                  >  2:278210/1‑251 (MQ=255)                                                               gTCGCAGAATAATCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCCTCGCCGCCCGCGGTCTGGATg                                                                                                                                                                              >  2:249059/1‑251 (MQ=255)                                                                   cAGAATAATCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAgc                                                                                                                                                                          <  2:289181/251‑1 (MQ=255)                                                                    aGAATAATCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGTTGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGCCAGCCTGAGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAgcg                                                                                                                                                                         >  1:305045/1‑251 (MQ=255)                                                                         aatCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTc                                                                                                                                                                    <  1:60464/251‑1 (MQ=255)                                                                          atCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGCGGTCTGGATGAAGCGCAGTCg                                                                                                                                                                   >  1:302300/1‑251 (MQ=255)                                                                          atCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCg                                                                                                                                                                   <  1:232884/251‑1 (MQ=255)                                                                          atCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCCTCGCCGCGCGAGGTCTGGATGAAGCGCAggcg                                                                                                                                                                   >  2:193442/1‑248 (MQ=255)                                                                              gTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAg                                                                                                                                                               <  1:26575/251‑1 (MQ=255)                                                                                   gcgcTCTTCTATATAAACAGGATGGCGTGCTGTTTGATGGCGCGATGGCGACGTGGCTGGACGCGGGAGGGTGTAACTTAGCCGACAGCCTGGGGGCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAACGCGCtt                                                                                                                                                          <  2:100667/251‑2 (MQ=255)                                                                                           cttGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCa                                                                                                                                                   >  1:257855/1‑250 (MQ=255)                                                                                                   aCAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTgg                                                                                                                                          <  1:80402/251‑1 (MQ=255)                                                                                                    cAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGGCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCCAAAGCGTggg                                                                                                                                         >  2:203760/1‑250 (MQ=255)                                                                                                        atatCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTg                                                                                                                                     <  2:133099/251‑1 (MQ=255)                                                                                                              ggctgctgGAGGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTGAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGAGCGt                                                                                                                               <  2:126156/250‑1 (MQ=255)                                                                                                                 tgctgGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGCGGGCTGGATGAAGCGCCGTAGGAAGCGCTGGTGGGCAAAAGCGGTGTGCTGGcgcggaa                                                                                                                             >  2:68670/1‑247 (MQ=255)                                                                                                                       aTGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGGGGAGCTGGCGCGTAAAGCCCgg                                                                                                                      >  2:105646/1‑250 (MQ=255)                                                                                                                         ggCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATGACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCGGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTTGCGCGTAGAGCCCGTGa                                                                                                                    <  2:172373/251‑1 (MQ=255)                                                                                                                           cgcgATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGa                                                                                                                    >  1:247356/1‑249 (MQ=255)                                                                                                                            gcgATGGCGACGGGGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGc                                                                                                                 <  2:272086/251‑1 (MQ=255)                                                                                                                                    gACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCCGGTTGGCAAAAGCGGGGAGCGGGCGCGTAAAGCCCGTGACGCGCATCTgg                                                                                                         >  2:300207/1‑251 (MQ=255)                                                                                                                                          gCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCag                                                                                                      >  1:121616/1‑248 (MQ=255)                                                                                                                                                  ggcgAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCCCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAg                                                                                           <  2:293704/250‑1 (MQ=255)                                                                                                                                                    cgAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGc                                                                                         <  1:135335/251‑1 (MQ=255)                                                                                                                                                      aGGGTGTCACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAACCCCGGAGCTTATCCGCTAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCACAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAggagg                                                                                       <  2:216631/251‑4 (MQ=255)                                                                                                                                                       gggTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCggg                                                                                      <  1:122546/251‑1 (MQ=255)                                                                                                                                                              aCTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCggggaagctt                                                                               >  1:178029/1‑244 (MQ=255)                                                                                                                                                                 gaGCCGACATCGTGTGGTAGGCCGAAGTGCTGGTAGAAATCCCGGAGCTTAGCCGCGAAGTGCATATTGATATCTACCGGGGGGAGGCGCAAGGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCTCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGGGAAGCGTGGAGGGCAGAGAACCCGCAGGCGGGAAccctgctt                                                                            <  2:157644/250‑8 (MQ=255)                                                                                                                                                                         aGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGCTGGAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGGGGAGCTGGCGCGGAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCCCCGGCggggaaccctctggggggggg                                                                    >  2:61256/1‑233 (MQ=255)                                                                                                                                                                                ggTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTTGCGGGATCCGTc                                                             <  1:248799/251‑1 (MQ=255)                                                                                                                                                                                ggTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAACGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCGGGGACGCCTCtggggtggggggccccg                                                             >  2:135122/1‑237 (MQ=255)                                                                                                                                                                                          aGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGGCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCGGGGACCCCTCTGGGGGCgggggcccgcggggccttcg                                                   >  1:20210/1‑234 (MQ=255)                                                                                                                                                                                                   aGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCCGTCGAAAGCGCt                                                                                                                                                           >  1:285880/1‑138 (MQ=255)                                                                                                                                                                                                   aGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCt                                                                                                                                                           <  2:285880/138‑1 (MQ=255)                                                                                                                                                                                                    gAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCg                                                                      <  2:271223/222‑1 (MQ=255)                                                                                                                                                                                                       aaCCCGGAGCTAATCCGCGAAGTGCATCTTGATTACTACCGGGCGTGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGc                                      <  2:125825/251‑1 (MQ=255)                                                                                                                                                                                                                cTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTTGATGAAGCGCAGTCGAAAGCGCTTATTGGCCAAAGCGTGGAGCTGGCTCGTAAAGCCCGTTAAGCGTATCTGGCAGAGAACCCGCCGGCGGGAACCCTTCTGGTTGCCGGGTCCCTCGGGCCCTCCGgcgccgctcttgccgcaggccc                             >  1:101267/1‑233 (MQ=255)                                                                                                                                                                                                                  tATCCGCGAAGTTCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCGGGCAGAGAACCCGCAGGCGGGAACGCGTCTGGTGGCGGGAGCCGTCGGGCCTTACGGCGCGTATCTGGCGGATGGctct                            <  2:25002/250‑1 (MQ=255)                                                                                                                                                                                                                   aTCCGCGAAGTGCATCTTGATTACTACCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGATGGCTCTGa                          <  1:117926/251‑1 (MQ=255)                                                                                                                                                                                                                                  cTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCGGGGAAGCCTCTGGGGGCGGGGGCCGGCGGGCCTTTCGGCGCGGATCTGGCggggggggctggataccggggggggtt           >  1:146057/1‑226 (MQ=255)                                                                                                                                                                                                                                      aTTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGGCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCGGGGACGCCTCTTGGGGCGGGGACCGTCGGGCCTTTCGGGGCGTATTTGGCGGGTGGGCCTGGATACCCTGGCggtttatat       >  1:24201/1‑243 (MQ=255)                                                                                                                                                                                                                                        tactacCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGATGGCTCTGAATACCGTGGCGATTATCAtt      <  1:163297/250‑1 (MQ=255)                                                                                                                                                                                                                                            acCGGGCGGAGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGATGGCTCTGAATACCGTGGCGATTATCATTGTAg  <  1:193442/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 GTTGTGCTATGGTGCTTATTTCCTTACTCAACCCCGAAACGCAAAACA‑GGAGCCAGAACATGTCGCAGAATAATCCGTTACGCGCTCTTCTTGATAAACAGGATATCCTGCTGCTGGATGGCGCGATGGCGACGGAGCTGGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGATGGCTCTGAATACCGTGGCGATTATCATTGTAG  >  NC_000913/276656‑277139

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: ‑