Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,969,554	T→G	7.0%	L72R (CTG→CGG)	PP_2598 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,969,554	0	T	G	7.0%	222.6 / 7.5	71	L72R (CTG→CGG)	PP_2598	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/31);  new base G (3/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.01e-01

GTGGCGCTGGCCGTCGTGCTGGGCATGCTGCTGCTCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGCAGGCCCAGGCCAAACCTGACTATCGCTACGATC  >  NC_002947/2969406‑2969701                                                                                                                                                      |                                                                                                                                                      gtggCGCTGGCCGTCGTGCTGGGCATGCTGCTGCTCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCt                                                                                                                                                       <  1:147060/149‑1 (MQ=255)  tggCGCTGGCCGTCGTGCTGGGCATGCTGCTGCTCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCt                                                                                                                                                       >  6:293681/1‑148 (MQ=255)   ggCGCTGGCCGTCGTGCTGGGCATGCTGCTGCTCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGa                                                                                                                                                     >  1:23790/1‑149 (MQ=255)    gcgcTGGCCGTCGTGCTGGGCATGCTGCTGCTCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGa                                                                                                                                                     <  5:214468/148‑1 (MQ=255)    gcgcTGGCCGTCGTGCTGGGCATGCTGCTGCTCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGa                                                                                                                                                     <  5:149661/148‑1 (MQ=255)       cTGGCCGTCGTGCTGGGCATGCTGCTGCTCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGa                                                                                                                                                 <  3:240804/149‑1 (MQ=255)            cgtcgtGCTGGGCATGCTGCTGCTCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAACCCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCgg                                                                                                                                            >  7:228279/1‑149 (MQ=255)                  gCTGGGCATGCTGCTGCTCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCaaa                                                                                                                                      <  8:232819/149‑1 (MQ=255)                    tGGGCATGCTGCTGCTCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATTGCGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAag                                                                                                                                     >  1:102000/1‑149 (MQ=255)                      ggCATGCTGCTGCTCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAagca                                                                                                                                   <  2:173997/148‑1 (MQ=255)                               ctgctCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCct                                                                                                                         >  4:340109/1‑149 (MQ=255)                               ctgctCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCct                                                                                                                         >  3:339873/1‑149 (MQ=255)                                 gctCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCctgc                                                                                                                       >  6:249941/1‑149 (MQ=255)                                 gctCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCGGCTCGATCAGgcgcctg                                                                                                                                      <  5:109304/134‑4 (MQ=255)                                 gctCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCGGCTCGATCAGgcgcctg                                                                                                                                      >  6:109304/1‑131 (MQ=255)                                 gctCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCGGCTCGATCAGgcgcctg                                                                                                                                      >  6:103942/1‑131 (MQ=255)                                 gctCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCGGCTCGATCAGgcgcctg                                                                                                                                      <  5:103942/134‑4 (MQ=255)                                  ctCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCctgc                                                                                                                       <  8:169471/148‑1 (MQ=255)                                   tCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCctgctg                                                                                                                     >  3:366610/1‑149 (MQ=255)                                       cgccgcCCTGGCCTGGCCGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTc                                                                                                                 <  4:96904/149‑1 (MQ=255)                                        gccgccCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGCCGCGCGGGCCCGGCGCCTGATCGCGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTca                                                                                                                >  3:65698/1‑149 (MQ=255)                                          cgccCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTcaca                                                                                                              <  4:366610/149‑1 (MQ=255)                                           gccCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCacaa                                                                                                             >  2:196509/1‑149 (MQ=255)                                            cccTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCacaa                                                                                                             >  3:222362/1‑148 (MQ=255)                                            cccTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCacaa                                                                                                             >  8:20532/1‑148 (MQ=255)                                              ctggcctggcAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAacac                                                                                                          <  6:202988/149‑1 (MQ=255)                                              ctggcctggcAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAacac                                                                                                          >  2:180726/1‑149 (MQ=255)                                                 gcctggcAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCGGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACacgac                                                                                                       >  6:235907/1‑149 (MQ=255)                                                 gcctggcAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACacgac                                                                                                       <  8:332650/149‑1 (MQ=255)                                                 gcctggcAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGACGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACacga                                                                                                        <  8:365079/148‑1 (MQ=255)                                                   ctggcAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCacaac                                                                                                            >  5:220003/1‑142 (MQ=255)                                                   ctggcAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCacaac                                                                                                            <  6:220003/142‑1 (MQ=255)                                                      gcAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTg                                                                                                   <  5:280067/148‑1 (MQ=255)                                                               ccGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑gg                                                                                       <  4:222362/149‑1 (MQ=255)                                                               ccGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑gg                                                                                       <  1:196509/149‑1 (MQ=255)                                                               ccGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑gg                                                                                       <  3:340109/149‑1 (MQ=255)                                                                  cgGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGca                                                                                     >  3:193319/1‑148 (MQ=255)                                                                       tt‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAgcgc                                                                               <  8:228279/149‑1 (MQ=255)                                                                           gCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCt                                                                             <  1:180726/148‑1 (MQ=255)                                                                           gCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTg                                                                            >  2:236098/1‑149 (MQ=255)                                                                             cAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGa                                                                          <  5:293681/149‑1 (MQ=255)                                                                              aaaCCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGAt                                                                         >  3:7699/1‑149 (MQ=255)                                                                                  ccTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGc                                                                     >  1:30567/1‑149 (MQ=255)                                                                                    tCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTg                                                                   >  4:13529/1‑149 (MQ=255)                                                                                    tCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTg                                                                   >  5:123019/1‑149 (MQ=255)                                                                                      gTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGaagatcgg                                                                   <  4:193319/147‑8 (MQ=255)                                                                                                   gCGATGCCGCCCAGCGGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCg                                                    <  5:235907/149‑1 (MQ=255)                                                                                                   gCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCtttt                                                           >  1:153727/1‑142 (MQ=255)                                                                                                   gCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATc                                                     <  2:30567/148‑1 (MQ=255)                                                                                                   gCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCACTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCtttt                                                           >  2:88430/1‑142 (MQ=255)                                                                                                     gATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGc                                                                                                                           >  1:324493/1‑78 (MQ=255)                                                                                                     gATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGc                                                                                                                           <  2:324493/78‑1 (MQ=255)                                                                                                         ccgccCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGgcgc                                               <  2:23790/148‑1 (MQ=255)                                                                                                         ccgccCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCt                                              >  4:70141/1‑149 (MQ=255)                                                                                                          cgccCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCt                                              >  3:269638/1‑148 (MQ=255)                                                                                                           gccCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTg                                             >  6:111295/1‑148 (MQ=255)                                                                                                                 cagGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACttgt                                      <  4:269638/149‑1 (MQ=255)                                                                                                                  agGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACttgtt                                     <  2:102000/149‑1 (MQ=255)                                                                                                                    gTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACttgttg                                    >  3:189232/1‑148 (MQ=255)                                                                                                                        ccAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGcaggc                               >  7:213920/1‑149 (MQ=255)                                                                                                                         cAAGACGATCGACGAGAGGGCCAGGCTCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGTCTCATAACACGACGGTGCACACCCACTGGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGcagg                                <  1:106809/149‑1 (MQ=255)                                                                                                                           acaCGATCGACCAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGATGGCTCACAACACGACGGTGAACACCCA‑‑GGCACAGCGCCTGGATCGGATGGCGCTTTTAAGATCGGGAGCTGAACTTGTTGcatgccca                            >  8:315813/3‑149 (MQ=255)                                                                                                                           aGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGcaggcc                              >  4:182392/1‑147 (MQ=255)                                                                                                                            gACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGcaggcccag                           <  3:182392/149‑1 (MQ=255)                                                                                                                                 tCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACttgt                                      <  4:125518/133‑1 (MQ=255)                                                                                                                                 tCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACttgt                                      >  3:125518/1‑133 (MQ=255)                                                                                                                                    acgaGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGCAGGCCCAGGCCAAAcc                   >  5:132835/1‑149 (MQ=255)                                                                                                                                           ggCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGCAGGCCCAGGCCAAACCTGACTAt            >  3:294354/1‑149 (MQ=255)                                                                                                                                            gctttgcgcCGGCTCGATCAGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGCAGGCCCAGGCCATACCTGACTATc           >  6:183832/6‑149 (MQ=255)                                                                                                                                            gCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGCAGGCCCAGGCCAAACCTGACTATc           <  1:154843/149‑1 (MQ=255)                                                                                                                                              cAGGCGCCTTATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGGGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGCAGGCCCAGGCCAAACCTGACTATCGc         <  8:261462/149‑1 (MQ=255)                                                                                                                                              cAGCCGCCTGATCGCGCCGTCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGCTCGGCTGGCGCTTTTAAGATCGGGCGCTGACCTTGTTCCAGGCCCAGGCCAATCCTGACTATCGc         >  8:68878/1‑149 (MQ=255)                                                                                                                                                 gcgcCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGCAGGCCCAGGCCAAACCTGACTATCGCTAc      >  7:17269/1‑149 (MQ=255)                                                                                                                                                      tGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGCAGGCCCAGGCCAAACCTGACTATCGCTACGATc  >  5:191365/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GTGGCGCTGGCCGTCGTGCTGGGCATGCTGCTGCTCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGCAGGCCCAGGCCAAACCTGACTATCGCTACGATC  >  NC_002947/2969406‑2969701

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