Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,846,683	0	G	T	33.3%	45.2 / 12.5	45	Q244K (CAA→AAA)	hycD	formate hydrogenlyase subunit HycD
Reads supporting (aligned to +/- strand):  ref base G (7/23);  new base T (15/0);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 4.95e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCGGCGAGATAATGTTGACCTAATTTTTCTTCAGACATGCTCAAACTCTCTTTAATCACGCCGCCAGCAAGGAGACGAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCT‑GCAGCTCCTGCTCGGCTTCCGCCAGGTCGAACGGCAGTTTGCCCATTTCGATAAAGGTGGCGAACGCACAGGCACAAAGCGCCAGTACCAGCGGGATGCTCTGGCTCAGCGGCCAGTGATAAACGGTGTCGGT  >  NC_000913/2846434‑2846932                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                               gcggcgAGATAATGTTGACCTAATTTTTCTTCAGACATGCTCAAACTCTCTTTAATCACGCCGCCAGCAAGGAGACGAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTGt                                                                                                                                                                                                                                                               <  2:120842/251‑1 (MQ=255)      gagaTAATGTTGCCCTAATTTTTCTTCAGACATGCTCAAACTCTCTTTAATCACGCCGCCAGCAAGGAGACGAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTGTCCCCa                                                                                                                                                                                                                                                          <  1:1648/251‑1 (MQ=255)        gaTAATGTTGACCTAATTTTTCTTCAGACATGCTCAAACTCTCTTTAATCACGCCGCCAGCAAGGAGACGAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAAG‑TTTCCATTTGTCCCCAC‑g                                                                                                                                                                                                                                                       >  1:247536/1‑251 (MQ=255)               ttGACCTAATTTTTCTTCAGACATGCTCAAACTCTCTTTAATCACGCCGCCAGCAAGGAGACGAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGGG‑TAGG‑GTTCCCTTTTTCCCCCC‑GGACTaa                                                                                                                                                                                                                                                 >  2:14651/1‑250 (MQ=255)                       aTTTTTCTTCAGACATGCTCAAACTCTCTTTAATCACGCCGCCAGCAAGGAGACGAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGGG‑AAGGATTTCCTTTTG‑CCCCACCGGAA‑AAAAACCCcg                                                                                                                                                                                                                                         >  1:159023/1‑250 (MQ=255)                          tttCTTCAGACATGCTCAAACTCTCTTTAATCACGCCGCCAGCAAGGAGACGAACGCTAAAAATGCAAAGCCAAACCCAACCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGAGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGGGTAAGG‑TTTC‑ATTTTTccccccaggatcaccaaccccccac                                                                                                                                                                                                                                      >  2:424853/1‑230 (MQ=255)                                cAGACATGCTCAAACTCTCTTTAATCACGCCGCCAGCAAGGAGACGAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCt                                                                                                                                                                                                                               <  2:179901/251‑1 (MQ=255)                                       gCTCAAACTCTCTTTAATCACGCCGACAGCAAGGAGACGAACGCTAAAAATGCAAATCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAAGCACCAGCGCCATCAGCAGTCCACCGTCGGTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAa‑ca                                                                                                                                                                                                                        <  1:293660/250‑1 (MQ=255)                                           aaaCTCTCTTTAATCACGCCGCCAGCAAGGAGACGAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGGAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACAAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑GTGCCATTTTTCCCCAC‑GACATAAACACCCccacgcgaaactgcacaaccccc                                                                                                                                                                                                                    >  2:16171/1‑233 (MQ=255)                                                 tctTTAATCACGCCGCCAGCAAGGAGACGAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTGt                                                                                                                                                                                                                                                               <  1:39697/203‑1 (MQ=255)                                                                       ggAGACGAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTTCATTTGTCCCCAC‑GGCATAAACACCCCCCCGAACATCTGCAA‑CACCACCAGC‑TGTTTCC‑GGCTGGTACCCCCTTTcc                                                                                                                                                                                     >  1:55899/1‑250 (MQ=255)                                                                           aCGAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGAc                                                                                                                                                                                 <  2:247536/251‑1 (MQ=255)                                                                             gAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGCTAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCCGCAGCAGTCCACCGGCGGAG‑GAGG‑GTTTCATTTTTTCCCCC‑CGAAAAAAAACCCCGCCGAACATCTGCAA‑CCCCACCACA‑TGTTTCA‑CGGTGGTCCCCCATTTCATGACGc                                                                                                                                                                               >  1:64337/1‑251 (MQ=255)                                                                              aaCGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGcc                                                                                                                                                                              <  2:106240/251‑1 (MQ=255)                                                                              aaCGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAAGGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCAGTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTGCA‑GGCTGATACCCCATTTCATGACGcc                                                                                                                                                                              <  1:300763/251‑1 (MQ=255)                                                                               aCGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑TAGG‑GTTCCCTTTTTCCCCCC‑GGACTAAACACCCCCCCCAACATCTGCCA‑CACCACCACC‑CGTTTCA‑CACTGATAccccccttacttacaccc                                                                                                                                                                             >  2:166109/1‑237 (MQ=255)                                                                                     aaaaTGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGCG‑AAGG‑GTTCCCTTTTTCCCCAC‑CGAAAAAAAACCCCCACGAACAACTGCAAACACC‑CCACCCGGTTTC‑‑GGCCGATACCCCATTTCATGACGCCCAAgccc                                                                                                                                                                       >  2:382307/1‑250 (MQ=255)                                                                                          gCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAgccgctgcc                                                                                                                                                                  <  2:354073/251‑1 (MQ=255)                                                                                              aGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCCGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑GGTCCATTTTTCCCCCC‑GGGATAAACACCCCGACCGACATCTGCAA‑AACCACCCGC‑TGTTTTA‑GGCGGAGACCCCCTTTCAtatcccccaaaccccccgcgctg                                                                                                                                                              >  1:216056/1‑229 (MQ=255)                                                                                                               gTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGGATAAAAACCCCGCCCCACATCTTCAA‑CACCACCAGC‑CGTTTTC‑GGCTGATACCCCATTTCATGACGCCCAAACCCCCGCCGCGGTATTTAGAAAGCGGGcc                                                                                                                                             >  1:237003/1‑251 (MQ=255)                                                                                                                tAATGCGCGGAGTATTATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑GTTCCCTTTTTCCCCCC‑GGGATAAAAACCCCCCCGAACATCTGCCA‑CACCCCCCGC‑TGTTTCA‑GGCTGATACCCCATTTCATTACCCCCAAGCCACCGCCGCCGTATTCAGAGAGCGGGCCt                                                                                                                                            >  1:356992/1‑251 (MQ=255)                                                                                                                tAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGGTTTCC‑TTTTTCCCCCC‑GGAATAAAAACCCCCACGAACATCTGCAA‑CACCACCCGC‑TGTTTCA‑GGCTGATAACCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAAAACGGGCCt                                                                                                                                            >  1:83337/1‑251 (MQ=255)                                                                                                                    gcgcGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTcc                                                                                                                                         <  1:16342/250‑1 (MQ=255)                                                                                                                           gTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑GTTCCATTTTTCCCCCC‑CGAATAAACACCCCCACGAACATCTGCCA‑CACCCCCAGC‑CGTTTCA‑GGCTGATACCCCATTTCCTGACGCCAAAACCCCTGCCCCCGTATTCAGAGAGCCGGCCTTCCC‑GCAGccc                                                                                                                                >  1:28973/1‑249 (MQ=255)                                                                                                                                    aGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCT‑GCAGCTCCTGCTCGGc                                                                                                                       <  1:402783/251‑1 (MQ=255)                                                                                                                                        gCAGCCGGGCCATTCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCT‑GCAGCTCCTGCTCGGCTTcc                                                                                                                   <  2:398385/251‑1 (MQ=255)                                                                                                                                        gCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGGTTTCC‑TTTTTCCCCCC‑GGGATACACACCCCGCCGAACATCTTCAA‑CACCACCAGC‑CGTTTCA‑GGCTGATACCCCCTTTCATTACGCCAAAACCGCCGCCCCCGTATTCAGAAAACGGGCCTTTCTCGCA‑CCCCTGCCCGGCTTcc                                                                                                                   >  1:421875/1‑251 (MQ=255)                                                                                                                                        gCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCCTTTTTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAAC‑TGTTTCA‑GGCTGATACCCCATTTCATGAACCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCT‑GCCACTCCTGCTCGGCTTcc                                                                                                                   >  1:81459/1‑251 (MQ=255)                                                                                                                                               ggCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAAG‑TTTCCCTTTTTCCCCCC‑CGAATAAACACCCCCACGAACATCTTCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAACCGCTTCCGCTGTATTCCGAGAACGGGGCTTCCT‑GCAGCTCCTGCTCGGCTTccccccggg                                                                                                            >  1:139707/1‑244 (MQ=255)                                                                                                                                                 ccATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCATCAGTCCACCGGCGGTG‑AAGT‑TTTCCATTTGTCCCCAC‑g                                                                                                                                                                                                                                                       <  2:292453/114‑1 (MQ=255)                                                                                                                                                 ccATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑GTTCCATTTTTCCCCAC‑g                                                                                                                                                                                                                                                       >  1:292453/1‑114 (MQ=255)                                                                                                                                                  cATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCT‑GCAGCTCCTGCTCGGCTTCCGCCAGGTCGa                                                                                                         <  2:114753/251‑1 (MQ=255)                                                                                                                                                  cATGCTATTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTTTCCCCAC‑AGGATAAACACCCCGACGAACATCTGCAA‑AACCACCAGC‑TGTTTCA‑GGCTGATAACCCCTTTCCTGAACACAAAACCGCCGCCGCTGTATTCAGAGAGCGGGGCCTCCT‑GCCACTCCTGCTCGGCTTCCCCCCGGTCgc                                                                                                         >  1:352761/1‑250 (MQ=255)                                                                                                                                                        gTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCT‑GCAGCTCCTGCTCGGCTTCCGCCAGGTCGAACgg                                                                                                     <  2:139707/249‑1 (MQ=255)                                                                                                                                                         ttttCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑GTTTCATTTGTCCCCCC‑GGGATAAACACCCCCACGGACATCTGCAA‑CACCACCCAC‑TGGTTCAAGCCTG‑TAACCCCTTTCATGACGCCAAAGCCGCTGCCCCTGTATTCATAGAGCGGGGCTTCCC‑GCAGCTCCTGCTCCGCTTCCCCCAGGTCGAACGGcgg                                                                                                  >  1:203775/1‑249 (MQ=255)                                                                                                                                                             cGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCTGCGGTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGTCTTCCT‑GCAGCTCCTGCTCGGCTTCCGCCAGGTCGAACGGCAGTTTg                                                                                              <  2:149066/251‑1 (MQ=255)                                                                                                                                                                        aTAACCGGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCTGCGGTG‑AAGG‑TCTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGTCTTCCT‑GCAGCTCCTGCTCGGCTTCCGCCAGGTCGATCGGCAGTTTGCCCATTTctat                                                                                   <  2:91874/251‑4 (MQ=255)                                                                                                                                                                            ccAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCT‑GCAGCTCCTGCTCGGCTTCCGCCAGGTCGAACGGCAGTTTGCCCATTTCGATAAAg                                                                               <  1:114075/251‑1 (MQ=255)                                                                                                                                                                                                    cGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCT‑GCAGCTCCTGCTCGGCTTCCGCCAGGTCGAACGGCAGTTTGCCCATTTCGATAAAGGTGGCGAACGCACAGGCACAAAgc                                                       <  1:14651/251‑1 (MQ=255)                                                                                                                                                                                                        ggCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCT‑GCAGCTCCTGCTCGGCTTCCGCCAGGTCGAACGGCAGTTTGCCCATTTCGATAAAGGTGGCGAACGCACAGGCACAAAGCGcc                                                    <  1:382307/250‑1 (MQ=255)                                                                                                                                                                                                                               gTCCACCGGCGGTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCT‑GCAGCTCCTGCTCGGCTTCCGCCAGGTCGAACGGCAGTTTGCCCATTTCGATAAAGGTGGCGAACGCACAGGCACAAAGCGCCAGTACCAGCGGGATGCTCTGGCTc                            <  1:268376/251‑1 (MQ=255)                                                                                                                                                                                                                                         ggTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCT‑GCAGCTCCTGCTCGGCTTCCGCCAGGTCGAACGGCAGTTTGCCCATTTCGATAAAGGTGGCGAACGCACAGGCACAAAGCGCCAGTACCAGCGGGATGCTCTGGCTCAGCGGCCAGt                  >  1:421230/1‑251 (MQ=255)                                                                                                                                                                                                                                           tG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCT‑GCAGCTCCTGCTCGGCTTCCGCCAGGTCGAACGGCAGTTTGCCCATTTCGATAAAGGTGGCGAACGCACAGGCACAAAGCGCCAGTACCAGCGGGATGCTCTGGCTCAGCGGCCAGTg                 >  2:228589/1‑250 (MQ=255)                                                                                                                                                                                                                                           gg‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCT‑GCAGCTCCTGCTCGGCTTCCGCCAGGTCGAACGGCAGTTTGCCCATTGCGATAAAGGTGGCGAACGCACAGGCACAAAGCGCCAGTACCAGCGGGATGCTCTGGCTCAGCGGCCAGTGa                <  2:55899/250‑1 (MQ=255)                                                                                                                                                                                                                                                 g‑TTTCCATTTGTCCCCAC‑GGAATAAACTCCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACTCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCT‑GCAGCTCCTGCTCGGCTTCCGCCAGGTCGAACGGCAGTTTGCCCATTTCGATAAAGGTGGCGAACGCACAGGCACAAAGCGCCAGTACCAGCGGGATGCTCTTGCTCAGCGGCCATTGATAAAc           <  1:312503/251‑1 (MQ=255)                                                                                                                                                                                                                                                            gTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCT‑GCAGCTCCTGCTCGGCTTCCGCCAGGTCGAACGGCAGTTTGCCCATTTCGATAAAGGTGGCGAACGCACAGGCACAAAGCGCCAGTACCAGCGGGATGCTCTGGCTCAGCGGCCAGTGATAAACGGTGTCGGt  <  2:28973/250‑1 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                               GCGGCGAGATAATGTTGACCTAATTTTTCTTCAGACATGCTCAAACTCTCTTTAATCACGCCGCCAGCAAGGAGACGAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGG‑TTTCCATTTGTCCCCAC‑GGAATAAACACCCCGACGAACATCTGCAA‑CACCACCAGC‑TGTTTCA‑GGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCT‑GCAGCTCCTGCTCGGCTTCCGCCAGGTCGAACGGCAGTTTGCCCATTTCGATAAAGGTGGCGAACGCACAGGCACAAAGCGCCAGTACCAGCGGGATGCTCTGGCTCAGCGGCCAGTGATAAACGGTGTCGGT  >  NC_000913/2846434‑2846932

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