Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	424,640	A→C	8.5%	W224G (TGG→GGG) ‡	PP_0349 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	424,640	0	A	C	8.5%	207.0 / 8.0	70	W224G (TGG→GGG) ‡	PP_0349	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base C (3/3);  total (34/36)
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.35e-01

GCCTGGAAACGCGATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCTTGTGGCGTACGGTGCCGAAT  >  NC_002947/424497‑424785                                                                                                                                                   |                                                                                                                                                    gCCTGGAAACGCGATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCg                                                                                                                                                >  7:37590/1‑149 (MQ=255)  ccTGGAAACGCGATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCg                                                                                                                                                <  4:89242/148‑1 (MQ=255)   cTGGAAACGCGATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GAGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCACCCATCGCAGCCAGGCTGTCGTTGAGCCAGGCCGc                                                                                                                                               >  4:331259/1‑148 (MQ=255)    tGGAAACGCGATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCaa                                                                                                                                             >  3:56309/1‑149 (MQ=255)    tGGAAACGCGATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCaa                                                                                                                                             >  3:366217/1‑149 (MQ=255)      gAAACGCGATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTACGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATg                                                                                                                                           <  1:263939/149‑1 (MQ=255)       aaaCGCGATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGa                                                                                                                                          <  8:37590/149‑1 (MQ=255)         aCGCGATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCATTGCGGCCTGGCGCAACGAca                                                                                                                                        <  7:339060/149‑2 (MQ=255)             gATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACGATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGt                                                                                                                                    <  5:318402/149‑1 (MQ=255)             gATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGt                                                                                                                                    <  7:153749/149‑1 (MQ=255)                gAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCAGGCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGa                                                                                                                                          >  1:265095/1‑141 (MQ=255)                gAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGa                                                                                                                                          <  2:265095/140‑1 (MQ=255)                   gTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGCAGGCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGa                                                                                                                                >  1:107253/1‑149 (MQ=255)                    tCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAg                                                                                                                             >  8:87086/1‑149 (MQ=255)                    tCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAg                                                                                                                             <  6:221746/149‑1 (MQ=255)                    tCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAg                                                                                                                             >  6:334561/1‑149 (MQ=255)                    tCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAg                                                                                                                             <  1:49299/149‑1 (MQ=255)                      cGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAgg                                                                                                                            <  3:294687/148‑1 (MQ=255)                       gggCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGc                                                                                                                          >  3:292862/1‑149 (MQ=255)                       gggCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCATTGCGGCCTGGCGCAACGAca                                                                                                                                        >  2:11270/1‑134 (MQ=255)                       gggCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGATGTCATTGCGGCCTGGCGCAACGAca                                                                                                                                        <  1:11270/135‑2 (MQ=255)                        ggCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCGGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCa                                                                                                                         <  3:331259/149‑1 (MQ=255)                              tGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGAt                                                                                                                    >  7:309317/1‑148 (MQ=255)                                 ccGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTAcc                                                                                                                >  6:315274/1‑149 (MQ=255)                                 ccGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTAcc                                                                                                                >  2:107712/1‑149 (MQ=255)                                  cGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTAccc                                                                                                               <  6:325563/149‑1 (MQ=255)                                   ggCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTAcccc                                                                                                              <  4:292862/149‑1 (MQ=255)                                      gcgGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGGCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGt                                                                                                           <  5:99204/149‑1 (MQ=255)                                          ccTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACg                                                                                                       >  4:61820/1‑149 (MQ=255)                                           cTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGc                                                                                                      <  3:82666/149‑1 (MQ=255)                                            tGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTAcc                                                                                                                >  5:173082/1‑138 (MQ=255)                                            tGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTAcc                                                                                                                <  6:173082/138‑1 (MQ=255)                                            tGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCt                                                                                                     >  5:287351/1‑149 (MQ=255)                                            tGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCATTGCGGCCTGGCGCAACGAca                                                                                                                                        >  4:319847/1‑113 (MQ=255)                                            tGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCATTGCGGCCTGGCGCAACGAca                                                                                                                                        <  3:319847/114‑2 (MQ=255)                                              ttCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCt                                                                                                     <  5:315274/147‑1 (MQ=255)                                                 ggCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCg                                                                                                >  6:330968/1‑149 (MQ=255)                                                   cAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGt                                                                                               >  7:29077/1‑148 (MQ=255)                                                       gTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCg                                                                                                                                  <  2:129606/109‑1 (MQ=255)                                                        tCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCg                                                                                                <  1:27064/142‑1 (MQ=255)                                                        tCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCg                                                                                                >  2:27064/1‑142 (MQ=255)                                                         ccA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTc                                                                                                          <  4:10468/131‑1 (MQ=255)                                                         ccA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTc                                                                                                          >  3:10468/1‑131 (MQ=255)                                                          ca‑gcagGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGcc                                                                                        >  1:1962/1‑148 (MQ=255)                                                               agGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCccaggc                                                                                   <  2:107253/149‑1 (MQ=255)                                                                          tGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑a                                                                      >  3:348930/1‑149 (MQ=255)                                                                               ccgccgcTCGTCGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGc                                                                   >  3:34830/5‑149 (MQ=255)                                                                                    gcTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTTCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑CCGCCGAt                                                               >  6:197497/1‑148 (MQ=255)                                                                                    gcTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGAt                                                               >  5:61348/1‑148 (MQ=255)                                                                                     cTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCa                                                             <  5:330968/149‑1 (MQ=255)                                                                                      tCGACGAGCGGCGGGGCGTCGAGGTAGGGGGCAACCATCGCAGACAGGCTGTCGTTGCTCCAGGCCTCTATGACCAGGTCGGAGAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCGGATCAg                                                            <  6:366474/149‑1 (MQ=255)                                                                                        gacgaGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAgg                                                           <  1:247922/148‑1 (MQ=255)                                                                                           gaGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGGACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTtgt                                                       >  7:243188/1‑149 (MQ=255)                                                                                             gcggcgGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGc                                                     <  3:186110/149‑1 (MQ=255)                                                                                             gcggcgGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGc                                                     <  3:61820/149‑1 (MQ=255)                                                                                              cggcggGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCa                                                    <  4:366217/149‑1 (MQ=255)                                                                                                gcggGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAgg                                                  <  4:56309/149‑1 (MQ=255)                                                                                                 cggGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCCGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGt                                                 >  3:58056/1‑149 (MQ=255)                                                                                                     gCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAg                                             >  5:153130/1‑149 (MQ=255)                                                                                                           cgGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAgcg                                       >  3:154156/1‑149 (MQ=255)                                                                                                                        ccATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGAt                          >  8:6936/1‑149 (MQ=255)                                                                                                                             gCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCtt                     <  6:61348/149‑1 (MQ=255)                                                                                                                               aGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCTtgt                   <  2:1962/149‑1 (MQ=255)                                                                                                                                gCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCTtgtg                  >  5:58787/1‑149 (MQ=255)                                                                                                                                gCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCAGCGGGATTTCTtgtg                  <  4:350224/149‑1 (MQ=255)                                                                                                                                     gCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCTTGTGGCGTa             >  2:22552/1‑149 (MQ=255)                                                                                                                                        gTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGGACTGATTACCCCGGTCACGCTGACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCTTGTGGCGTACgg          >  1:98940/1‑149 (MQ=255)                                                                                                                                         tCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCAAGCTCACCGTCAATGCCCAGGCCAGGGTGACCTTCCGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCTTGTGGCGTACg           <  1:107712/149‑1 (MQ=255)                                                                                                                                           gTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCTTGTGGCGTACGGTg        <  4:154156/148‑1 (MQ=255)                                                                                                                                                gcCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCTTGTGGCGTACGGTGCCGAAt  <  4:361210/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GCCTGGAAACGCGATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCTTGTGGCGTACGGTGCCGAAT  >  NC_002947/424497‑424785

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