Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,185,238	0	A	G	24.6%	140.3 / 13.8	69	S22P (TCC→CCC)	thiM	hydroxyethylthiazole kinase
Reads supporting (aligned to +/- strand):  ref base A (18/34);  new base G (17/0);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.15e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTA  >  NC_000913/2184993‑2185484                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         caccGCTACTGGATCAAGGTTCCAGGGTGTTTGAGAGCTTTTGGCTTGCGCAACGGCAGCCCGCATCGTCTGAGCGCGTGGCGGCCTCAGTTTGCCAACGTTAATCAACAAGGAACTGGCGATAGCCGCAAACTGACGGGCCTCTTCGTTTTCGATAACCAGCGCTGGCGATGCACCGAGCGCCAGCAGGGTATTGGCGGGAAA‑GAT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGtt                                                                                                                                                                                                                                                     <  2:386460/251‑1 (MQ=255)     gCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGG‑TTGCACCACATCCTTGGGCA‑TGCAGTGCGCAAG‑AGGGGGATGTTgggg                                                                                                                                                                                                                                                 >  1:60727/1‑249 (MQ=255)        aCTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑aaa                                                                                                                                                                                                                                             <  1:326288/251‑1 (MQ=255)              tCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTGAA‑GGG‑TTGGACCACAGCCTTGGGCAGTG‑AGTGCACAAG‑AGGGGGAGGTTGGGG‑GAAAAgggg                                                                                                                                                                                                                                       >  1:344732/1‑249 (MQ=255)               cAAGCGTCCAGGGTGGTTGAAAGCTTTGTGCTTGCTCAACGGCAGCCCGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGAGAGCCGCAAACTGACTGGCCTCTGCGGTTTCGATAACCAGCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGgtgt                                                                                                                                                                                                                                      <  1:138853/251‑1 (MQ=255)                                 gagaGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑AGT‑TTGGCCCCCATCCTTGGGCA‑TGCAGTGGACAAGGAGGGGA‑TGTTGGGG‑GAAAAGgtggctcggtttgggtttcgtt                                                                                                                                                                                                                    >  1:89786/1‑232 (MQ=255)                                            gcttgctCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCca                                                                                                                                                                                                         <  2:84163/250‑1 (MQ=255)                                               tgctCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCcagc                                                                                                                                                                                                       <  1:79463/249‑1 (MQ=255)                                                  ccAACGGCCGCACGCATCGCCTGAGCGGGTGGGTGCGGCAGTGTTCCCACGTGAAGCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTTGCGGGAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAT‑AGGTGAATGTTGGTG‑AAAAAGGTGTAACGCGTACGCAGATTGCTCT‑GAACCCAGCAGGTc                                                                                                                                                                                                  <  2:20104/250‑1 (MQ=255)                                                   cAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTc                                                                                                                                                                                                  <  1:473489/250‑1 (MQ=255)                                                    aaCGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGG‑TTTCACCCCAACATTGGTCA‑TGCAGTGCCCAAG‑AGGGGGATGGTGGTG‑AAAAAGGTGTAACCCGCGCGCAGATTTCGCC‑GAAACCCGCAgggcgg                                                                                                                                                                                                >  1:69813/1‑247 (MQ=255)                                                          agcaCGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGG‑TTGCACCACATCATTGGTCA‑TGCAGGGCACAAG‑AGGGGGATGTTGGGG‑AAAAAGGGGTAACGCGTGCGCAGATTGGGCT‑GAACCCAGCAGGGCGGCTTGc                                                                                                                                                                                           >  1:192385/1‑250 (MQ=255)                                                            aaCGCATCGCCTGGGCGCGGGGCGGCGTCAGTGTGCAAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTTTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATa                                                                                                                                                                                        <  2:56738/250‑1 (MQ=255)                                                             aCGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGGAAA‑GGTGTGGCACCCCTCCATTG‑TCA‑TTCCGGGGACAAG‑AGGGGAATTGTGGGG‑GAAAAGGTGGAACGCGTGCGCAGATTGTGCT‑GTACCCAGAGGGTCGagtttcaaat                                                                                                                                                                                       >  2:374608/1‑242 (MQ=255)                                                               gCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGGAAA‑GGG‑TTGCACCACACTATTTGTCA‑TTCCCGGCGCGAG‑AGGGGGAAGTTGGGGGAACAAGGGGGAA‑GCGTGCGCGAATTGCGCT‑GAACCCAGCAGGGCGGCTTGCATAGtt                                                                                                                                                                                     >  1:13900/1‑251 (MQ=255)                                                                        gAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCt                                                                                                                                                                              <  2:434895/249‑1 (MQ=255)                                                                                      gTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTgaagaa                                                                                                                                                              <  1:422084/251‑1 (MQ=255)                                                                                      gTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGACAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGG‑TTTCACCACATCATTGGTCA‑TGCAGGGCACAAG‑AGGGGGATGGTGGGG‑GAAAAAGTGGAAAGCGTGCGCCGGTTGCGCC‑GAACCCCGCCGGGCGGCTTGCATTGTTTTCTCCTTCCATaaacggaaga                                                                                                                                                               >  1:487274/1‑242 (MQ=255)                                                                                        cAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGG‑TTGCACCACATCATTGGTCA‑TGCCGGGCACAAG‑AGGGGAATTTTGGGG‑AAAAAAGTGGAAAGCGTGCGCAGATTGCGCT‑GAACCCAGCCGGTCGGCTTGCCAAGGTTGCTCCTTCCCTAACGGGAAAAAGc                                                                                                                                                            >  1:68963/1‑251 (MQ=255)                                                                                        cAGTGTGCCAACGTAAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGc                                                                                                                                                            <  2:363300/251‑1 (MQ=255)                                                                                         aGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACCGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGG‑TTGCACCACATCATTGGGCA‑CGCCGAGCACCAC‑GGGGGGATGTTGGGG‑GAAAAGGTGTTAAGCGGGCGCAGATTGCGCTTGAACC‑ACCAGGGCGGCTTGCATAGGTTGCTCCTTCCACAACggggcgaagg                                                                                                                                                            >  1:278271/1‑241 (MQ=255)                                                                                             tgCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGG‑TTGCACCACATCCTTGGTGA‑TGCCGTGCACAAG‑AGGGGAATGGTGGTG‑GAAAAGGTGTAACGCGGGCGCAGATTGCGCT‑GAACCCCGCAGGGCGACCTTCATAGTTTGCTCCTGCCATAAAGTGAAGAAGCAATg                                                                                                                                                        >  1:295206/1‑250 (MQ=255)                                                                                                      tAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGG‑TTGGAACACATCCTTGGTCA‑TTCAGTGCACAAG‑AGGGGGAAGTTGGTG‑AAAAAAGTGTAACGCGTGCGCCGAATGGGCT‑GAACCCAGCAGGGCCGCTTGCATAGTTTGCTCCTGCCATAActggcagaaggcacgacccggtggg                                                                                                                                              >  1:94714/1‑227 (MQ=255)                                                                                                       aaTCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACCTCATTGGTCA‑TGCAGGGCCCAAG‑AGGGGGATGTTGGGG‑AAAAAGGGGTAACGCGGGCGCAGAtt                                                                                                                                                                                                                      >  1:402379/1‑179 (MQ=255)                                                                                                       aaTCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGAtt                                                                                                                                                                                                                      <  2:402379/179‑1 (MQ=255)                                                                                                           aacaaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGGAAA‑GGG‑TTGCGCCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAAAGTTGGTG‑GAAAAGGGGTAAAGCGTGCGCAGATTGCGCT‑GGAACCAGCAGGGCGACTTGCATAGTTTGCTCCCGCCCTAAAGTGGAGAAGCCATGACCTGGTGGGCCGTg                                                                                                                                         >  1:29618/1‑251 (MQ=255)                                                                                                             caaGGCCCTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCG‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGAc                                                                                                                                       <  2:192385/251‑1 (MQ=255)                                                                                                              agggCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACt                                                                                                                                      <  1:196661/249‑1 (MQ=255)                                                                                                              aaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCAACACAACATTGGTCA‑CGCCGTGCACAAG‑AGGGGGGGGTTGGGG‑AAAAAGGGGTAACGCGTGCGGAAATTGCGGT‑GAACCCAGGAGGGGGACTTGGATAGTTTGCTCCTGCCAGAAAGGGAAGAAGCAATGAACTGGGGGGCCGTggtt                                                                                                                                      >  1:467588/1‑248 (MQ=255)                                                                                                                   aCTGGCGATAGCCGCAAACTGACCGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCt                                                                                                                                 <  1:428170/251‑1 (MQ=255)                                                                                                                      ggCGATAGCCGCAAACTGGCTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGTGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATGGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GACCCCAGCAGGTCGACTTGCATAGTT‑GCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTAc                                                                                                                               <  1:29852/249‑1 (MQ=255)                                                                                                                      ggCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACg                                                                                                                              <  2:452878/251‑1 (MQ=255)                                                                                                                        cgagaGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGGCCGTGGCTTCCCTACGCt                                                                                                                            <  2:415631/247‑1 (MQ=255)                                                                                                                        cGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCt                                                                                                                            <  2:194976/251‑1 (MQ=255)                                                                                                                         gATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTgcgc                                                                                                                                                                                                                  <  2:31476/165‑1 (MQ=255)                                                                                                                         gATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGG‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAAGTGTAACGCGTGCGCAGATTgcgc                                                                                                                                                                                                                  >  1:31476/1‑165 (MQ=255)                                                                                                                               cGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGG‑TTGCACCACATCATTGGGTA‑TGCAGTGCACAAG‑AGGGGGATGTTGGGG‑AAAAAAGGGTAACGCGTGCGCAGATTTGGCC‑GAAACCAGCAGGGCGACTTGCATAGTTTGCTCCCTCCCTAACGTGAAGAAGCAATGACCCGGTGGGCCGTGACttaccccccccgggctt                                                                                                                      >  2:319558/1‑235 (MQ=255)                                                                                                                                    actgactgGCCGTTTCGGTGTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAg                                                                                                                <  2:192588/251‑1 (MQ=255)                                                                                                                                     ctgactgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGTGTTG‑ACCACCTCATTTGGCA‑TGCAGGGCACAAG‑AGGGGGGTGTTGGGG‑AAAAAGGGGTAACGCGTGCGCAGATTTCGCT‑GAACCCACCAGGTCGGCTTGCCTAGTTTTCTCCTGCCATAAAGTGCAGAAGGAATGACCAGGTGGGCCGTGACTTCCCCACGCTGGCATTAttccg                                                                                                                >  2:445528/1‑246 (MQ=255)                                                                                                                                          tgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGG‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGGATGTTGGGG‑AAAAAGGTGTAACGCGGGCGCAGATTGCGCT‑GAACCCAGCAGGTCGGCTTTCATAGGTTGCTCCTGCCCTAACGTGGAGAAGGAATGACCTGGGGGGCCGTGACTTCCCCACGCCGGCATTATCCAGATCAgg                                                                                                          >  1:130565/1‑251 (MQ=255)                                                                                                                                           ggCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAg                                                                                                           <  2:411806/249‑1 (MQ=255)                                                                                                                                                  tCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATAc                                                                                                    <  1:243876/249‑1 (MQ=255)                                                                                                                                                  tCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACg                                                                                                   <  1:480421/250‑1 (MQ=255)                                                                                                                                                     gTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCAACCCATCATTGGTGA‑TGCAGTGGACAAG‑GGGGGGATGTTGGTG‑AAAAAAGTGTAACGCGTGCGCAGATTGCGCT‑GAAACCAGCAGGGCGGCTTGCATTGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTTCGCTGGCATTATCCAGATCAGGAGGTACggggt                                                                                               >  1:108111/1‑249 (MQ=255)                                                                                                                                                         cGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGAAAA‑GGG‑TTGGACCACATCCTTGGGCA‑TGCAGGGGACAAG‑AGGGGGATGTTGGGG‑GAAAAGGTGTAACGCGTGCGCCGATTGCGCT‑GAACCCAGCAGGGCGACGTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGGCCGTGGCTTCCCTCCGACGNCATTATCCAGATGAGGTGGTACGGGTAtttt                                                                                           >  2:232850/1‑250 (MQ=255)                                                                                                                                                           aTAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTtctc                                                                                         <  2:411825/251‑1 (MQ=255)                                                                                                                                                              aCCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGGT‑TT‑CCCCACATCATTGGTCA‑TGCAGTGGACCAG‑AGGGGGATGTTGGTG‑AAAAAAGGGTAACGCGGGGGCCGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAAAGTGAAAAAGCAATGACCTGGGGGTCCGTGAATTTCCTACGCTGGCATTATCCAGATCAGGGGATACGGGGATTTCTCAGc                                                                                      >  1:315234/1‑251 (MQ=255)                                                                                                                                                                      tGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGG‑TTTCACCCCATCATTTGGCA‑TGCAAGGAACACG‑AGGGGGATGGTGGGG‑GAAAAGGTGTAACGCGGGCGCAGAATTCGGT‑GAACCCCGCCGGGCGGATTGCCTTGTTTGCTCCCTCCATAACCGGAAGAAACAATGACCTGGTGGGCCGTGGCTTCCCTACCCCGGCCTTATCCAGATCAGGGGGTACGGGGTTTTTTCAGCCtttccgc                                                                              >  1:214766/1‑246 (MQ=255)                                                                                                                                                                                   ggcgcCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCAcccc                                                                 <  2:425394/250‑1 (MQ=255)                                                                                                                                                                                    gcgcCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGt                                                                                                                                                                                                                                                  <  2:49966/75‑1 (MQ=255)                                                                                                                                                                                       ccAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGt                                                             <  1:146654/251‑1 (MQ=255)                                                                                                                                                                                           cAAGGTATTGGCGGTAAA‑GGG‑TTGCACCACATCATTGGTCA‑TGCAGGGCACAAA‑GGGGGAATGTTGGGGGAAAAAGG‑GTAACGCGTGCGCAAAATGCGCT‑GAACCCAGCAGGTCGACCTGTATAGTTTTCTCCTGCCCTAACGTGAAAAAGCAATGACCTGGGGGTCCGGGACTTCCCTACGCTGGCATTATCCAGATCAGGGGATACGGGGATTTTTCAGCCTTTACGCAGAAGGGCACCCCGAGTCGtt                                                         >  1:461296/1‑251 (MQ=255)                                                                                                                                                                                                    ggCGGGAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCTAGTCGTTTGGTTGCg                                                 <  2:315234/250‑1 (MQ=255)                                                                                                                                                                                                       ggTAAA‑GGA‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGCCTTCCCTACGCTGGCATTATCCAGCTCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAACGGCACCCCGCGTCGCTTGGTTGCgatg                                              >  2:265432/1‑250 (MQ=255)                                                                                                                                                                                                                  ttGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGAGCAGATTGCGCT‑GCAACCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGCGTCGTTTGGTTGAGATGCTAAATCgac                                    >  2:402227/1‑249 (MQ=255)                                                                                                                                                                                                                         acaTCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGcggagg                             >  2:198546/1‑251 (MQ=255)                                                                                                                                                                                                                             catTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATgg                         >  1:21862/1‑251 (MQ=255)                                                                                                                                                                                                                                 ggTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAAt                     <  1:112974/251‑1 (MQ=255)                                                                                                                                                                                                                                   tCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGt                   >  1:212512/1‑251 (MQ=255)                                                                                                                                                                                                                                   tCA‑TGCAGTGAACAAG‑AGGGGAATGTTGGTG‑AAAACGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGTGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGt                   >  2:471366/1‑251 (MQ=255)                                                                                                                                                                                                                                    cA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCCATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTc                  >  1:72880/1‑251 (MQ=255)                                                                                                                                                                                                                                         cAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAgc              >  2:379857/1‑251 (MQ=255)                                                                                                                                                                                                                                         cAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAgc              <  1:340373/251‑1 (MQ=255)                                                                                                                                                                                                                                         cAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAgc              <  2:68963/251‑1 (MQ=255)                                                                                                                                                                                                                                               acaAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCAcgc        <  2:69813/251‑1 (MQ=255)                                                                                                                                                                                                                                                  ag‑agGGGATTGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCAAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGtt     >  2:287891/1‑251 (MQ=255)                                                                                                                                                                                                                                                   g‑agGGGAATGTTGGTG‑AAAAAGGTTTAAAGCGTGCGCAGATTGCGCT‑GAACCCAGGAGAACGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCAGTGACTTCCCGACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGAACCCAGCGTCGTTTGGTTGCGATGATCAATCGCCGAGGAGTATGGTAATGTCCAGCGCACGCGTTg    >  2:277944/1‑251 (MQ=255)                                                                                                                                                                                                                                                      ggggAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTa  <  1:208980/251‑1 (MQ=255)                                                                                                                                                                                                                                                      ggggAATGTTGGTG‑AAAAAGGAGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGCTGGAAATGTACAGCGCACGCGTTGt   >  2:316629/1‑250 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         CACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAA‑GGT‑TTGCACCACATCATTGGTCA‑TGCAGTGCACAAG‑AGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCT‑GAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTA  >  NC_000913/2184993‑2185484

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