Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,684,779	0	A	C	45.6%	38.2 / 12.0	50	G243G (GGT→GGG)	glyA	serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (10/15);  new base C (1/20);  total (11/39)
Fisher's exact test for biased strand distribution p-value = 6.17e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAG‑CGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTG  >  NC_000913/2684557‑2685013                                                                                                                                                                                                                                |                                                                                                                                                                                                                                           gtggtTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAG‑CGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGATCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTAGAGCTCTTCGCTACCACCTTTCGCCCGGGTCa                                                                                                                                                                                                                            >  2:256889/1‑239 (MQ=255) gtggtTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAG‑CGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGGGCATCAACGGACCGCCCCGAACACCAGGGAAAACGGCAGAGTTCAGTATTTTGTAAACCTCCTCCCTTCCACCTTTCACCACTATCa                                                                                                                                                                                                                            >  2:1753693/1‑239 (MQ=255)                   ggAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCACCGCTTTAGCGTTTTTAG‑CGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCg                                                                                                                                                                                                          <  1:1550106/239‑1 (MQ=255)                       aCCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAG‑CGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAg                                                                                                                                                                                                      <  1:186238/239‑1 (MQ=255)                                      cgcTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAG‑CGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGGACAGCCCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGGTTTGTGAggtgt                                                                                                                                                                                       >  2:1293018/1‑235 (MQ=255)                                        tccggggAACTCTTCTACCATCGCTTTAGGTTTTTTAG‑CGCCCTGCTGCTGTTAATTTTTGACCCAGGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGGGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGtggtgg                                                                                                                                                                                     <  2:514467/234‑1 (MQ=255)                                        cTCGAGTAACACTTCTGCCATCGCTTGAGCGTTTTTAG‑CGGCCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAGCGGCAGAGTTCAGTTTTTTGTACAGCTCTCCGATACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGtggtgg                                                                                                                                                                                     <  2:9724/239‑1 (MQ=255)                                          cGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAG‑CGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGgta                                                                                                                                                                                   <  1:296039/239‑1 (MQ=255)                                             ggACCACTTCTACCTCCGCTTTGTCGTTTTTAC‑CGACCAGCTGCGTGTAAGTTTTGACCGCAGCCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTGGTACAGCTCTTCGCTCCCCCCTTCCGCCGGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGgtagt                                                                                                                                                                                 <  2:1286753/238‑1 (MQ=255)                                             gaaacctttCTACGTTCGCTTT‑GCCTGTTTAGTCGCCCGCCTGCTGGGAAGTTTTGAACCCAGGCTCCATCGCTTCTTTCAGACCCACCGATTTCCCGGCGATTACGTGCTTCAACGGACCGCCCTGACCCCCAGGGAAAACGGCCGAGTCAATTTTTTTGTCCAGCTCTTCGCTCCCCCCTTTCGCCGGGATCAGGCCCCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGgtagt                                                                                                                                                                                 <  2:137365/231‑1 (MQ=255)                                                  ctttctcacATCACTTTGTCTTTTTTAG‑CGCCCTGCTGTTGGTTAGTTTTAACGTCAGGCACCATCGCTTCTTTCAGAGCAACCGCTTTCCCGGCGATTCCGTGCATTAACGGACCGCCCTGCCCACCAGGGAAAACGGCAGGGTTCAGTTTTTTGTACAGCTCTTCGCTACCTCCTTTCGCCGGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTaacaac                                                                                                                                                                           <  2:1688094/231‑1 (MQ=255)                                                  ccttttgccACCGCGTTGGCTTTTTTAG‑CGTCCTGCTGCAGATAATTTTCAAACGCAGTCCCCATCGTTTCTTTCAGAGCAACCGCTTTGCCGGCGATTACCTGAATCAACGGCCCCCCCTACCCACCAGGGAAAGCGCCAAAGTTCAGTTTTTTGTCCAGCTCTTCGCTACCACCTTTCTCCAGGATCGGCCCCCCGCGGGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTaacaac                                                                                                                                                                           <  2:409775/232‑1 (MQ=255)                                                    ttatcccATCGCTTTAGGTTTTTTAG‑CGACCTGTTCCTGGTACGTTTTGAACTAAGGCTCCATCGCTTCTTTCAGAGCACCCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCCCCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTGCGCTCCCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGt                                                                                                                                                                         <  2:1649788/234‑1 (MQ=255)                                                    ttCTCCCATCGCTTTAGCGTTTTTAG‑CGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCACCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACg                                                                                                                                                                          <  1:822962/238‑1 (MQ=255)                                                        tccATCGCTTGTGCGTTTTTAG‑CGCCCTGCCTGTGGGAACTTGTGACCTGCGGCTCCATCGCTTCTTACAGAGCAACGGCTTTGCCGGCAATTACGTGCATCAACGGACCGCCCTCCCCACCGGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCGGGGTCGGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGc                                                                                                                                                                     <  2:2101254/238‑1 (MQ=255)                                                                                                            ggCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCCCCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTCCCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTg                                                                                                                  <  1:1426417/239‑1 (MQ=255)                                                                                                              cTCCATCGCTTTTTTCAGAGCAACCGCTTTACCGGGGATTACGTGCATCAACGGACCGCCCTGACCCCCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTTTTCGCTCCCACCTTTCGCCAGGATCAGGCCGCCGCGCGGCCCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTgc                                                                                                                 <  1:1164605/238‑1 (MQ=255)                                                                                                               tCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGCCCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTTCCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTaa                                                                                                                                                                               <  1:1169226/175‑1 (MQ=255)                                                                                                               acaaTCGCTTCTAGCAGAGCACCCGATTGACCGTCGATTACGTGCATCCACGGCCCGCCCTGACCTCCAGGAAAACCGGCGGGGTTAAGTTTTTTGTGCAGCTCTGCGCTCCCCCCTTCCGCGAGGATAGGGCCGCCGGGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTgcg                                                                                                                <  2:14351/235‑1 (MQ=255)                                                                                                                    cGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCatat                                                                                                          <  1:1239500/239‑1 (MQ=255)                                                                                                                      cttcttTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGCCCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAgtg                                                                                                                                                                                         <  1:1514691/158‑1 (MQ=255)                                                                                                                      cttcttTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGGACAGCTCTTCGCTACACCCTTTAGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGGGAGTGGAGGGTGGAACAACACGGGGCGGACGAGCCTGGTTCAGGTAGATGCCAGGCGCCAACAGGGCCGGCACGttagcccgaata                                                                                                        >  2:883626/1‑228 (MQ=255)                                                                                                                                         tttCCCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTTTACAGCTCTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAAc                                                                                                                                 <  1:929829/195‑1 (MQ=255)                                                                                                                                         tttACCGGCGATTACGTGCATCAACGGCCCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCGCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAgtg                                                                                                                                                                                         <  1:1098702/139‑1 (MQ=255)                                                                                                                                                 cGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTggtggt                                                                                                                                                                                    <  1:1348153/136‑1 (MQ=255)                                                                                                                                                      aCGTGCATCACCGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCGCCTTCCGCCAGGAGCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGcc                                                                                                                                          <  1:926903/173‑1 (MQ=255)                                                                                                                                                            aTCAACGGACCGCCCTTCCCACCAGGAAAAACGGCAGGTTTCAGTTTTTTGTCCAGCTTTTCCCTACCACCTTTCGCCAGGATCAGGCCGCCGCGGGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACg                                                                  <  1:1461851/239‑1 (MQ=255)                                                                                                                                                              cAACGGACCGCCCTGACCCCCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGa                                                                                                                                                                                                              <  1:520255/97‑1 (MQ=255)                                                                                                                                                              cAACGGACCGCCCTGACCCCCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTtgt                                                                                                                                                                                              <  1:1234847/113‑1 (MQ=255)                                                                                                                                                                      cccccTGACCACCAGGAAAAACGGCAGAGTGCAGTTTTTTGTACAGCTCTTCGCTCCCACCTTTCGCCAGGATCGGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGcc                                                        <  1:1345603/237‑1 (MQ=255)                                                                                                                                                                      cGCCCTGACCCCCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGcc                                                        <  1:1154966/239‑1 (MQ=255)                                                                                                                                                                              ccaccaGGAAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTccac                                                <  1:1195143/239‑1 (MQ=255)                                                                                                                                                                               ccccggggAAAACGGGAGAGTTCAGTTTTTTTTACAGCTATGCGTTACCACCTTTCGCCGGGATCGGGCCGCCTCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTccacc                                               <  1:1753693/234‑1 (MQ=255)                                                                                                                                                                                    gggAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACg                                                                  <  1:1559820/215‑1 (MQ=255)                                                                                                                                                                                    gggAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTg                                                                                                                  <  1:1127896/167‑1 (MQ=255)                                                                                                                                                                                      gAAAACGGCAGAGTTCAGTTTTTTGTACAGGTTTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACg                                            <  1:1623032/235‑1 (MQ=255)                                                                                                                                                                                        aaaCGGCAGATTTCAGTTTTTTGTACATCTCTTCGCTACCCCCTTTCCCCAGGATCAGGCCGCCGCGCGGCCCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGa                                       <  1:1293018/238‑1 (MQ=255)                                                                                                                                                                                         aaCGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGTTCCCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTa                                                                                             <  1:1428828/183‑1 (MQ=255)                                                                                                                                                                                              cAGAGTTCAGTTTTTTGTACAGCTCTTCGCTTCCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAgtg                                                                                                                                                                                         <  1:1505679/86‑1 (MQ=255)                                                                                                                                                                                                  gTTCATTTTTTTGTACACCTCTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCATCACGCCGGAATATGCagag                             <  1:544161/238‑1 (MQ=255)                                                                                                                                                                                                             tGTACAGCTCTTCGCCACCACCt                                                                                                                                                                                                                                         >  2:656125/1‑23 (MQ=255)                                                                                                                                                                                                                aCAGCTCTTCGCTACCACCTTTc                                                                                                                                                                                                                                      >  2:445602/1‑23 (MQ=255)                                                                                                                                                                                                                  aGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCa                                                                                                                               <  1:946800/124‑1 (MQ=255)                                                                                                                                                                                                                  aGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCa                                                                                                                               >  2:946800/1‑124 (MQ=255)                                                                                                                                                                                                                    ctctTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGCCGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATCTGCAGAGAAACCACCGATAATCAttt          >  1:802154/1‑239 (MQ=255)                                                                                                                                                                                                                    ctctTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCAttt          >  1:1958312/1‑239 (MQ=255)                                                                                                                                                                                                                    ctctTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGAAGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCCCGCCGGAATATGCAGAGAAACCACCGATAATCAttt          >  1:1755612/1‑239 (MQ=255)                                                                                                                                                                                                                    ctctTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGGTTTGGGAGTGGTGGTAGTaa                                                                                                                                                                               >  2:1150490/1‑74 (MQ=255)                                                                                                                                                                                                                         cGCTCCCGCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGc     <  1:1688658/239‑1 (MQ=255)                                                                                                                                                                                                                            tACCCCCTTTCGCCAGGATCGGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTtg  <  1:1835764/239‑1 (MQ=255)                                                                                                                                                                                                                              ccaccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCg                                         >  2:1828232/1‑198 (MQ=255)                                                                                                                                                                                                                              ccaccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCg                                         <  1:1828232/198‑1 (MQ=255)                                                                                                                                                                                                                              ccaccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAata                                    <  1:1172958/203‑1 (MQ=255)                                                                                                                                                                                                                              ccaccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAata                                    >  2:1172958/1‑203 (MQ=255)                                                                                                                                                                                                                              ccaccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGAAGCGAAACAACCGATAGTCATTTTCgg      >  2:730956/1‑233 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                           GTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAG‑CGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTG  >  NC_000913/2684557‑2685013

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑