Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	854,311	A→G	26.1%	E222G (GAG→GGG)	ybiR →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	854,311	0	A	G	26.1%	60.8 / 2.3	39	E222G (GAG→GGG)	ybiR	putative transporter
Reads supporting (aligned to +/- strand):  ref base A (9/18);  new base G (5/6);  total (14/25)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGC‑TCCTGCTCCTGTGCTG‑GTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGC‑TGGTGTGGAGTTGTCTGGGGCTGTATATCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGACGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAACGTGAGTCATCTATC‑TGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAAC  >  NC_000913/854089‑854534                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                 tCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGC‑TCCTGCTCCTGTGCTG‑GTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGCTTGGAAACCGCGGC‑TGGGGGGGAGTTGGCTGGGGCTGTAATCCGGCTTTCTGACGGCGCTGGGGGTCAAAAAAGAGGtgg                                                                                                                                                                                                                             >  1:71118/1‑227 (MQ=255) tCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCGGACGC‑TCCTGCTCCTGTGCTG‑GTGTTGTTTCCCTGGAAAGGCGATGCAATTCCATCCGGGGGTGCAAACACCGGAGTGGAAACCGCGGC‑TGGGGTGGGGTTGACTGTGGCTGTATAACGGCTTTCCGACGGCGCGGGAGTTAACACAAGAGgtgt                                                                                                                                                                                                                             >  1:500022/1‑228 (MQ=255) tCTGGGGACGTTCTGGTCTTTCGGTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGC‑TCCTGCTCCTGTGCTG‑GTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGC‑TGGTGTGGAGTTGTCTGGGGCTGTATATCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGAGCtgt                                                                                                                                                                                                                             <  2:16371/228‑1 (MQ=255)   tGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGC‑TCCTGCTCCTGTGCTG‑GTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGC‑TGGGGGTGAGTTGTCTGGGGCTGGATTCTGCCTTTCTGACGGCCGTGGAGGTCAAACAAAGGCTgggg                                                                                                                                                                                                                           >  1:492527/1‑228 (MQ=255)     gggACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTTACGC‑TCCTGCTCCTGTGCTG‑GTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGC‑TGGTGTGGAGTTGTCTGGGGCTGTATATCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGAGCTGTgggg                                                                                                                                                                                                                         >  2:529339/1‑228 (MQ=255)           ttCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGC‑TCCTGCTCCTGTGCTG‑GTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGCGTGGAAACCGCGGC‑GGCGATGGAGTTGTATGGGGCTGGATATCGTATTTCTGACGGCGCCGGAGTTCGAAAAAGAGCTGTGGGGAATggg                                                                                                                                                                                                                   >  1:144558/1‑227 (MQ=255)           ttCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGC‑TCCTGCTCCTGTGCTG‑GTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGGGCAAACACCGCAGTGGAAACCGCGGC‑TGGGGTGGAGTTGTCTGGGGCTGCCTCACTCCTTTCTGACGCCGCTGGGGGTCAAAGAAGGGCGGTggggctgggg                                                                                                                                                                                                                   >  1:2503/1‑222 (MQ=255)               ggTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGC‑TCCTGCTCCTGTGCTG‑GTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGC‑TGGTGTGGAGTTGTCTGGGGCTGTATATCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGGGCTGTGGGGACTGGTGCtt                                                                                                                                                                                                               >  2:286197/1‑228 (MQ=255)                                          cAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGC‑TCCTGCTCCTGTGCTG‑GTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGC‑TGGTGCGGGGTTGTCTGGGGCTGTATATCGTCTTTCTGACGGCGCTGGAGGTCAAACAAGGGCGTGGGGGACTGTGGATTGTggcggggtgggtttggctgctgggc                                                                                                                                                                                    >  1:415614/1‑208 (MQ=255)                                                 gtccgtggggtgccgcAAGGATGCTTACGC‑TGTGGGTCCTGTGCGG‑GGGTGGTTTCCCGGGAAAGGCGATGCAATACGATACGGGGGTGCAAACACAGGAGGGGAAACCGCGGC‑TGGTGTGGAGTTGTCTGGGGATGTATATCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTcgcg                                                                                                                                                                             <  1:529339/215‑1 (MQ=255)                                                 cACCGCTGGCTGGCGCAATGATTCTGACGC‑TCCTGCTCCTGGGCTG‑GTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGC‑TGGTGTGGAGTTGTCTGGGGCTGTATATCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTcgcg                                                                                                                                                                             <  2:492527/228‑1 (MQ=255)                                                     gttggttgggcgatgggtggtGACCC‑TTCTGGCCCTGTGCGG‑GTGTTGGTTCCCTGGAAAGGCGATGCAATGCCATACGGGGGTGCAAACACCGGAGTGGAATGCGCGGC‑TGGTGTGGAGTTGTCTGGGGCTGTATATCGTCTTTCTGGCGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCgtggt                                                                                                                                                                         <  1:527936/208‑1 (MQ=255)                                                       tggttgggGCAAGGATGGGGGCGC‑GCCTGCTCCGGGGCGT‑GGGTTGTTTCCCTGGAAAGGTGATGGAATACCGGAGGGGGGTGCAAACACCGGAGTGGAAACCGCGGC‑TGGTGTGGAGTGGGCTGGGGCTGTAGAGCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGc                                                                                                                                                                       <  1:222053/228‑1 (MQ=255)                                                       tggctggcGCAATGACGCTGACGC‑TCCTGCTTCTGTGCTG‑GTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGC‑TGGTGTGGAGTTGTCTGGGGCTGTATATCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGc                                                                                                                                                                       <  2:500022/228‑1 (MQ=255)                                                        gggtgtggcAATGATGGCGACGCGTCCT‑CCCCTGTTCTGTGTGTT‑TGTCACTGGAAAGGCGATGCAATGCCATACGGGGGTGCAAACCCCTGAGTGGAAACCGCGGC‑TGGTGTGGGAGTGTCTGGGGCTGTATATCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGTTTTGCGCTGCTGGCACGTCGCGTGGTGCt                                                                                                                                                                      <  1:35296/221‑1 (MQ=255)                                                         gctggcGCAATGATGCTGACGC‑TCCTGCTCCTGTGCTG‑GTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGC‑TGGTGTGGAGTTGTCTGGGGCTGTATATCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGGGCTGTGGGGACTGGTGATTGTGCCGGCAGGCTTTGCGCTGCTGGCACgcgccgcgggggtc                                                                                                                                                                     >  2:547178/1‑215 (MQ=255)                                                                  atgatgCTGACGC‑TCCTGGTCCTGTGCTG‑GTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACTCCGGAGTGGAAACCGCGGC‑TGGTGTGGAGTTGTCTGGGGCTGTATATCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCAGTGtgga                                                                                                                                                             <  2:71118/227‑1 (MQ=255)                                                                                 cccgctcctgTGCTG‑GTGTAGTTTCCCTGGAAAGTCGATGCAATTACATACGGGGGTGCAAACACCGGAGGGGAAACCGCGGC‑TGGTGTGGAGTTGTCTGGGGCTGTATATCGTCTTTCTGACGGCGATGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGACgctgctgct                                                                                                                                              <  2:495968/225‑1 (MQ=255)                                                                                        ctgTGCTG‑GTGTTGTTTCCCTGGAAAGGCGCTGCAATACCATACGTGGGTGCAAACACCGGAGTGGAAACCGCGGC‑CGGAGTGGAGTTGTCTGGGGCTGTATATCGTCTTTCTGCCGGCGCTGGGGTTCAAACAAGAGCTGTGGGGACTGGTACTTTTGGCGGCAGCCTTTGCCCTCCGTGCCCGCTGGGTGGGGTTCGGTGTGGACTGGTCGCTGCTGCCggggggt                                                                                                                                       >  1:446147/1‑223 (MQ=255)                                                                                              tG‑GTTTTGTTTTCCTTGAAAGGGGGTGGATGTGCTTACGGGGGTGCAATCGCCGGGGTGGAAGCGCGGGG‑TGGTGGGGGGTTGTCTGGGGCTGTGTGTCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGGGCTGTGGGGACTGGTGTTTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGGCGCTGCTGCTGGTGTTTATGGCg                                                                                                                                 <  1:314656/228‑1 (MQ=255)                                                                                                           ccTGGAAAGGCGATGCGATGCCTGACGGGGGTGCAAACACCGGAGTGGAAACCGCGGC‑TGGTGTGGAGTTGTCTGGGGCTGTATAGCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGACGCTGCTGCTGGTGTTTATGGCGATGTTTATCGa                                                                                                                      <  2:219255/227‑1 (MQ=255)                                                                                                           ccTGGAAAGGCGAGTCAATACCATGCGGGGGTGCAAACACCGGAGTGGAAACCGCGGC‑TGGTGTGGAGTTGTCTGGGGCTGTATATCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGACGCTGCTGCTGGTGTTTATGGCGATGTTTATCGa                                                                                                                      <  2:415614/227‑1 (MQ=255)                                                                                                            cTGGGAAGGCGATGGAATTCCATAGGGGGGTGCAAACACCGGAGTGGAAACCGCGGT‑TGGTGTGGAGTTGTCTGGGGCTGTATATCGTCTTTCTGACGGCGCGGGAGTTCAAACAAGGGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGACGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACg                                                                                                                    <  1:667758/228‑1 (MQ=255)                                                                                                                       ttggaatgccttgcGGGGGTGCAAACCCCGGGGTGGAAACCGGGGC‑TGGTGTGGGGTTGTCTGGGGCTTTATATCGTCTTTCGGACGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGACGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTg                                                                                                         <  1:449292/215‑1 (MQ=255)                                                                                                                                  tgggggggTGCAAATGCGGGGGGGGAAACCGGGGC‑TGGTGGGGAGTTGTCTGGGGCTGTATATCGTCTTTCTGGCGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGACGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTc                                                                                               <  1:286197/224‑1 (MQ=255)                                                                                                                                  tACGGGGGTGCAAACACCGGAGTGGAAACCGCGGC‑TGGTGTGGAGTTGTCTGGGGCTGTATATCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGTCCGGCTTTGCTCTGCTGGCCCGCGGCGGGGGGCTCGATGTGGCCCGGAGGCTGCCGCTGGGGGTTTTGGGGCTGGTTTTGTTCGTCCACTTTCTGAgccccctccc                                                                                              >  1:537738/1‑221 (MQ=255)                                                                                                                                          tGCAAACACCGGAGTGGAAACCGCGGC‑TGGTGTGGAGTTGTCTGGGGCTGTATATCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGACGCTGCTGCTGGGGTTTATGGCGATGTTTATTGAGGGCCATTTACTGACCCAGCTTTCAGCGGTGc                                                                                      >  2:150298/1‑228 (MQ=255)                                                                                                                                               acacCGGAGTGGAAACCGGGGT‑TGGTGTGGGGTTGTCGGGGGCTGTATATCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGACGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGc                                                                                 <  2:446147/228‑1 (MQ=255)                                                                                                                                                   cGGGGTGGAAATCGCGGCGGGGTGTGG‑GTGGTCTGGGGCTTTATATCGGCTTTCGGACGGGGGTTGAGTTCAAACAAGGGGCGGGGGGGCTGGTGATTGTGGCGGCAGGGTTTGCGCTTCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGACGCTGCTGCTGGTGTTTATGGCGATATTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCgtg                                                                              <  1:328427/227‑1 (MQ=255)                                                                                                                                                                          tgggggTTGGCTGGGGATGTATAGCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGACGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAACGTGAGTCatctatc                                                       <  2:524155/223‑1 (MQ=255)                                                                                                                                                                                ttGTTTGGGGTTGTATTTCGTGTTTCTGGCGGCGGTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGACGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAACGTGAGTCATCTATC‑TGAAc                                                 <  1:150298/227‑1 (MQ=255)                                                                                                                                                                                       ggggTGTTTTTCGTTTTTCGGGCGGCGCTGGAGTTCAAACAAGGGCTGTGGGGACTGGTGATTGGGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGACGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAACGTGAGTCATCTATC‑TGAACCCGGGTTa                                         <  1:547178/228‑1 (MQ=255)                                                                                                                                                                                                    tCTTTCTGACGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCCGGGTGGACTGGACCCTGCTGCTGGTGTTTATGGCGATGTTTTTCGGCGGTCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAACGTGTGTCCTCTTTC‑TTAACCCGGGTTATGGTTAACGGCa                             >  2:223788/1‑227 (MQ=255)                                                                                                                                                                                                      tttCGGGCGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGGGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGACGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAACGTGAGTCATCTATC‑TGAACCCGGGTTATGGTTAACGGCAAt                           <  2:214888/227‑1 (MQ=255)                                                                                                                                                                                                                cgcTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCCGGGGGGGCTGGGCGGCGCTGCTGGTGGTTTTGGGGGGGTTTTTCTGCGTCCCTTTACTGACCCACCCTTCAGCGTGTGAAGGCGTGTTGGGGAAAGTGAGGACTCCATCTTGACCCCGGGTT‑TGGTTTACGGCAATTGGGTTTTCGc                >  1:536047/1‑228 (MQ=255)                                                                                                                                                                                                                    gggtgtCAAACGAGGGGTGGGGGGACTGGTGAATGTGGCGGCAGGCTTTGAGCTGCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGACGCTGCTGCTGGTGTTTATGGCGATGATTAGCGACGTCCATTTACTGACCCAGCTTCCGGCGTTGCAAGGCGTGTTGGGTAACGTGAGTCATCTATC‑TGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGt            <  2:123037/223‑1 (MQ=255)                                                                                                                                                                                                                    ggATTTCAGACAAGTGGTGTGGGTATTGGTGATTGTGGCGGCAGGCTTTGCGCGGCTTGCACGTCGGGTGGTGCTCAGTGTGGACTGGACGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAACGTGAGTCATCTATC‑TGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGt            <  2:503840/228‑1 (MQ=255)                                                                                                                                                                                                                        ttCAAACAAGAGCTGGGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGACGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAACGTGAGTCATCTATC‑TGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGAt         <  2:537738/227‑1 (MQ=255)                                                                                                                                                                                                                               aaGGGCTGTGGGGACTGGTGATTGTGGCGGCCGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGACGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAACGTGAGTCATCTATC‑TGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAAc  <  2:536047/227‑1 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                 TCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGC‑TCCTGCTCCTGTGCTG‑GTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGC‑TGGTGTGGAGTTGTCTGGGGCTGTATATCGTCTTTCTGACGGCGCTGGAGTTCAAACAAGAGCTGTGGGGACTGGTGATTGTGGCGGCAGGCTTTGCGCTGCTGGCACGTCGCGTGGTGCTCAGTGTGGACTGGACGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAACGTGAGTCATCTATC‑TGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAAC  >  NC_000913/854089‑854534

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