Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,216,083	A→G	5.3%	L208P (CTG→CCG)	PP_4596 ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,216,083	0	A	G	5.3%	222.0 / 2.1	75	L208P (CTG→CCG)	PP_4596	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (33/38);  new base G (2/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GCGCTCAAGGCCACGCTCTTTCAGCGGGTTGAAGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCTTCGTT  >  NC_002947/5215941‑5216231                                                                                                                                                |                                                                                                                                                      gcgcTCAAGGCCACGCTCTTTCAGCGGGTTGAAGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGt                                                                                                                                                 <  4:44233/149‑1 (MQ=255)  cgcTCAAGGCCACGCTCTTTCAGCGGGTTGAAGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTg                                                                                                                                                <  2:20744/149‑1 (MQ=255)   gcTCAAGGCCACGCTCTTTCAGCGGGTTGAAGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTgg                                                                                                                                               <  8:325611/149‑1 (MQ=255)     tCAAGGCCACGCTCTTTCAGCGGGTTGAAGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCGGATCGTGGaa                                                                                                                                             >  4:373278/1‑149 (MQ=255)     tCAAGGCCACGCTCTTTCAGCGGGTTGAAGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCGGATCGTGGaa                                                                                                                                             >  4:384098/1‑149 (MQ=255)       aaGGCCACGCTCTTTCAGCGGGTTGAAGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAaca                                                                                                                                           >  2:308797/1‑149 (MQ=255)        aGGCCACGCTCTTTCAGCGGGTTGAAGCACACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGGTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAacac                                                                                                                                          <  6:387711/149‑1 (MQ=255)         ggCCACGCTCTTTCAGCGGGTTGAAGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAacac                                                                                                                                          >  7:223168/1‑148 (MQ=255)          gCCACGCTCTTTCAGCGGGTTGAAGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAAcacca                                                                                                                                        >  3:353249/1‑149 (MQ=255)            cACGCTCTTTCAGCGGGTTGAAGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACg                                                                                                                                      >  4:176282/1‑149 (MQ=255)                     tCAGCGGGTTGAAGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGccc                                                                                                                                <  1:6600/146‑1 (MQ=255)                       aGCGGGTTGAAGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGATCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGcc                                                                                                                           <  1:320921/149‑1 (MQ=255)                            gTTGAAGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCAccc                                                                                                                      >  6:158274/1‑149 (MQ=255)                                aaGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCGGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGa                                                                                                                   <  3:373278/148‑1 (MQ=255)                                aaGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCGGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGa                                                                                                                   <  3:384098/148‑1 (MQ=255)                                aaGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑g                                                                                                                 >  1:303629/1‑149 (MQ=255)                                aaGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑g                                                                                                                 >  7:222518/1‑149 (MQ=255)                                 aGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑gg                                                                                                                >  5:355985/1‑149 (MQ=255)                                  gCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑ggg                                                                                                               >  6:381372/1‑149 (MQ=255)                                    aGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGtt                                                                                                             >  7:374128/1‑149 (MQ=255)                                      acgacCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTc                                                                                                           <  5:381372/149‑1 (MQ=255)                                           cTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACAACGGTGCCCTTGCCCACCCAGA‑GGGATTCCGACa                                                                                                      >  6:271670/1‑149 (MQ=255)                                            tGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACAt                                                                                                     >  7:255362/1‑149 (MQ=255)                                                 gcCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCTCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCcg                                                                                                >  6:220538/1‑149 (MQ=255)                                                   cACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCcgcg                                                                                              >  5:147681/1‑149 (MQ=255)                                                    aCGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCgcgc                                                                                             <  2:149083/149‑1 (MQ=255)                                                    aCGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCgcgc                                                                                             >  2:19237/1‑149 (MQ=255)                                                     cGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGcgcc                                                                                            <  3:377752/149‑1 (MQ=255)                                                      gCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGcgccg                                                                                           <  4:367515/149‑1 (MQ=255)                                                      gCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGCCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGcgccg                                                                                           >  8:293426/1‑149 (MQ=255)                                                          gACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTgg                                                                                       >  3:241884/1‑149 (MQ=255)                                                          gACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTgg                                                                                       >  7:3742/1‑149 (MQ=255)                                                          gACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTgg                                                                                       >  1:357518/1‑149 (MQ=255)                                                          gACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTg                                                                                        <  2:171462/148‑1 (MQ=255)                                                                ttCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGt                                                                                 >  7:318641/1‑149 (MQ=255)                                                                   gcgcAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTgg                                                                               <  8:3742/148‑1 (MQ=255)                                                                   gcgcAGGGGTTCGAGCATGCGCGGGTGGGTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTgg                                                                               <  8:222518/148‑1 (MQ=255)                                                                    cgcAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGca                                                                             >  4:3220/1‑149 (MQ=255)                                                                      cAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGcaca                                                                           >  2:190944/1‑149 (MQ=255)                                                                       aGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACAt                                                                          >  8:139430/1‑149 (MQ=255)                                                                          ggTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGATGAGTTTCCGACATGCCCCCCCCGATGG‑ATCGTTGGACATGtt                                                                       >  1:328702/1‑149 (MQ=255)                                                                           gTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTc                                                                      <  8:318641/149‑1 (MQ=255)                                                                           gTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCAGGCTTCGTGGCACATGTTc                                                                      <  3:209070/149‑1 (MQ=255)                                                                            ttCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTCTTGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCAACCCCGA‑GGGTTTCCGACATGCCCCCGCCGCTGGCTTCGTGGCACATGTTc                                                                      >  1:57927/1‑149 (MQ=255)                                                                              cGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGa                                                                    <  7:23508/148‑1 (MQ=255)                                                                              cGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGa                                                                    <  2:357518/148‑1 (MQ=255)                                                                               gAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAg                                                                   >  6:136130/1‑148 (MQ=255)                                                                               gAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAg                                                                   <  3:176282/148‑1 (MQ=255)                                                                               gAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGGCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAg                                                                   <  8:364462/148‑1 (MQ=255)                                                                                      cgcgGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCgg                                                           >  2:394210/1‑149 (MQ=255)                                                                                      cgcgGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCACATCGTGGAACACCACGGTGCCCTTGCCTACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCgg                                                           <  1:308797/149‑1 (MQ=255)                                                                                         gggTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGACTATCGCGTCGGCCCGCTCCAGATCGTGGACCACCCCGGTGCCCTTGCCCCCCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGGACATGTTGGAGGACTGCGGGat                                                        >  4:294928/1‑148 (MQ=255)                                                                                         gggTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGaa                                                        <  4:191328/149‑1 (MQ=255)                                                                                           gtggtTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAgt                                                      >  2:363205/1‑149 (MQ=255)                                                                                            tggttggtGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAgtt                                                     <  5:158274/149‑1 (MQ=255)                                                                                                 ggtGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGttgttg                                                 <  7:293426/148‑1 (MQ=255)                                                                                                 ggtGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGttgttg                                                 <  5:220538/148‑1 (MQ=255)                                                                                                    gCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGc                                             <  3:3220/149‑1 (MQ=255)                                                                                                     ccGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGcc                                            <  2:303629/149‑1 (MQ=255)                                                                                                      cGGGGGTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCg                                           <  2:57927/149‑1 (MQ=255)                                                                                                       ggggTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTgg                                                                               >  8:303284/1‑112 (MQ=255)                                                                                                       ggggTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTgg                                                                               <  7:303284/112‑1 (MQ=255)                                                                                                       ggggTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGt                                          <  8:255362/149‑1 (MQ=255)                                                                                                          gTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAgg                                       >  6:326177/1‑149 (MQ=255)                                                                                                                 cccGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGgcacgc                                  <  1:19237/148‑1 (MQ=255)                                                                                                                  ccGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGgcacgcac                                <  6:355985/149‑1 (MQ=255)                                                                                                                   cGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGgcacgcac                                >  7:71828/1‑148 (MQ=255)                                                                                                                   cGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGgcacgcac                                <  8:71828/148‑1 (MQ=255)                                                                                                                     aTGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGaa                             <  7:139430/149‑1 (MQ=255)                                                                                                                        aCGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACa                           <  6:147681/148‑1 (MQ=255)                                                                                                                             tatCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGt                     <  1:363205/149‑1 (MQ=255)                                                                                                                                 gcgTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTaaa                  <  5:156740/148‑1 (MQ=255)                                                                                                                                 gcgTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAg                 <  4:241884/149‑1 (MQ=255)                                                                                                                                         cAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAg                 >  7:44382/1‑141 (MQ=255)                                                                                                                                         cAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAg                 <  8:44382/141‑1 (MQ=255)                                                                                                                                                aGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCTTCGtt  <  1:394210/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                      GCGCTCAAGGCCACGCTCTTTCAGCGGGTTGAAGCAGACGACCTGGGCGCCACGCTTGACCGCTTCGCGCAGGGGTTCGAGCATGCGCGGGTGGTTGGTGCCGGGGTTC‑TGCCCGATGACGAATATCGCGTCGGCCAGCTCCAGATCGTGGAACACCACGGTGCCCTTGCCCACCCCGA‑GGGTTTCCGACATGCCCGCGCCGCTGGCTTCGTGGCACATGTTCGAGCAATCCGGGAAGTTGTTGGTGCCGTAGGCACGCACGAACAGCTGGTAAAGGAACGCCGCTTCGTT  >  NC_002947/5215941‑5216231

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