Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,174,605	0	T	C	25.0%	127.1 / 10.1	57	G187G (GGA→GGG)	ycfT	inner membrane protein YcfT
Reads supporting (aligned to +/- strand):  ref base T (35/7);  new base C (0/14);  total (35/22)
Fisher's exact test for biased strand distribution p-value = 2.00e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCAGGCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCA‑GGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCA  >  NC_000913/1174356‑1174831                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                     cagcagGCTGGGGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGTCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCAtt                                                                                                                                                                                                                                      <  2:163040/250‑1 (MQ=255) cagcagGCGGGTGGAGCGCATACCGAATCGTTTCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCAtt                                                                                                                                                                                                                                      <  2:426902/250‑1 (MQ=255)  agcagGCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACACCAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GGTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTcc                                                                                                                                                                                                                                    <  2:38146/251‑1 (MQ=255)       gcgggtggAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAAGAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCccacc                                                                                                                                                                                                                               <  1:256771/248‑1 (MQ=255)          ggtggAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACg                                                                                                                                                                                                                            >  1:325332/1‑251 (MQ=255)                  gcATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAAGCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCCCCACGGC‑GTGGGa                                                                                                                                                                                                                   <  1:305967/251‑1 (MQ=255)                     tACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACg                                                                                                                                                                                                                >  2:15755/1‑251 (MQ=255)                       ccGAATCGTTGCTCGTATTTTTAAAAGAGCTTTATGAGCACCACAACCGATACCAGCGACAACAACAGCGAGATATTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTGAAC‑A‑CAGGTCATTATG‑GTTGCGCCATACC‑AGGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCGCGGC‑GTGGTAACGaa                                                                                                                                                                                                              <  2:424898/251‑1 (MQ=255)                         gAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TNTTCATTCCCCACCACGGC‑GTGGGAACGAAAt                                                                                                                                                                                                            >  2:208425/1‑251 (MQ=255)                         gAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAAt                                                                                                                                                                                                            >  1:150765/1‑251 (MQ=255)                              gTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAAcc                                                                                                                                                                                                       >  1:162993/1‑251 (MQ=255)                                tGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑TTGGGAACGAAATTAACCGc                                                                                                                                                                                                     >  1:22017/1‑251 (MQ=255)                                       aTTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTc                                                                                                                                                                                              >  2:54659/1‑251 (MQ=255)                                                    cTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTacaaac                                                                                                                                                                                  >  1:114051/1‑250 (MQ=255)                                                        aTGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGc                                                                                                                                                                             >  1:51395/1‑251 (MQ=255)                                                          gATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCTAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCg                                                                                                                                                                            >  1:155190/1‑250 (MQ=255)                                                                   aaTCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGc                                                                                                                                                                   >  1:431092/1‑250 (MQ=255)                                                                               cgacaccggaggcGAGATAGAAAACAACAAGGCACCGACCGCAAGAACGGAAAGCAAAGAAGCCATCAGCAGAGGGCGGCGCAACGGCCCCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCAAACCC‑AGGCGCCAAGGCTGTAATAAGG‑CAAA‑TG‑GCGGATCACAC‑TGTTCATCCCCCCCCCCGGC‑GGGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGc                                                                                                                                                      <  1:29440/240‑1 (MQ=255)                                                                               cGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACNAACAAGGCGAATAGTGGCAGCGCCAGACGGc                                                                                                                                                      >  2:212780/1‑251 (MQ=255)                                                                                gACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGCCGGCt                                                                                                                                                     >  2:361185/1‑251 (MQ=255)                                                                                cacaacaGCAGCGAGATAGAAAACACCCAGGCACCGCCCGCAAAAACGGTAAGCAAAGAAGCCACCAGCAGATGGCGGCGCACCGCCACCT‑TTTAACCA‑CAGGTCTTTATT‑GTTCCGCCAAACC‑ATCCGCCAAGGCTGTAAAAAGG‑AAAA‑TT‑GCGAATCACCC‑TGTTCATCCCCCCCCCGGGC‑GTGGGAAGAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCt                                                                                                                                                     <  1:146640/250‑1 (MQ=255)                                                                                    acaGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGACATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaa                                                                                                                                                  >  2:305701/1‑250 (MQ=255)                                                                                       gcagcGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtt                                                                                                                                              >  1:332898/1‑251 (MQ=255)                                                                                        cggggagaaaaaaaacaccaagcccccgaccgCCAGAAGGGAAAGCAAAAAAGCCAACAGAAGAGGGCGGCACACCGCCACCTTTATACC‑A‑CAGTTCATTTTT‑GTTGCGCAAAACC‑AGGCGCAAAGGCTGTAATAAGG‑CAAA‑GT‑GCGGACCACAC‑TCTTCATCCCCCACCAGGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAttt                                                                                                                                             <  2:244140/226‑1 (MQ=255)                                                                                        cagcGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAttt                                                                                                                                             >  2:378774/1‑251 (MQ=255)                                                                                         agcGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGCACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtttt                                                                                                                                            >  1:76283/1‑251 (MQ=255)                                                                                             agaTAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACAt                                                                                                                                        >  1:161293/1‑251 (MQ=255)                                                                                               aTAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGTAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAACTTTTACATAc                                                                                                                                      >  1:465800/1‑251 (MQ=255)                                                                                               aTAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGAAACGGCACCTCTTTAAC‑A‑CCGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTAAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAc                                                                                                                                      >  1:465820/1‑251 (MQ=255)                                                                                                 agaaaaCACCCAGGCCCCGACCGCCAAAACGGTAAGCAAAGAAGCCACCAGCAGATGGCGGCGCAACGGCACCTCTTAACC‑A‑CAGGTCTTTTTT‑GTTGCGCCAACCC‑ATGCGCCAAGGCTGAAAAAAGG‑CAAA‑TG‑GCGGACCCCAC‑TTTTCTTCCCCCCCCCCGGC‑GGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCa                                                                                                                                    <  1:378774/248‑1 (MQ=255)                                                                                                      acaacCAGCCCCCGCCCGCCAGAACGTTAAGCAAAAAAGCCATCAGCAGATGGCGGCGCACCGGCCCCTCTTTAAC‑A‑CAGGTTATTATT‑GTTGCGCAAAACC‑ATGCGCCAAGGCTGTAAAAAGG‑AAAA‑TT‑GCGGATCACAC‑TGTTCATCCCCCCCCAGGGC‑GGGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGaaa                                                                                                                               <  1:208425/251‑1 (MQ=255)                                                                                                        aacCAGGCACCGACCGCCAGAACGGTAAGCAAAGCAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACgcaaa                                                                                                                              >  2:360963/1‑246 (MQ=255)                                                                                                                        caaaaagggaaaacAAAGAACCCACCGGCAGGTGGCGGCGCACCGGCCCCTCTTAAAC‑A‑AAGGT‑ATTATTTGTTGCGCAAACCC‑ATGCGCCAGGGCTGTAATAAGG‑CAAA‑TG‑GCGGACCACCT‑TGTTCTTCCCCCACCACGGC‑GGGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATAc                                                                                                             <  2:401636/238‑1 (MQ=255)                                                                                                                        caaGAACGGAAAGCAAAGAACCCATCAGCAGAGGGCGGCGCAACGGCCCCTCTTAAAA‑A‑CAGGTCATTATT‑GTGGCGCCAAACC‑ATGCGCCAAGGCTGTAAAAAGG‑AAAA‑TT‑GCGGATCACAC‑TGTTCACCCCCCCCCACGGG‑GGGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATAc                                                                                                             <  2:161293/249‑1 (MQ=255)                                                                                                                                    gCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCAACTCTTTAAC‑A‑CAGGAATTTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTAATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCCGACGGCTAAAAATTTTACATACCACGAAAAAGACAATTAACGCCTACATATACCACaac                                                                                                 >  2:464879/1‑250 (MQ=255)                                                                                                                                     aaaagacgccaccagtagtgggcggcGAAACGGCACCT‑TTTAAC‑AACAGGTTTTTTGT‑GTGGCGCAAAACC‑ATGCCCCAGGGCTGTAA‑AAGGTCAAT‑TT‑CCGGATCACTC‑TTTTCTTCCCCCCCCACGGC‑GGGGAAACAAATTTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGc                                                                                                <  1:373255/231‑1 (MQ=255)                                                                                                                                          aagtcatcacatgtgggctgcccAACGCCACCTCTTAAAA‑A‑CAGGTTATTATT‑GTTGCGCCAATCC‑CTGCGCCAAGTTTG‑AAAAAAG‑CAAA‑TTTGCG‑ATCCCCCCTGTCCATCCCCCCCCACGGC‑GTGGAAACGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTATCGCATACAGATACCACAAGCTGGTg                                                                                           <  2:325332/229‑1 (MQ=255)                                                                                                                                              cATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTgg                                                                                       >  1:129647/1‑251 (MQ=255)                                                                                                                                                  agcggATGGCGGCCCAACGCCCCCTTTTTAAC‑C‑CAGGTAATTTTT‑GTGGC‑CCAAACCCATCCCCCAGGGCTGTAATAAGG‑CAAT‑TG‑GCGGATCCCCC‑TGTTCATTCCCCCCCTCGGC‑GTGGGAAGAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGt                                                                                   <  1:56608/247‑1 (MQ=255)                                                                                                                                                     agagggcgcgacagcgtcccctttttAAC‑A‑CNGGACATTATG‑GTTGCGCCAAACC‑ACGCGCAAAGGCTAAAAGAAA‑‑CAAAGTG‑GCGGATCCCAC‑TTTTCCTCCCCCACCACGGG‑GGGGAAAAAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGat                                                                                <  2:204858/228‑1 (MQ=255)                                                                                                                                                         ggcggcGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCAAACAGATACCACAAGCGGGTGCTGGCGGCGatcatc                                                                            >  1:207870/1‑251 (MQ=255)                                                                                                                                                            ggcgCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCCCAAGCTGGTGCTGGCGGTGATCAAcccc                                                                         >  1:212281/1‑250 (MQ=255)                                                                                                                                                                    ggCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCCTGCA‑GGa                                                                 >  1:361057/1‑250 (MQ=255)                                                                                                                                                                      cACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCA‑GGAAcc                                                              >  1:136902/1‑250 (MQ=255)                                                                                                                                                                       ccctttttAAC‑A‑CAGTTCATTATT‑GTTGCGCCACACC‑AGCCCCCAGGGCTTTAATAAGG‑CAAA‑TT‑GCGAACCCAAC‑TGTTCATTCCCCCCCCCGGC‑GTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCA‑GGAACTc                                                             <  1:54659/246‑1 (MQ=255)                                                                                                                                                                               taaacaagggTCATTATG‑GTTGCCCCAAACA‑AGGCCCCAAGGGTTAAAAAAGG‑AAAA‑TT‑GAGGACCCCAC‑TGTTCCTCCCCCCCCACGGC‑GGGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCA‑GGAACTCAcc                                                          <  2:332898/240‑1 (MQ=255)                                                                                                                                                                                         gtatttaggttggcgcCAAACC‑AGGCCCCAGGGCTGAAAAAAGG‑AAAA‑GT‑GGGCATCCCCC‑GGTTCAT‑CCCCCCCCCGGCCGTGGAAACGAATTTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCA‑GGAACTCACCGGTGGAATCGGc                                              <  2:134994/239‑1 (MQ=255)                                                                                                                                                                                          atttttt‑gttgCGCTAACCC‑ACGCGCAAAGGCTGTAATAAGG‑CAAA‑TG‑GCGGACCCCCC‑TGTTCATCCCCCACCGCGGG‑GGGGAAAAAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCA‑GGAACTCACCGGTGGAATCGGCa                                             <  2:155190/245‑1 (MQ=255)                                                                                                                                                                                           attatt‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑AAAA‑TT‑GCGGATCACAC‑TTTTCTTCCCCCCCCGCGGG‑GTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCA‑GGAACTCACCGGTGGAATCGGCAt                                            <  1:361185/250‑1 (MQ=255)                                                                                                                                                                                            ttatt‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCCTCCCGTGCA‑GGAACTCACCGGGGGAATCGGCATa                                           >  1:282191/1‑250 (MQ=255)                                                                                                                                                                                             tatt‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATTCCGGGCAAGGCAC‑CACCGGCGGGATCGGCATaa                                          >  2:25822/1‑250 (MQ=255)                                                                                                                                                                                                       gcCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCA‑GGAACTCACCGGGGGGATCGGCATAAGCGGCATTgg                                >  1:278031/1‑251 (MQ=255)                                                                                                                                                                                                        caaaaCC‑AGGCGACAAGGCTATAAAAAGG‑AAAA‑TT‑GCGTATCCCAA‑TTTTCAGTCCCCCCCGTGGC‑GGGTAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCA‑GGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGa                               <  1:307805/249‑1 (MQ=255)                                                                                                                                                                                                         cAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGAAAATCCCGTGCA‑CGAAATCACCGGGGGAAACGGGATaaaggacatttggg                              >  2:341871/1‑239 (MQ=255)                                                                                                                                                                                                                     ccAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCA‑GGAAATCCCCCGGGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTa                    >  1:397692/1‑250 (MQ=255)                                                                                                                                                                                                                           cTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCCTGCA‑GGAACTCAACGGTGGAATCGGCATAAGCGGCATTGGAGGCATTTCTTAAATcgcg             >  1:409543/1‑251 (MQ=255)                                                                                                                                                                                                                             ccgatcggg‑CAAA‑TG‑GCGGACCACAC‑TGTTCATCCCCCCCCCCGGC‑GGGGGAACAAAATTACCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCA‑GGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCt           <  2:173059/244‑1 (MQ=255)                                                                                                                                                                                                                                    gg‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGGGCA‑GGAAATCACCGGTGGGATCGGCAAAAGCGGCATTGGAGGCATTTCTTAAATCGCGCt           >  1:173059/1‑244 (MQ=255)                                                                                                                                                                                                                                        aaa‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGGA‑GGAACTCACCGGGGGAATCCGCATAAGCAGCATTAGAGGCATTGCTTAAATCGCGCTCAGGTgccc  >  1:58139/1‑249 (MQ=255)                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                     CAGCAGGCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAAC‑A‑CAGGTCATTATT‑GTTGCGCCAAACC‑ATGCGCCAAGGCTGTAATAAGG‑CAAA‑TT‑GCGGATCACAC‑TGTTCATTCCCCACCACGGC‑GTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCA‑GGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCA  >  NC_000913/1174356‑1174831

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