Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,969,563	C→A	6.4%	P75Q (CCG→CAG)	PP_2598 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,969,563	0	C	A	6.4%	234.7 / 2.3	72	P75Q (CCG→CAG)	PP_2598	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (36/31);  new base A (3/2);  total (39/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 = 1.00e+00

CGTCGTGCTGGGCATGCTGCTGCTCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGCAGGCCCAGGCCAAACCTGACTATCGCTACGATCCACGCCTGC  >  NC_002947/2969417‑2969710                                                                                                                                                    |                                                                                                                                                      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                                                                                                                                               >  6:103942/1‑131 (MQ=255)                      gctCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCGGCTCGATCAGgcgcctg                                                                                                                                               >  6:109304/1‑131 (MQ=255)                      gctCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCGGCTCGATCAGgcgcctg                                                                                                                                               <  5:109304/134‑4 (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                                                                                                                      >  8:20532/1‑148 (MQ=255)                                 cccTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCacaa                                                                                                                      >  3:222362/1‑148 (MQ=255)                                   ctggcctggcAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAacac                                                                                                                   >  2:180726/1‑149 (MQ=255)                                   ctggcctggcAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAacac                                                                                                                   <  6:202988/149‑1 (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                                                                                                                     <  6:220003/142‑1 (MQ=255)                                        ctggcAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCacaac                                                                                                                     >  5:220003/1‑142 (MQ=255)                                           gcAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTg                                                                                                            <  5:280067/148‑1 (MQ=255)                                                    ccGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑gg                                                                                                <  1:196509/149‑1 (MQ=255)                                                    ccGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑gg                                                                                                <  3:340109/149‑1 (MQ=255)                                                    ccGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑gg                                                                                                <  4:222362/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                                                                            >  5:123019/1‑149 (MQ=255)                                                                         tCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTg                                                                            >  4:13529/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)                                                                                                                                             aTCGAGCCGGCGCAAAGCAGCACTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGCAGGCCCAGGCCAAACCTGACTATCGCTACGATCCAc        <  1:88430/149‑1 (MQ=255)                                                                                                                                               cGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGCAGGCCCAGGCCAAACCTGACTATCGCTACGa             >  2:203559/1‑142 (MQ=255)                                                                                                                                               cGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGCAGGCCCAGGCCAAACCTGACTATCGCTACGa             <  1:203559/142‑1 (MQ=255)                                                                                                                                                gAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGCAGGCCCAGGCCAAACCTGACTATCGCTACGATCCACGc      <  3:13529/148‑1 (MQ=255)                                                                                                                                                   ccGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGCAGGCCCAGGCCAAACCTGACTATCGCTACGATCCACGCCTgc  >  7:343033/1‑149 (MQ=255)                                                                                                                                                   ccGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGCAGGCCCAGGCCAAACCTGACTATCGCTACGATCCACGCCTgc  >  1:82341/1‑149 (MQ=255)                                                                                                                                                   ccGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGCAGGCCCAGGCCAAACCTGACTATCGCTACGATCCACGCCTgc  <  4:339873/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CGTCGTGCTGGGCATGCTGCTGCTCGGCGCCGCCCTGGCCTGGCAGGGTTACCGCGGAATT‑CGCCAAACCCTCGTGGCCGCCGCTGGCGATGCCGCCCAGCAGGTTGCCAAGACGATCGACGAGCGGGCCAGGCGCCTGATCGAGCCGGCGCAAAGCAGCATTCGCCTGCTGGCTCACAACACGACGGTGCACACCCA‑‑GGCACAGCGCCTGGATCGGCTGGCGCTTTTAAGATCGGGCGCTGAACTTGTTGCAGGCCCAGGCCAAACCTGACTATCGCTACGATCCACGCCTGC  >  NC_002947/2969417‑2969710

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