Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,786,763	T→G	6.3%	L59R (CTG→CGG)	trxC →	thioredoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,786,763	0	T	G	6.3%	271.2 / 5.0	87	L59R (CTG→CGG)	trxC	thioredoxin
Reads supporting (aligned to +/- strand):  ref base T (46/35);  new base G (4/2);  total (50/37)
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

CCACTGCAACGGCCTCAACCGCCTGCCCGCCGAGCGCCTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGG  >  NC_002947/5786617‑5786906                                                                                                                                                   |                                                                                                                                                ccACTGCAACGGCCTCAACCGCCTGCCCGCCGAGCGCCTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCgctgc                                                                                                                                               >  5:348804/1‑149 (MQ=255)    cTGCAACGGCCTCAACCGCCTGCCCGCCGAGCGCCTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTgg                                                                                                                                            >  5:389520/1‑149 (MQ=255)    cTGCAACGGCCTCAACCGCCTGCCCGCCGAGCGCCTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTgg                                                                                                                                            >  8:262975/1‑149 (MQ=255)    cTGCAACGGCCTCAACCGCCTGCCCGCCGAGCGACTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTgg                                                                                                                                            >  7:331533/1‑149 (MQ=255)      gCAACGGCCTCAACCGCCTGCCCGCCGAGCGCCTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTc                                                                                                                                          >  1:134711/1‑149 (MQ=255)       cAACGGCCTCAACCGCCTGCCCGCCGAGCGCCTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTc                                                                                                                                          >  4:18876/1‑148 (MQ=255)        aaCGGCCTCAACCGCCTGCCCGCCGAGCGCCTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGa                                                                                                                                        <  5:319437/149‑1 (MQ=255)            gCCTCAACCGCCTGCCCGCCGAGCGCCTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATc                                                                                                                                    >  6:382059/1‑149 (MQ=255)             ccTCAACCGCCTGCCCGCCGAGCGCCTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCt                                                                                                                                   <  6:58015/149‑1 (MQ=255)              cTCAACCGCCTGCCCGCCGAGCGCCTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTg                                                                                                                                  <  4:230328/149‑1 (MQ=255)                 aaCCGCCTGCCCGCCGAGCGCCTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGc                                                                                                                               <  7:262975/149‑1 (MQ=255)                 aaCCGCCTGCCCGCCGAGCGCCTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCGGCTGGTc                                                                                                                                          >  6:348309/1‑138 (MQ=255)                 aaCCGCCTGCCCGCCGAGCGCCTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCGGCTGGTc                                                                                                                                          <  5:348309/138‑1 (MQ=255)                 aaCCGCCTGCCCGCCGAGCGACTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTggg                                                                                                                                <  8:331533/148‑1 (MQ=255)                           ccgccgAGCGCCTGGGCAATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTCCCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGCCCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGc                                                                                                                     >  3:376973/1‑149 (MQ=255)                              ccgAGCGCCTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGc                                                                                                                  <  6:40315/149‑1 (MQ=255)                              ccgAGCGCCTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGc                                                                                                                  <  1:237287/149‑1 (MQ=255)                                gAGCGCCTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGccc                                                                                                                >  4:423737/1‑149 (MQ=255)                                     ccTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGc                                                                                                                               <  3:143606/129‑1 (MQ=255)                                     ccTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGc                                                                                                                               >  4:143606/1‑129 (MQ=255)                                      cTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCaa                                                                                                          <  6:348804/149‑1 (MQ=255)                                        gggCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCAGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGt                                                                                                        >  7:327671/1‑149 (MQ=255)                                                cgcCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGcc                                                                                                <  4:279929/149‑1 (MQ=255)                                                 gcCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGccc                                                                                               >  3:248711/1‑149 (MQ=255)                                                 gcCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGccc                                                                                               >  5:254470/1‑149 (MQ=255)                                                 gcCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGccc                                                                                               >  6:287137/1‑149 (MQ=255)                                                 gcCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGccc                                                                                               >  1:207312/1‑149 (MQ=255)                                                   cGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGccccc                                                                                             <  6:169398/149‑1 (MQ=255)                                                   cGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGccccc                                                                                             <  2:207312/149‑1 (MQ=255)                                                          cGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTc                                                                                       <  6:254470/148‑1 (MQ=255)                                                           ggTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCg                                                                                      >  1:277295/1‑148 (MQ=255)                                                           ggTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGa                                                                                     <  8:327671/149‑1 (MQ=255)                                                            gTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGaa                                                                                    >  6:444033/1‑149 (MQ=255)                                                               gCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAagag                                                                                 <  2:134711/149‑1 (MQ=255)                                                               gCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAagag                                                                                 >  5:20749/1‑149 (MQ=255)                                                                 tGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGgc                                                                               >  7:288481/1‑149 (MQ=255)                                                                  gCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAagag                                                                                 >  6:266496/1‑146 (MQ=255)                                                                  gCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAagag                                                                                 <  5:266496/146‑1 (MQ=255)                                                                    aaaaaGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGgccgc                                                                            >  6:280126/1‑149 (MQ=255)                                                                      aaaGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGccc                                                                          >  8:405443/1‑149 (MQ=255)                                                                            tgtgCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGc                                                                    >  5:1416/1‑149 (MQ=255)                                                                                cTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGc                                                                >  6:319486/1‑149 (MQ=255)                                                                                   ttGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTg                                                              <  3:18876/148‑1 (MQ=255)                                                                                    tGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGCCTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCCCCCGGAAGCTTGCTGGc                                                            >  1:402973/1‑149 (MQ=255)                                                                                      tgCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCg                                                          <  3:91351/148‑1 (MQ=255)                                                                                      aGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGcc                                                           <  6:389520/148‑1 (MQ=255)                                                                                        cacaCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCg                                                          >  4:91351/1‑147 (MQ=255)                                                                                        cACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCt                                                        >  6:389078/1‑149 (MQ=255)                                                                                           ccccTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTg                                                                 >  8:339870/1‑137 (MQ=255)                                                                                           ccccTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTg                                                                 <  7:339870/137‑1 (MQ=255)                                                                                               ttCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCt                                                 >  3:411221/1‑149 (MQ=255)                                                                                                  gAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCg                                                                                      <  8:391068/109‑1 (MQ=255)                                                                                                  gAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCg                                                                                      >  7:391068/1‑109 (MQ=255)                                                                                                  gAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTgg                                               <  5:389078/148‑1 (MQ=255)                                                                                                  gAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGc                                              >  8:39030/1‑149 (MQ=255)                                                                                                  gAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGc                                              <  2:28865/149‑1 (MQ=255)                                                                                                      ttACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAg                                          >  3:6479/1‑149 (MQ=255)                                                                                                         ccGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCt                                        <  6:20749/148‑1 (MQ=255)                                                                                                         ccGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTc                                       >  4:271870/1‑149 (MQ=255)                                                                                                          cGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTc                                       >  6:433581/1‑148 (MQ=255)                                                                                                           gAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGa                                     <  3:423737/149‑1 (MQ=255)                                                                                                            aaGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGa                                     >  7:36123/1‑148 (MQ=255)                                                                                                                 aGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCg                                >  1:21862/1‑148 (MQ=255)                                                                                                                 aGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCg                                >  1:166300/1‑148 (MQ=255)                                                                                                                 aGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCg                                >  1:289668/1‑148 (MQ=255)                                                                                                                  gCTACGCCAGGCAGCTCAAAGGCTACCTTGCGCTGGTGGTTGACATCTTGGCCGACTGGTGCGGCCCGTGCAACTCATTCGGCCCCACCTTCGAAGTGGGCGGCCGGAAGGTTGGGGGCCGGTGCCGGCTGGCCAAGCTCGACAGCGaa                              <  2:402973/149‑1 (MQ=255)                                                                                                                  gCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCg                                >  1:171774/1‑147 (MQ=255)                                                                                                                  gCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGaa                              <  7:89319/149‑1 (MQ=255)                                                                                                                    tACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGc                            >  8:24148/1‑149 (MQ=255)                                                                                                                    tACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGc                            >  2:149371/1‑149 (MQ=255)                                                                                                                    tACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGc                            >  3:256200/1‑149 (MQ=255)                                                                                                                      cGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGc                            <  4:256200/147‑1 (MQ=255)                                                                                                                      cGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCa                          <  3:271870/149‑1 (MQ=255)                                                                                                                       gccagccaGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCaa                         <  1:149371/149‑1 (MQ=255)                                                                                                                        ccagccagATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAAc                        <  5:319486/149‑1 (MQ=255)                                                                                                                           gccagATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGt                     <  3:220661/149‑1 (MQ=255)                                                                                                                              agATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAAc                  >  2:310702/1‑149 (MQ=255)                                                                                                                               gATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAAcc                 >  3:400434/1‑149 (MQ=255)                                                                                                                                aTCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAAc                  <  6:60515/147‑1 (MQ=255)                                                                                                                                aTCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAAc                  >  5:60515/1‑147 (MQ=255)                                                                                                                                    atgtcGACCAGCAGCGGCAGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCAACAGCGAAGCCAACCGTAACCTCg              >  3:242182/5‑147 (MQ=255)                                                                                                                                    aaGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCa            <  6:1416/149‑1 (MQ=255)                                                                                                                                      gtcGACCAGCAGCGGCAGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAgg          >  3:257497/3‑149 (MQ=255)                                                                                                                                       tcGACCAGCAGCGGCAGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTgg                                               >  5:244036/2‑111 (MQ=255)                                                                                                                                       tcGACCAGCAGCGGCAGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTgg                                               <  6:244036/110‑1 (MQ=255)                                                                                                                                          aCCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCa       >  7:159227/1‑148 (MQ=255)                                                                                                                                              gCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCACTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTgg  >  4:176251/1‑149 (MQ=255)                                                                                                                                               ccGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTgg  <  4:248711/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                CCACTGCAACGGCCTCAACCGCCTGCCCGCCGAGCGCCTGGGCGATGCGCCGAAATGCGGTCGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGG  >  NC_002947/5786617‑5786906

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