Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,008,355	C→T	19.9%	L27F (CTC→TTC)	yedE →	UPF0394 family sulphur transport domain‑containing inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,008,355	0	C	T	19.9%	85.0 / 4.3	48	L27F (CTC→TTC)	yedE	UPF0394 family sulphur transport domain‑containing inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (11/26);  new base T (8/2);  total (20/28)
Fisher's exact test for biased strand distribution p-value = 8.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTAAAGCCCTGATGTCGCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGTTACTTTAAAATTATCCATCTGG‑AAGG‑ATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGC‑TTTGCCGCAGCGCTGTGGGCCAACAATGTCAA  >  NC_000913/2008127‑2008581                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            aTTAAAGCCCTGGTGGTCCCAGCCTTTCAAGGTTAGGAGGAAAATTTTATCTAACTATTATTTTTGTGAATTCTTACCCGGAAAATAACTTTGCATGGTGTTTGATTCGCTGCTCGTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCAGGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCtactt                                                                                                                                                                                                                                    <  2:1057446/238‑1 (MQ=255) aTTAAAGCCCTGATGTCGCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCtactt                                                                                                                                                                                                                                    <  2:672498/238‑1 (MQ=255)                  gCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTTTCTTCTTACTC‑TTTTTGCAttcttg                                                                                                                                                                                                                 >  1:427014/1‑235 (MQ=255)                   cccgctttttAAGGTTTAGAGAGAAATTTTATTTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTgg                                                                                                                                                                                                                <  2:376611/233‑1 (MQ=255)                    cAGCCTTTCAAGGTTCAAGGAGAAATTTGATCTAATTATGATTTTGGGGAATTATTACCCCGAAAATAAATTTGCATGGTGTATTATTCGCGCCTCTTTTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCATCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTGCTGGc                                                                                                                                                                                                               <  1:166176/239‑1 (MQ=255)                                       aagaAATTTGATCTAACTTTTTTTTTTGTAAATTATTCCCCCAAAAATAAATTTGCAGGGTGTATTATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCATAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCCGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACgg                                                                                                                                                                                          <  2:427014/238‑1 (MQ=255)                                                                  tggAATTTTTTCCCCGAAAATAAATTTGCATGGTGTTTGGTTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTAAATTCAGGTCAGGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCCCC‑TTTTGGGCTG‑GCACGGGTGAATTTACCCGTTGGGGCGGCCAGc                                                                                                                                                               <  2:126459/237‑1 (MQ=255)                                                                  gtgAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAATTTCTGGGCCCCCATCCTCGGGGTCATCGCGGAGGGTTTTTTCTTCTCTTCACC‑TTTTTGCATTCCACGCCCCTTTTTGCGCT‑‑TCACCGGTGAATTTACACGTTGGgggggcaga                                                                                                                                                                >  2:637901/1‑235 (MQ=255)                                                                  gtgAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGTTTTTCTCCCTTTCATATTTTTGGCTTTCCGGGC‑CC‑TTTTGGGGTT‑TCACGGGGGAATTTATCCCTTGgggggggcga                                                                                                                                                                >  2:960071/1‑234 (MQ=255)                                                                   acagttatttCCCCAAAAAAAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAACCCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGc                                                                                                                                                               <  1:419318/234‑1 (MQ=255)                                                                                  aaaaTAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGGTTTTCTCTCCTTACTC‑TTTTTGCCTTTATGTCACC‑CTTTTGGCGT‑TCACGGGTGAATTTACACCTTGGGGGGGCCCACTcccccgaaagagg                                                                                                                                                 >  2:419318/1‑226 (MQ=255)                                                                                     ataaatTTGCATGGGGTATGTTGCGCGTCTCTTGTTTAAAACCTTTTGAGTTAAGATGTCACATGCATGGCAGGGCAGCTATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGt                                                                                                                                           <  2:777991/239‑1 (MQ=255)                                                                                       aaatTTGCATGGTGTATTATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTcc                                                                                                                                         <  1:404088/239‑1 (MQ=255)                                                                                         atTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTTTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTTGGCTG‑TCCCGGGTGAATTTACCCGTTTGGGGGGCCAGGCCCTTCCACTGTTC‑GGCGTCCAt                                                                                                                                       >  1:1012735/1‑239 (MQ=255)                                                                                            ggctgggttttttgttCGGGACTTTTTTTTCAAACCTTCGTAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTGC‑GGCGTCCATG‑Ct                                                                                                                                   <  2:18041/225‑1 (MQ=255)                                                                                             gCATGGTGTATGGTTCGCGACTTTTGTTTCAACCCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTg                                                                                                                                  <  1:960071/239‑1 (MQ=255)                                                                                                        gtttCCCGATTTTTGTTTAAAACCTTCTGAGTTCAGAGGTTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGt                                                                                                                      <  2:249898/237‑1 (MQ=255)                                                                                                             gcgACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTTGCATTACTGGCCCC‑TTTTTGGCTG‑TACCGGGTGAATTTTACCGTTGGGGCGGCCAGCTCCTGCAAACTTTC‑GGGGTCCTTG‑CTGCAGAAG‑GGGGGTTCtt                                                                                                                 >  1:470161/1‑239 (MQ=255)                                                                                                             gcgACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTTGCATTACTGGCACC‑CTTTTGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTT‑GGGGTCCATG‑CTGAAGAGT‑GGGGGTACtt                                                                                                                 >  1:249898/1‑239 (MQ=255)                                                                                                             gcgACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGTCACC‑TTTTGGGGTT‑TCCCGGGGGGATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTTAAGAGT‑GGGGGGACtt                                                                                                                 >  1:297011/1‑239 (MQ=255)                                                                                                             gcgACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACCA‑TTTTGGCATTACTGGCACC‑TTTTTGGCTG‑TCACAGGTGAATTTACCCGTTGGGGCGGCCAACTCCTTCAACTTTTC‑GCCGTCCCTG‑CGGTGAAGT‑ggggcgggct                                                                                                                 >  1:18041/1‑233 (MQ=255)                                                                                                             gcgACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCTCTACTTACTA‑TTTTTGGATTACTTGCACC‑TTTTGGGCTG‑TCACGGGGGAATTTACCCGTTGGGGGGGCCCGACCCTTCAACTGTTC‑GGCGCTCATG‑CTGGAACGC‑GGGGTTaact                                                                                                                 >  1:735501/1‑236 (MQ=255)                                                                                                             gcgACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTTTCTACTTTCCA‑TTTTTGCCTTACTTGCACC‑TTTTTGGCTGTTCACGGG‑GAATTTTCCCCTTGGGGGGGCCCACTCCCGCCACTGTTC‑GGCGGCCATT‑CTGGAAAGT‑GGGGGTTCtt                                                                                                                 >  1:759803/1‑239 (MQ=255)                                                                                                                     tGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGCATTCTCTCTACTTACTA‑TTTTGGCCTTACTGGCACC‑TTTTGGGCTG‑TCACAGGTGAATTTACCCCTTTGGGGGGCCCGCGCCTGCAACCGTTC‑GGCGTCCATG‑CTGAACAGT‑TGGGGTACTTTAAAATTa                                                                                                         >  1:873783/1‑239 (MQ=255)                                                                                                                       tttCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGTTACTTTAAAATTATc                                                                                                       <  1:421773/239‑1 (MQ=255)                                                                                                                       tttCAAACCTTCTGAGTTCAGAGGCTACATTCATGGCATGGCAGCAATTCAAACGCGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGTTACTTTAAAATTATc                                                                                                       <  1:82920/239‑1 (MQ=255)                                                                                                                                    gAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGTTACTTTAAAATTATCCATCTGG‑AAGG‑At                                                                                        <  1:639496/239‑1 (MQ=255)                                                                                                                                    gAGTTCAGAGGCTACATTCATGCCATGGCAGCAATTCAAACTCGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGTTACTTTAAAATTATCCATCTGG‑AAGG‑At                                                                                        <  1:189492/239‑1 (MQ=255)                                                                                                                                            agGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATtct                                                                                                                                                                                                                                            <  2:707543/91‑1 (MQ=255)                                                                                                                                                 aCATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTACTTACTT‑TTTTTGCCTTACTGGCACC‑TTTTTGGCCG‑TCACGGGTGGATTTTACCGTTGGGGCGGGCAGCCCCTGCCACCGTTC‑GGCCGCCCTG‑CTTAAAGGTGGGGGTTACTTAAAATTTATCAATCGTG‑AGGA‑ATCGCC‑TTAAcccc                                                                           >  2:169329/1‑238 (MQ=255)                                                                                                                                                      catgtcatgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTACTTACTA‑TTTTGGCCTTACTGGCACC‑TTTTTGGCTG‑TCACGGGGGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑TGGGTCCATG‑CTGTAAAGG‑GGGGTTCACTTAAAATTATCCACCTGG‑AAGG‑ATCGCCCTTAAACCCCactc                                                                      >  1:892815/1‑236 (MQ=255)                                                                                                                                                       tttttctgGCAGCAATTCAAACCCGCTTTGTTGTTTTAATTCTTGGCCCCCATCCCTGCGGTCATCGCGGCGGGTTTTCTCTCTGCTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGTTACTTTAAAATTATCCATCTGG‑AAGG‑ATCGCCATTAACCCGCATCGa                                                                     <  2:470161/233‑1 (MQ=255)                                                                                                                                                       ttttcctgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGTCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTTAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGTTACTTTAAAATTATCCATCTGG‑AAGG‑ATCGCCATTAACCCGCATCGa                                                                     <  2:1012735/233‑1 (MQ=255)                                                                                                                                                       ggtataagagctggaTTCAAACGCGCCTGGTTGATTAAATTCTGGGCCCCCATCCCGGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCGCC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGTTACTTTAAAATTATCCATCTGG‑AAGG‑ATCGCCATTAACCCGCATCGa                                                                     <  1:876503/225‑1 (MQ=255)                                                                                                                                                       cgtgtattgcagcaATTAAAACCCGGCTTGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGGGGCGGGTATTCTCTCTACGTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGTTACTTTAAAATTATCCATCTGG‑AAGG‑ATCGCCATTAACCCGCATCGa                                                                     <  1:637901/231‑1 (MQ=255)                                                                                                                                                            atgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTTTCTACTTACTA‑TTTTGGGCTTACTGGCACC‑TTTTTGGCTG‑TCACCGGGGAATTTTCCCCTTGGGGCGGCCCGCCTCTTCACATGTTCTGGCGT‑CATGCCTAAAGAGG‑GGGGGCTTTTTAAAATT‑TCCCTCTGG‑GAGG‑GACGCCCATAAccccgcctccgcggg                                                                >  2:1022665/1‑227 (MQ=255)                                                                                                                                                               gcagcaATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTACTTACTC‑TTTTTGCATTACTTGCACC‑TTTTGGGCTG‑TCACGGGTGGATTTACCCGTTGGGGCGGCCAGCCCCTGCAAATGTTT‑GGGGGCCCTG‑CCTAAGGGG‑GGGGGTGCTTTTAAAATTTCCACTTGG‑GAGG‑ATGGCCCATAACCCCCaccccaggggggg                                                             >  2:267874/1‑227 (MQ=255)                                                                                                                                                                 agcaATTCAAACCCGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGTTACTTTAAAATTATCCATCTGG‑AAGG‑ATCGCCATTAACCCGCATCGACGGgatgat                                                            <  1:169329/238‑1 (MQ=255)                                                                                                                                                                   caATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑CTTTTGGCTG‑TCACGGGGGCATTTACCCGTTGGGGCGGCCCGCTCCTGCCACTTTTT‑GGCGTCCATC‑CTGAAGAGT‑GGGGTTACCTTTAAATTATCCATCTGG‑GAGG‑GTCGCCCTTAACCCCCATCGACGGggtgctgt                                                          >  1:214417/1‑237 (MQ=255)                                                                                                                                                                      ttCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGTTACTTTAAAATTATCCATCTGG‑AAGG‑At                                                                                        <  1:372072/205‑1 (MQ=255)                                                                                                                                                                                ccTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTTTCTCTACTTACTA‑TTTTGGCCTTACTTGCACC‑CTTTTGGGTG‑TCCAGGGGGAATTTTCCCGCTGCGGCCGCCCCCCCCCGCCACTGGTG‑GGCGCCCATG‑CTGATGACG‑GGGGGTTCTTTTAAATTTTTCCTTTGG‑TAGGCAGCG‑CATTATACCGCATCCTCGGGATGATGTCCCTCTGTATtttt                                            >  2:689504/1‑239 (MQ=255)                                                                                                                                                                                           aaaTTCTGGGCCCCCATCCCGGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCAAC‑TTTTGGGCTG‑GCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGTTACTTTAAAATTATCCATCTGG‑AAGG‑ATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGC‑tt                                <  1:1028428/239‑1 (MQ=255)                                                                                                                                                                                                   ggCCCCCATCCCTGCGGTCATCGCGGCGGGTATTTTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTTGGCTG‑TCCTGGGTGGATTTACCCGTTTGGGCGGCCACCCCCTGCCACTGTTC‑CGCGTCCATG‑CTTGAGGGG‑GGGGGTGCTTTTAAATTATCCTTCTGG‑AAGG‑ATCCCCCTTAACCCGCCACCACGGGGGGGTGGTCCCCGGGATGTTTTGGGGCTGCGTTTCCCGc                          >  2:129995/1‑238 (MQ=255)                                                                                                                                                                                                   ggCCCCCATCCCTGCGGTCATCGCGGCGGGGATTATCTCTACTTACTC‑TTTTCGCCTTACTGGCACC‑TTTTGGGCTG‑TTACGGGTGAATTTATCCGCTGGGGCGGCCAGCTCCTGCCACTGTTC‑GGCTGCCCTC‑CTGAAGCGT‑GGGGTTACTTTAAAATTATCCCTCTGGAAACG‑AT‑GCCATTAAACCCCTCCCACGGGGTGATGAACACCGGGATTTTTTGCGGGCGC‑TTTGCCGccg                        >  2:406990/1‑237 (MQ=255)                                                                                                                                                                                                     cccccA‑CCCTGCGGTCATCGCGGCGGGTATTCTTTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGTTACTTTAAAATTATCCATCTGG‑AAGG‑ATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGC‑TTTGCCGCAgcgc                     <  1:932884/239‑1 (MQ=255)                                                                                                                                                                                                      ccccATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGTTACTTTAAAATTATCCATCTGG‑AAGG‑ATCGCCATTAACCCGCATCGACggg                                                                 <  1:305807/196‑1 (MQ=255)                                                                                                                                                                                                              cTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCCCC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGTTACTTTAAAATTATCCATCTGG‑AAGG‑ATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGC‑TTTGCCGCAGCGCTGTGGGcc             <  1:267874/239‑1 (MQ=255)                                                                                                                                                                                                                gcgttcctcggtgcggGTTTTTTCTTTTCTTCCTT‑TTTTGGCATTACTGGCACC‑TTTTGGGCTT‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGTTACTTTAAAATTATCCATCTGG‑AAGG‑ATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGC‑TTTGCCGCAGCGCTGTGGGCcaa           <  2:873783/227‑1 (MQ=255)                                                                                                                                                                                                                         gcggggGGTTTTTTCTCTTCTTACTT‑TTTTGGCATTCCTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGTTACTTTAAAATTATCCATCTGG‑AAGG‑ATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGC‑TTTGCCGCAGCGCTGTGGGCCAACAATGTCaa  <  2:759803/239‑1 (MQ=255)                                                                                                                                                                                                                                     ctctctACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTgg                                                                                                                                                                                                     <  2:355813/39‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            ATTAAAGCCCTGATGTCGCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTA‑TTTTGGCATTACTGGCACC‑TTTTGGGCTG‑TCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTC‑GGCGTCCATG‑CTGAAGAGT‑GGGGTTACTTTAAAATTATCCATCTGG‑AAGG‑ATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGC‑TTTGCCGCAGCGCTGTGGGCCAACAATGTCAA  >  NC_000913/2008127‑2008581

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