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

TCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTG  >  NC_000913/2684543‑2685013                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           tcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTaccacc                                                                                                                                                                                                                                          <  1:31294/239‑1 (MQ=255)          aaCAGGTGGTTATCAGTGCCGCCGGAAACCCCTTTGTAGCCGCGCTCGAGGAAAACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAAAAGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGc                                                                                                                                                                                                                                    <  1:1807362/236‑1 (MQ=255)                    tATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGgccg                                                                                                                                                                                                                       <  1:813346/239‑1 (MQ=255)                                    aaCCACTTTGTACCCGCGCTCGAGGACCCTTTCTACCATCGCTTTAGCGTTTTTAGCGACCTTCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCACCCGCTTTACCGGCGATTACGTGCATCAACGGCCCGCCCTGACCCCCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGcc                                                                                                                                                                                                        <  2:194400/238‑1 (MQ=255)                                                 ccGCGCCGGGGGAACACTTCTCCCATCGCTTTGGCGTTTTTGGCGCCCTGCTGCTGTTAATTTTTGACCTCAGGCACCATCGCTTCTTCCAGAGCACCCGCTTTCCCGGCGATTCCGTGCATCAACGGACCGCCCTGACCCCCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAg                                                                                                                                                                                           <  2:791254/238‑1 (MQ=255)                                                  cgcgCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACGCGGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGCCCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAgtg                                                                                                                                                                                         <  2:1473229/239‑1 (MQ=255)                                                  cccgcTGGAGGAACCCTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGACCCCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACGGCTCTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAgtg                                                                                                                                                                                         <  2:121792/237‑1 (MQ=255)                                                           ggAACTGTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGTTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGgtagt                                                                                                                                                                                 <  1:1542561/238‑1 (MQ=255)                                                                tcTTCCACCATCTCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCGGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGCCCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTCCCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTaacaac                                                                                                                                                                           <  2:1014732/238‑1 (MQ=255)                                                                         aTCGCTTTAGCGTTTTTAGCGACCGGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATg                                                                                                                                                                  <  2:347490/239‑1 (MQ=255)                                                                                    gTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCtt                                                                                                                                                                                                                                        <  2:1750780/158‑1 (MQ=255)                                                                                              aCCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCa                                                                                                                                                                                                       <  1:859863/181‑1 (MQ=255)                                                                                                ctgctgcggGTAATTTTTGAACTCAGGACCCATCGCTTCTTTCAGAGCACCCGCTTTCCGGACGTTTACGTGCATCAACGGCCCGCCCTGCCCCCCAGGGAAAACGGCAGAGTGCATTTTTTTGTCCAGCTCTTCGCTCCCCCCTTCCGCCAGGATCAGGCCGCCGCGGGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGc                                                                                                                                           <  2:1814725/239‑1 (MQ=255)                                                                                                             gTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCCCCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCCCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGCTTTTTGAGTGGTGGTAGTaacaac                                                                                                                                                                           >  1:1486816/1‑194 (MQ=255)                                                                                                              ttttGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTCCCCCCTTTCGCCAGGAGCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAgtg                                                                                                                                                                                         <  1:1452056/179‑1 (MQ=255)                                                                                                              ttttGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAgtg                                                                                                                                                                                         <  1:681459/179‑1 (MQ=255)                                                                                                              ttttGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATg                                                                                                                                                                  <  1:1708103/202‑1 (MQ=255)                                                                                                                  gAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCg                                                                                                                         <  1:853097/239‑1 (MQ=255)                                                                                                                                        ttCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGaa                                                                                                   <  1:1447622/239‑1 (MQ=255)                                                                                                                                                   cGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGc                                                                                                                                                                                                                                    <  1:365690/99‑1 (MQ=255)                                                                                                                                                    gCTTTCCCGGGGATTACATGCATCAACGGACCGCCCTGCCCCCCAGGGAAACGGGCAGAGTTCGGTTTTTTGTACAGCTCTTCGCTACCCCCTTTCGCCAGAATCAGGCCGCCGCGCGACCCCGCCAGGGTTTtgtg                                                                                                                                                                                             <  2:454212/137‑1 (MQ=255)                                                                                                                                                        tACCGGCGATTCCGTGCATCAACGGACCGCCCTGACCACCAGGAAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATg                                                                                   <  1:263713/239‑1 (MQ=255)                                                                                                                                                                           ccttatccggCCCTGGCCACCAGGAAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTCCCCCCTTTCGCCAGGATCAGGCCGCCGGGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCa                                                                <  1:1100751/230‑1 (MQ=255)                                                                                                                                                                           cAAGGGACCGCCCTGACCACCAGGAAAAGCGGCAGAGTTCATTTTTTTGTACAGCTCTTCGCTACCCCCTTTCCCCAGGATCGGGCGGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAAc                                                                                                                                                            <  1:1636992/147‑1 (MQ=255)                                                                                                                                                                            aaCGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACTTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGa                                                                                                                                                                                            <  1:297718/114‑1 (MQ=255)                                                                                                                                                                               ggACCGCCCTGACCCCCAGGGAAAACGGCAGAGTTCATTTTTTTGTACAGCTCTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCAtttt                                                            <  1:1092922/239‑1 (MQ=255)                                                                                                                                                                                gACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGtggtgg                                                                                                                                                                                     <  1:650593/117‑1 (MQ=255)                                                                                                                                                                                   cGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGcc                                                        <  1:679857/239‑1 (MQ=255)                                                                                                                                                                                     cccTGACCACCAGGGAAAACGGCAGTGTTCAGTTTTTTGTCCAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGTCGCCACGCGGACCCGCCAGGGTTTTGTGAGTGGAGGCAGTAACAACGTAAGCATAGGGAACCGGGTTCcgg                                                                                                                                                  >  2:1636992/1‑147 (MQ=255)                                                                                                                                                                                     cccTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGTTCCCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAgtg                                                                                                                                                                                         <  1:21544/108‑1 (MQ=255)                                                                                                                                                                                     cccTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGaa                                                                                                   <  1:675965/194‑1 (MQ=255)                                                                                                                                                                                           cccccaGGGAAAGCGGCAGAGTTCGGTTTTTTGTGCAGCTCTTCGCTACCCCCTTTCGCCAGGAGCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTccac                                                <  1:846126/236‑1 (MQ=255)                                                                                                                                                                                           cccccaGGGAAAACGGCAGAGTTCAGTTTTTTGTACACCCCTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGNa                                                                                                                                                                                                              <  2:921326/78‑1 (MQ=255)                                                                                                                                                                                            caccaGGAAAAACGGCAGAGTTCATTTTTTTGTACAGCTCTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAg                                                                                                                                                                                           <  2:193561/99‑1 (MQ=255)                                                                                                                                                                                                   gAAAACGGCGGGGTTCAGTTTTTTGTACAGCTCTTCGATACCGCCTTTCGCCAGGGTCGGGCGGCCGCGGGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAGCGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATg                                                                                   <  1:826787/196‑1 (MQ=255)                                                                                                                                                                                                                                 ctctTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCg                                                                                                                                                                                                          <  1:1247412/47‑1 (MQ=255)                                                                                                                                                                                                                                 ctctTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCg                                                                                                                                                                                                          >  2:1247412/1‑47 (MQ=255)                                                                                                                                                                                                                                 ctctTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCAttt          >  1:272804/1‑239 (MQ=255)                                                                                                                                                                                                                                 ctctTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCAttt          >  2:901306/1‑239 (MQ=255)                                                                                                                                                                                                                                 ctctTCGAGACCACCTTTGGCCGGGATCGGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAg                                                                                                                                               <  2:1730771/106‑1 (MQ=255)                                                                                                                                                                                                                                  tctTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCAttt          <  1:343809/238‑1 (MQ=255)                                                                                                                                                                                                                                    ttCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTtgt                                                                                                                                                                                              <  1:248218/56‑1 (MQ=255)                                                                                                                                                                                                                                    ttCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTtgt                                                                                                                                                                                              >  2:248218/1‑56 (MQ=255)                                                                                                                                                                                                                                      cGCTACCACCTTTCGCCGGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAgtg                                                                                                                                                                                         <  2:1253743/59‑1 (MQ=255)                                                                                                                                                                                                                                      cGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAgtg                                                                                                                                                                                         >  1:1253743/1‑59 (MQ=255)                                                                                                                                                                                                                                      cGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGt                                                                                              >  1:701662/1‑150 (MQ=255)                                                                                                                                                                                                                                      cGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGt                                                                                              <  2:701662/150‑1 (MQ=255)                                                                                                                                                                                                                                      cGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCg       >  1:1448144/1‑237 (MQ=255)                                                                                                                                                                                                                                      cGCTACCACCTCTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCg       >  1:1748271/1‑237 (MQ=255)                                                                                                                                                                                                                                         tACCCCCTTTCGCCAGGATCAGACCGCGGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTtg  <  1:1244374/239‑1 (MQ=255)                                                                                                                                                                                                                                         tACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGa                                                                                                                                                                                            <  2:451828/53‑1 (MQ=255)                                                                                                                                                                                                                                           ccaccTTTCGCTAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGTGCATGAGGAACCGGGTTCGGGTAGACGCcagcag                                                                                                                                     >  1:1730771/1‑106 (MQ=255)                                                                                                                                                                                                                                           ccaccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGaa                                                                                                   >  2:795750/1‑140 (MQ=255)                                                                                                                                                                                                                                           ccaccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGaa                                                                                                   <  1:795750/140‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           TCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTG  >  NC_000913/2684543‑2685013

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