Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,869,613	G→C	5.7%	G41R (GGG→CGG) ‡	PP_5144 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,869,613	0	G	C	5.7%	284.9 / 5.5	88	G41R (GGG→CGG) ‡	PP_5144	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (39/43);  new base C (2/3);  total (42/46)
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.07e-01

ACCAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTGCA  >  NC_002947/5869468‑5869761                                                                                                                                                   |                                                                                                                                                      aCCAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTggg                                                                                                                                                     <  6:71236/148‑1 (MQ=255)   cAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTg                                                                                                                                                  >  7:232239/1‑149 (MQ=255)   cAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTg                                                                                                                                                  >  7:309854/1‑149 (MQ=255)    aaaaaaCACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGAACGCCCTGACGCTGAGCGTGTCGGC‑CAACCGCGACGAATGCTATGCCCGCCCCCCGCAGCCCCTGGAGGCCTGGTAGGATGAACCTGGGGTg                                                                                                                                                  >  8:166928/1‑148 (MQ=255)     aaaaaCACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCa                                                                                                                                                >  3:329502/1‑149 (MQ=255)      aaaaCACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCAt                                                                                                                                               <  6:1681/149‑1 (MQ=255)       aaaCACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATg                                                                                                                                              >  7:15000/1‑149 (MQ=255)       aaaCACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCCCAGGTGCAtc                                                                                                                                              <  3:321867/149‑2 (MQ=255)         acacAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCCCAGGTGCAtcct                                                                                                                                            <  7:9041/149‑4 (MQ=255)          cacaGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGACGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCg                                                                                                                                           >  2:95139/1‑149 (MQ=255)           acaGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCgg                                                                                                                                          >  2:980/1‑149 (MQ=255)            caGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCggg                                                                                                                                         <  6:365305/149‑1 (MQ=255)             agagGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCTCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGCATGCACCCCAGGTGCAtcctccca                                                                                                                                        <  1:42895/149‑8 (MQ=255)               agGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGgcgc                                                                                                                                      >  4:86484/1‑149 (MQ=255)                   ccAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGAcc                                                                                                                                  <  6:128641/149‑1 (MQ=255)                   ccAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGAcc                                                                                                                                  >  6:271191/1‑149 (MQ=255)                      aCCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTg                                                                                                                                >  8:70973/1‑148 (MQ=255)                        cTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCTCGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTTCATGCCGGGCGCGACCTGGag                                                                                                                             <  6:210293/149‑1 (MQ=255)                           tgtgCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGAGTGCATGCCGGGCGCGACCTGGagg                                                                                                                            <  7:126524/147‑1 (MQ=255)                           tgtgCCTGATCGTAGTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGagg                                                                                                                            <  5:50211/147‑1 (MQ=255)                            gtgCCTGATAGTATTTGCCTGGAGGCCGGGGCAAGCCCTGCCGCTGATCGTGGCGGC‑AAACTGCGACGAATTATATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGaggc‑ag                                                                                                                        >  7:153893/1‑149 (MQ=255)                              gCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGaggc‑aggc                                                                                                                      >  3:44945/1‑149 (MQ=255)                               ccTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGt                                                                                                                                                   <  7:10810/120‑1 (MQ=255)                               ccTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGt                                                                                                                                                   >  8:10810/1‑120 (MQ=255)                                 tGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCgg                                                                                                                    >  2:320247/1‑148 (MQ=255)                                        aTTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACctggc                                                                                                            >  4:28603/1‑149 (MQ=255)                                         tttGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCACAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACctggc                                                                                                            >  1:274960/1‑149 (MQ=255)                                             ccTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGc                                                                                                       <  2:5312/149‑1 (MQ=255)                                             ccTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGc                                                                                                       <  1:304023/149‑1 (MQ=255)                                              cTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCg                                                                                                      <  8:88517/149‑1 (MQ=255)                                                 gcggcCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGt                                                                                                     <  4:197975/147‑1 (MQ=255)                                                 gcggcCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTa                                                                                                    <  8:309854/148‑1 (MQ=255)                                                 gcggcCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAg                                                                                                   >  4:331960/1‑149 (MQ=255)                                                  cggtCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAgg                                                                                                  <  1:95139/149‑1 (MQ=255)                                                      cGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCca                                                                                              <  4:295994/149‑1 (MQ=255)                                                       ggggCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTTGGGTGCATGCCGGGCGCGCCCTGGATGC‑AGGCGGTACCTGGCTGGGCGTAGGCCca                                                                                              >  2:329692/1‑148 (MQ=255)                                                       ggggCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCcac                                                                                             <  8:15000/149‑1 (MQ=255)                                                           cACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCCTACGCCcagccag                                                                                         <  7:51445/149‑6 (MQ=255)                                                                 cTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCtt                                                                                   <  3:331960/149‑1 (MQ=255)                                                                   gCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTg                                                                                 >  3:97075/1‑149 (MQ=255)                                                                      gCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCcg                                                                              >  2:19675/1‑149 (MQ=255)                                                                      gCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCcg                                                                              >  7:34756/1‑149 (MQ=255)                                                                          aTCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGgcgc                                                                                     >  4:79431/1‑138 (MQ=255)                                                                          aTCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGgcgc                                                                                     <  3:79431/138‑1 (MQ=255)                                                                             gTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCGCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGa                                                                        <  4:329502/148‑1 (MQ=255)                                                                             gTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGAc                                                                       >  1:86533/1‑149 (MQ=255)                                                                             gTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGAc                                                                       <  8:232239/149‑1 (MQ=255)                                                                                gcggc‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCaa                                                                    <  3:86484/149‑1 (MQ=255)                                                                                       aaCCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCcg                                                              <  5:125444/149‑1 (MQ=255)                                                                                          cgcgACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGa                                                           >  1:61708/1‑149 (MQ=255)                                                                                          cgcgACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGa                                                           >  4:214512/1‑149 (MQ=255)                                                                                            cgacgaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGc                                                                                <  1:334478/126‑1 (MQ=255)                                                                                            cgacgaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGc                                                                                >  2:334478/1‑126 (MQ=255)                                                                                              acgaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCg                                                       <  6:180808/149‑1 (MQ=255)                                                                                               cgaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCgg                                                      <  7:70973/149‑1 (MQ=255)                                                                                               cgaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCgg                                                      <  5:271191/149‑1 (MQ=255)                                                                                                gaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGAc                                                          >  3:4106/1‑144 (MQ=255)                                                                                                gaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGAc                                                          <  4:4106/144‑1 (MQ=255)                                                                                                gaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCggg                                                     >  3:362516/1‑149 (MQ=255)                                                                                                 aaTTCTCTGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGc                                                    <  5:143932/149‑1 (MQ=255)                                                                                                     cTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAg                                                 <  4:97075/148‑1 (MQ=255)                                                                                                      tATGCCCGCCCCACTCCGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGATGC‑AGTCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCTCGCTTACCAACATCCGCGACCCGGGCCAGgc                                               <  7:166928/149‑1 (MQ=255)                                                                                                          cccgcccCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAg                                           >  3:30232/1‑149 (MQ=255)                                                                                                          cccgcccCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTAGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAg                                           >  3:104694/1‑149 (MQ=255)                                                                                                           ccgcccCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAgg                                          <  1:19675/149‑1 (MQ=255)                                                                                                           ccgcccCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAgg                                          >  4:45801/1‑149 (MQ=255)                                                                                                            cgcccCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑GGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAggg                                         >  2:339793/1‑149 (MQ=255)                                                                                                                 cACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGc                                                    >  7:156263/1‑133 (MQ=255)                                                                                                                 cACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGc                                                    <  8:156263/133‑1 (MQ=255)                                                                                                                    tCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCctcgct                                 >  2:306587/1‑149 (MQ=255)                                                                                                                       gCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCtcgctcgc                              <  4:362516/149‑1 (MQ=255)                                                                                                                         cccTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTcgcg                             >  3:316567/1‑148 (MQ=255)                                                                                                                         cccTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCgg                            <  2:86533/149‑1 (MQ=255)                                                                                                                          ccTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCggg                           >  8:60509/1‑149 (MQ=255)                                                                                                                             ggcggcGGGGGAGGATGCACCGGGGGTGCATGCGGGGCGCGACCTGGCGGCTCCCCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGt                         <  1:320247/149‑1 (MQ=255)                                                                                                                                  tcTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGc                     <  3:316159/146‑1 (MQ=255)                                                                                                                                   cTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCtgg                  >  2:119659/1‑149 (MQ=255)                                                                                                                                    tGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGc                     >  4:316159/1‑145 (MQ=255)                                                                                                                                      ggaggaTGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCtggtgg               <  4:44945/149‑1 (MQ=255)                                                                                                                                       cggcaTGCACCCCAGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCTAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCtggtgg               >  2:299810/5‑148 (MQ=255)                                                                                                                                        ggcaTGCACCCCAGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCtggtgg               >  7:189890/4‑147 (MQ=255)                                                                                                                                        aggaTGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAgg                                          >  6:13125/1‑120 (MQ=255)                                                                                                                                        aggaTGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAgg                                          <  5:13125/120‑1 (MQ=255)                                                                                                                                         ggaTGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAg            <  1:119659/149‑1 (MQ=255)                                                                                                                                         ggaTGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAg            <  1:92821/149‑1 (MQ=255)                                                                                                                                          gaTGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑GGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGc           <  1:339793/149‑1 (MQ=255)                                                                                                                                              cACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTCGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGttt       <  4:104694/149‑1 (MQ=255)                                                                                                                                               aCCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTc      >  5:337640/1‑149 (MQ=255)                                                                                                                                                ccTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCt     <  3:28603/149‑1 (MQ=255)                                                                                                                                                   ggggTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTg    >  3:174337/1‑147 (MQ=255)                                                                                                                                                   ggggTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGAGGCAGCGTTTCTGCa  <  4:30232/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      ACCAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGC‑CAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGC‑AGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTGCA  >  NC_002947/5869468‑5869761

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