Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,684,779	0	A	C	36.5%	81.0 / 18.3	76	G243G (GGT→GGG)	glyA	serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (11/36);  new base C (2/25);  total (13/63)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCT  >  NC_000913/2684552‑2685011                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         agcAGGTGGTTATCGGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCGTCCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGAAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGcc                                                                                                                                                                                                                                 <  2:1911352/235‑1 (MQ=255)      gtggtTATCAGTGCCGCCGGAAACAACTTTGTAGCCGCGCTAGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCCTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGA‑CAGCCCTGACCACCAGGGAAAAAGGCAGAGCTCAGGTTTTTGTACAGCTCTT‑CGCTACCCACTTTCGCCAGGATCa                                                                                                                                                                                                                          >  2:1636566/1‑239 (MQ=255)                           aaCCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCa                                                                                                                                                                                                     <  1:416375/239‑1 (MQ=255)                                 tttGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CTCTACCACCt                                                                                                                                                                                                                                       <  2:805421/199‑1 (MQ=255)                                    gTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGATTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTtgt                                                                                                                                                                                            <  1:1658407/239‑1 (MQ=255)                                    gTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTtgt                                                                                                                                                                                            <  1:1984514/239‑1 (MQ=255)                                         cgcgcgcgGGGAACCATTCTCCCTTCGTTTTAGCTTTTTTACGGCCCGGCTGCTGTTAAGTTTTGAACCCAGGCTCCATCGCTTTTTTAAAACCACCCGCTTACCCGGCGATTAGGTGAATCAACGGA‑CCGCCCTGACCACCGGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAgtg                                                                                                                                                                                       <  2:33847/239‑1 (MQ=255)                                         cgcgCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAGGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAGCCGCTTTACCGGCGATTACGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAgtg                                                                                                                                                                                       <  1:1178966/239‑1 (MQ=255)                                         cgcgCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAg                                                                                                                                                                                         <  1:434823/237‑1 (MQ=255)                                               cGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGCCCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGgta                                                                                                                                                                                 <  1:1754269/239‑1 (MQ=255)                                                             cccATCGTTTAGGCGTTTTTAGCGCCCGGCTGCTGGTAAGTTTTGAACTCAGGCCCCATCGCTTCTTGCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGC‑CCGCCCGGACCCCCGGGGAAACCGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGc                                                                                                                                                                   <  2:709223/238‑1 (MQ=255)                                                                  cGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGACCTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGgtagt                                                                                                                                                                               <  1:438458/222‑1 (MQ=255)                                                                      ttAGCGTTTGTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCACCCGCTTTACCGGCGATTCCGGGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAAc                                                                                                                                                          <  2:96422/239‑1 (MQ=255)                                                                             tttgAGCGCCCTGTGGCCGGTAATTTTGAAACTCAGGCTCCATCGCTTCTTCCAGGGCAACCGCTTTACCGGCGATTACGGGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCATTTTTTTGTACAGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTc                                                                                                                                                   <  2:1948249/239‑1 (MQ=255)                                                                                aaGCGCCCTGCTGCTGGTAATTTTTGACCTCAGGCTCCACCGCTTCTTACAGAGCAACCGCTTTCCCGGCGATTCCGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTCCCCCCTTCCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCggg                                                                                                                                                <  2:253767/238‑1 (MQ=255)                                                                                    gACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATGGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGA‑CCGCTCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTa                                                                                                                                              <  1:960822/237‑1 (MQ=255)                                                                                    gACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTa                                                                                                                                              <  1:2244622/237‑1 (MQ=255)                                                                                        tgctgctgGTAAGTTTGGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGcc                                                                                                                                        <  1:1436644/239‑1 (MQ=255)                                                                                           tgctgGGAAGTGTGGAACTCAGGCTCAATCGCTGTTTTCAGGGCAACCGCTTTCCCGGCGATTACGTGCATCAAGGGA‑CCGCCCGGACCACCAGGGAAAACGGCGGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCcag                                                                                                                                      <  2:638962/238‑1 (MQ=255)                                                                                                  aatttttGACCTCAGGCTCCACCGCTTCTTTCAGAGCAACCGTTTTACCGGCGATTACGTGCATCAGCGGC‑CCGCCCTGCCCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAAc                                                                                                                               <  2:525544/235‑1 (MQ=255)                                                                                                     ttttGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCACCGGA‑CCGCCCTGACCCCCAGGGAATACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGtggtgg                                                                                                                                                                                   <  1:1785033/183‑1 (MQ=255)                                                                                                     ttttGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGC‑CCGCCCTGCCCCCCGGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTCCCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGtttt                                                                                                                                                                                              <  1:734431/172‑1 (MQ=255)                                                                                                            cTCAGGCTCCTTCGCTTCTTTCAGAGCACCCGCTTTACCGGCGATTACGTCCATCAACGGC‑CCGCCCTGACCCCCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAAcc                                                                                                                                                         <  2:681526/202‑1 (MQ=255)                                                                                                             tCAGGCTCCAGCGCTTCTTCCAGAGCAACCGCTTTACCGGCGATTACGGGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAAc                                                                                                                   <  1:416361/239‑1 (MQ=255)                                                                                                                  cTCCATCGCTTCTTTCAGAGCAACCTCTTTCCCGGCTATTACGTCCATCACCTGA‑CCGCCCTGACCACCAGGGAAATCGGCAGAGTTCATTTTTTTGTACAGCTCTT‑CGCTCCCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACg                                                                                                                  <  2:1376787/235‑1 (MQ=255)                                                                                                                  cTCCATCGCTTCTTCCAGAGCACCAGCTTTACCGCCGATTACGTGCATCAACGGA‑CCGCCCTGACCACAAGGAAAACCGGCAGAGTTCATTTTTTGGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTgc                                                                                                               <  2:392546/238‑1 (MQ=255)                                                                                                                   tCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTgcgc                                                                                                             <  1:1970284/239‑1 (MQ=255)                                                                                                                   ccactCGCTTCTTTCAGACCAACCCTTTTACCGGCAATTACGGGCATCAACGGC‑CCCCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACCGCTCTT‑CGCTCCCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTgcg                                                                                                              <  2:1705945/234‑1 (MQ=255)                                                                                                                       tCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGAATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACCGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGACAGGGTTTTGTGAGTGGTGGTAGTAACCACGCGAGCATGAGGAACCGGGTTCGGGTAGCCGCCAGCCGCCACCAGTCCCGCAACGTGCGCCa                                                                                                           >  2:1272955/1‑237 (MQ=255)                                                                                                                         gCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAg                                                                                                                                                                                         <  2:1694934/157‑1 (MQ=255)                                                                                                                          tttcttcCACAGCAACCGTTTTACCGGCGATTACGT‑CCTACACGGATCCGCCCTGCCCACCAGGGAAACCGGCAGAGTTCAGTTTTTTGTACAGCTCTG‑CGCTCCCTCCTTCCGCCAGTATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCa                                                                                                      <  2:1742196/238‑1 (MQ=255)                                                                                                                                       aaCCGCTTTACCGGCGATTCCGTGCATCAACGGA‑CCGCCCTGACCACCAGGAAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAg                                                                                         <  1:1841107/239‑1 (MQ=255)                                                                                                                                            cTTTACCGGCGATTACGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAg                                                                                                                                                                                         <  1:1663383/138‑1 (MQ=255)                                                                                                                                             tttACCGGCGATTCCGTGCATCAACGGA‑CCGCCCTGACCCCCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCa                                                                                       <  1:1444766/235‑1 (MQ=255)                                                                                                                                             tttACCGCCGATTCCGGGCATCAACGGA‑CCGCCTGGACCACAAGGAAAACCGGAAGATTTCAGTTTTTTTTACAGCTTTT‑CGCTACCCCCTTCCGCCGGGATCAGGCCGCCGCGGGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCa                                                                                       <  2:2217071/235‑1 (MQ=255)                                                                                                                                                     cGATTACGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCt                                                                               <  1:2195316/235‑1 (MQ=255)                                                                                                                                                           cGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGt                                                                                                                                                                                                 <  1:780604/115‑1 (MQ=255)                                                                                                                                                             tGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCcgccgc                                                                                                                                                                                                                 <  1:147445/97‑1 (MQ=255)                                                                                                                                                                aTCAACGGC‑CCGCCCTGACCCCCAGGCAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACg                                                                <  2:2027133/239‑1 (MQ=255)                                                                                                                                                                  cAACGGC‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCGCCGa                                                                                   <  1:57551/218‑1 (MQ=255)                                                                                                                                                                   aaCGGA‑CCGCCCTGACCACCGGGGAAAACGGAAGAGTGCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCGGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCGCCGATGCTGTCAGCGATTTCACGCAt                                                             <  1:1500024/239‑1 (MQ=255)                                                                                                                                                                        a‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTCCCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTc                                                         <  1:134048/238‑1 (MQ=255)                                                                                                                                                                           cGCCCTGACCCCCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTTTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGcc                                                      <  1:1801284/239‑1 (MQ=255)                                                                                                                                                                              ccTGCCCACCAGGGAAAACGGCAGAGTGCAGTTTTTTGTACAGCTCTC‑CGCTACCTCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCa                                                    <  1:1451010/238‑1 (MQ=255)                                                                                                                                                                                 gACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCCCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCACCGATTTCACGCATTTTCCCCCAGTcc                                                >  2:16286/1‑239 (MQ=255)                                                                                                                                                                                  accaccAGTGAAATCGGCAAATTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTcca                                               <  1:1636566/239‑1 (MQ=255)                                                                                                                                                                                   ccaccaGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTccac                                              <  1:549372/239‑1 (MQ=255)                                                                                                                                                                                    caccaGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGcc                                                      <  1:1482733/230‑1 (MQ=255)                                                                                                                                                                                           gAAACCGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACg                                          <  1:1279479/235‑1 (MQ=255)                                                                                                                                                                                             aaaCGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTCCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGAATTCACGCATTTTCGCCCAGTCCACCACGCCGGa                                     <  1:1481120/238‑1 (MQ=255)                                                                                                                                                                                                       gTTCATTTTTTTGTACCGCTCTT‑CGCTACCCCCTTGCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGc                               <  1:1409646/234‑1 (MQ=255)                                                                                                                                                                                                       gTTCAGTTTTTTGTACACCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCa                              <  1:1910438/235‑1 (MQ=255)                                                                                                                                                                                                        ttCAG‑TTTTTGTACCACCCTTCCCCTCCCCCCTTTCGCCGGGTCCAGGCCCCCGCGCGGCCCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAgaga                          <  2:620840/238‑1 (MQ=255)                                                                                                                                                                                                              tttttGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAacca                     <  1:2040683/237‑1 (MQ=255)                                                                                                                                                                                                                ttttTACGGCTGTT‑CGCAACCACTTGCCGCCAGGGTCAGCCCGCGGCGCGAGACCGCCAGGGTTTTGAGAGTGGTGGTAGTAACAACGTGAGCAGGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGa                 <  2:1465967/239‑1 (MQ=255)                                                                                                                                                                                                                tttGTACAGCTCTT‑CGCTaccac                                                                                                                                                                                                                                         >  2:1983388/1‑23 (MQ=255)                                                                                                                                                                                                                  tGTACAGCTCTT‑CGCTACCACCt                                                                                                                                                                                                                                       >  2:150215/1‑23 (MQ=255)                                                                                                                                                                                                                    tACAGCTCTT‑CGCTACCCCCTTTCGCCAGGATCAGCCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAAt             <  1:1983388/239‑1 (MQ=255)                                                                                                                                                                                                                     aCAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAAt             <  1:649972/238‑1 (MQ=255)                                                                                                                                                                                                                       aGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCa                                                                                                      >  2:1139948/1‑147 (MQ=255)                                                                                                                                                                                                                       aGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCa                                                                                                      <  1:1139948/147‑1 (MQ=255)                                                                                                                                                                                                                        gCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGtggtgg                                                                                                                                                                                   >  2:1176398/1‑69 (MQ=255)                                                                                                                                                                                                                        gCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGtggtgg                                                                                                                                                                                   <  1:1176398/69‑1 (MQ=255)                                                                                                                                                                                                                         ctttt‑CGCAACCACTTTTGGCCAGGATCGGGCCGCCGGGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCa                                                                                                           <  2:1619231/137‑1 (MQ=255)                                                                                                                                                                                                                         ctctT‑CGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAg                                                                                         <  2:1070470/158‑1 (MQ=255)                                                                                                                                                                                                                         ctctT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTaacaac                                                                                                                                                                         <  2:2020703/78‑1 (MQ=255)                                                                                                                                                                                                                         ctctT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTaacaac                                                                                                                                                                         >  1:2020703/1‑78 (MQ=255)                                                                                                                                                                                                                         ctctT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAg                                                                                         >  1:1070470/1‑158 (MQ=255)                                                                                                                                                                                                                         ctctT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCAttt        >  1:1899827/1‑239 (MQ=255)                                                                                                                                                                                                                         ctctT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCAttt        >  2:367294/1‑239 (MQ=255)                                                                                                                                                                                                                            tt‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCa                                                                                                           >  2:1011233/1‑137 (MQ=255)                                                                                                                                                                                                                            tt‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCa                                                                                                           <  1:1011233/137‑1 (MQ=255)                                                                                                                                                                                                                            tt‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAac                       >  2:2182989/1‑221 (MQ=255)                                                                                                                                                                                                                            tt‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAac                       <  1:2182989/221‑1 (MQ=255)                                                                                                                                                                                                                               cGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGc                                                                         <  2:445938/169‑1 (MQ=255)                                                                                                                                                                                                                               cGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGc                                                                         >  1:445938/1‑169 (MQ=255)                                                                                                                                                                                                                                 cTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGa                                                                                                                                                                                                            <  2:941068/36‑1 (MQ=255)                                                                                                                                                                                                                                    ccaccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAAc                                                                                                                               >  1:1680064/1‑110 (MQ=255)                                                                                                                                                                                                                                    ccaccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAAc                                                                                                                               <  2:1680064/110‑1 (MQ=255)                                                                                                                                                                                                                                    ccaccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACg                                                                                                                  >  2:2166244/1‑123 (MQ=255)                                                                                                                                                                                                                                    ccaccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACg                                                                                                                  <  1:2166244/123‑1 (MQ=255)                                                                                                                                                                                                                                    ccaccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGaa                                                                                                 >  1:1619231/1‑140 (MQ=255)                                                                                                                                                                                                                                    ccaccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCt  >  1:36171/1‑235 (MQ=255)                                                                                                                                                                                                                                    ccaccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCt  >  1:806265/1‑235 (MQ=255)                                                                                                                                                                                                                                     caccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGc                                                                         >  1:2064512/1‑163 (MQ=255)                                                                                                                                                                                                                                     caccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGc                                                                         <  2:2064512/163‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         AACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGA‑CCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTT‑CGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCT  >  NC_000913/2684552‑2685011

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