Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,893,842	A→G	7.9%	F406S (TTC→TCC) ‡	yegH ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,893,842	0	A	G	7.9%	171.0 / 3.0	62	F406S (TTC→TCC) ‡	yegH	putative transporter
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base G (2/3);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTGGCTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGC  >  NC_002947/3893695‑3893988                                                                                                                                                     |                                                                                                                                                   aCCTGGCTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGaa                                                                                                                                                   <  6:40468/149‑1 (MQ=255)  ccTGGCTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAAc                                                                                                                                                  <  4:71712/149‑1 (MQ=255)      gCTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTc                                                                                                                                                <  7:230387/147‑1 (MQ=255)       cTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGtt                                                                                                                                             <  4:234387/149‑1 (MQ=255)          aGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGa                                                                                                                                           <  1:4467/148‑1 (MQ=255)           ggTTCAGCGCACCGCTGCCCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGAc                                                                                                                                          >  4:50154/1‑148 (MQ=255)             ttCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACGAGTTCGGTGACTTCACCGGGCT‑GTTGACCat                                                                                                                                       <  6:248545/149‑2 (MQ=37)             ttCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGAcccc                                                                                                                                       >  1:196068/1‑147 (MQ=255)                 gcgcACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAg                                                                                                                                   >  5:148927/1‑149 (MQ=255)                 gcgcACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAg                                                                                                                                   <  2:97609/149‑1 (MQ=255)                      ccGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAat                                                                                                                              <  5:174338/149‑1 (MQ=255)                       cGCTGTCCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGTCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAata                                                                                                                             <  1:278875/149‑1 (MQ=255)                       cGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAata                                                                                                                             >  2:125459/1‑149 (MQ=255)                              cacaAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTc                                                                                                                       >  3:97941/1‑148 (MQ=255)                               acaAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGa                                                                                                                     >  3:46548/1‑149 (MQ=255)                                       ttCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTc                                                                                                                       <  1:129148/139‑1 (MQ=255)                                       ttCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTc                                                                                                                       >  2:129148/1‑139 (MQ=255)                                         cccccTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACcgcg                                                                                                            >  8:139618/1‑148 (MQ=255)                                         cccccTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACcgc                                                                                                             <  1:125459/147‑1 (MQ=255)                                          ccccTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCgcgc                                                                                                           <  4:24557/148‑1 (MQ=255)                                          ccccTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCa                                                                                                          >  3:17313/1‑149 (MQ=255)                                             ctcctGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCAttt                                                                                                       <  5:11266/149‑1 (MQ=255)                                                ctGAACGATGCCCGTGCCTTCAACCTGACTGTCATCCGGCAGCTCACCGTCAATCCACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTAACCGAACTCGTTGAGCACTAAGGCAATATGGGTCGACTGACCGCGCATTTGc                                                                                                     <  7:240975/148‑1 (MQ=255)                                                ctGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGc                                                                                                     <  7:139618/148‑1 (MQ=255)                                                   aaCGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTc                                                                                                   <  4:46548/147‑1 (MQ=255)                                                           ccGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTg                                                                                                                                            <  2:153873/98‑1 (MQ=255)                                                           ccGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTg                                                                                                                                            >  1:153873/1‑98 (MQ=255)                                                            cGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGa                                                                                        <  6:148927/149‑1 (MQ=255)                                                                  ttCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCCGGGCATTGAGGATGc                                                                                  >  1:230579/1‑149 (MQ=255)                                                                   tCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTcc                                                                                                  >  4:91496/1‑132 (MQ=255)                                                                   tCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTcc                                                                                                  <  3:91496/132‑1 (MQ=255)                                                                       ccTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACAACGAAGGCAATATGGGTAGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAg                                                                             >  8:171507/1‑149 (MQ=255)                                                                         tCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATAGAGGATGCTGAAGct                                                                           >  1:14115/1‑149 (MQ=255)                                                                           aCTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGctct                                                                         >  4:72927/1‑149 (MQ=255)                                                                                               cGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCACGTTCAACCGAGcc                                                     >  4:303387/1‑149 (MQ=255)                                                                                                     tCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGc                                                >  7:202623/1‑148 (MQ=255)                                                                                                       gACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCAt                                             <  8:202623/149‑1 (MQ=255)                                                                                                        tcTCAAGGATGGCGGTCTTGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCAt                                             <  2:232411/147‑1 (MQ=255)                                                                                                          tCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCt                                          >  1:61058/1‑149 (MQ=255)                                                                                                            aaGGATGTCGGTCATGGTCACC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGcgaa                                        >  1:232411/1‑146 (MQ=255)                                                                                                            aaGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCttt                                        >  5:300284/1‑149 (MQ=255)                                                                                                             aGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTc                                       <  4:97941/149‑1 (MQ=255)                                                                                                              ggATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTc                                       >  1:93312/1‑148 (MQ=255)                                                                                                               gATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGGCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCa                                     <  5:320619/149‑1 (MQ=255)                                                                                                                 tGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAgg                                   <  2:61058/149‑1 (MQ=255)                                                                                                                  gTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTc                                                                                                   >  2:9829/1‑84 (MQ=255)                                                                                                                     cttcgtGGTCAACGAGTTCGGTGACTTCACCGGGCT‑GTTGACCATGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCa                                     <  7:131350/138‑1 (MQ=25)                                                                                                                     cttcgtGGTCAACGAGTTCGGTGACTTCACCGGGCT‑GTTGACCATGAAGGCAATATGGGTCGACTGACCGCGAATTTGCTCCAGGGCATTGAGGATGCTGAAGATCTCCAGAAGGTTCAACGGCGCCCGGGCCATGATTTCCa                                     >  8:131350/6‑143 (MQ=21)                                                                                                                      ttcgtGGTCAACGAGTTCGGTGACTTCACCGGGCT‑GTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCggg                              >  5:31429/5‑149 (MQ=37)                                                                                                                      gTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCggg                              <  3:72927/149‑1 (MQ=255)                                                                                                                      gTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCggg                              >  7:255538/1‑149 (MQ=255)                                                                                                                          tGGTCAAC‑ACCCCGGTGAAGTCACCGAACTCGTTCACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctgg                          <  2:196068/149‑1 (MQ=255)                                                                                                                           ggTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctggc                         >  1:143670/1‑149 (MQ=255)                                                                                                                            gTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctggct                        >  6:300247/1‑149 (MQ=255)                                                                                                                             tCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGgctggctg                       >  7:52300/1‑149 (MQ=255)                                                                                                                                     cccGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCgaga                <  5:300247/149‑1 (MQ=255)                                                                                                                                      ccGGTGAAGTCACCGAACTAGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGagag               <  4:17313/149‑1 (MQ=255)                                                                                                                                       cGGTGACTTCACCGGGCT‑GTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagc              <  6:31429/148‑1 (MQ=255)                                                                                                                                       cGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGagag               >  3:278522/1‑148 (MQ=255)                                                                                                                                             aGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcc        <  6:300284/149‑1 (MQ=255)                                                                                                                                                 aCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGTGCTGGCTGCCCGAGagcagttcctt      >  7:50959/1‑147 (MQ=255)                                                                                                                                                  ccGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcag   <  8:255538/149‑1 (MQ=255)                                                                                                                                                    gAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATCCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcag   >  7:108936/1‑147 (MQ=255)                                                                                                                                                     aaCTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcagc  <  7:197684/147‑1 (MQ=255)                                                                                                                                                     aaCTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcagc  >  8:197684/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                   ACCTGGCTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGC  >  NC_002947/3893695‑3893988

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