Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,673,177	A→C	12.3%	D164A (GAT→GCT)	PP_2340 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,673,177	0	A	C	12.3%	256.4 / 26.7	97	D164A (GAT→GCT)	PP_2340	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (41/42);  new base C (5/7);  total (47/50)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

GTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCAC  >  NC_002947/2673033‑2673323                                                                                                                                                  |                                                                                                                                                   gTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGcc                                                                                                                                                <  5:42466/149‑1 (MQ=255)      cTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATggg                                                                                                                                           <  8:279225/149‑1 (MQ=255)       tGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATggg                                                                                                                                           >  4:117448/1‑148 (MQ=255)         gAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATggg                                                                                                                                           <  3:117448/146‑1 (MQ=255)         gAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCa                                                                                                                                         <  4:244399/148‑1 (MQ=255)              cTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCAtcgcgc                                                                                                                                   <  7:146867/149‑1 (MQ=255)              cTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCAtcgcgc                                                                                                                                   <  7:342957/149‑1 (MQ=255)              cTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTCCCTGCGCGCTGCCCATGGCAtcgcg                                                                                                                                    >  4:327986/1‑148 (MQ=255)                 tGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCATCGCGCa                                                                                                                                  <  3:386822/147‑1 (MQ=255)                 tGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCATCGCGCAg                                                                                                                                 <  6:50457/148‑1 (MQ=255)                 tGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAg                                                                                                                                 <  1:187193/148‑1 (MQ=255)                   gCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTCCCTGCGCGCTGCCCATGGCATCGCGCAggaa                                                                                                                              <  3:327986/149‑3 (MQ=255)                           ttGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑AACGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGc                                                                                                                      <  1:335939/149‑1 (MQ=255)                              ctgcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCTGCGTTATATCGCCCCCAGCTTCCTGCGCTATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCg                                                                                                                    <  7:259509/148‑1 (MQ=255)                               tgcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCggg                                                                                                                  >  7:387423/1‑149 (MQ=255)                                gcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAgc                                                                                                                                       <  6:332014/127‑1 (MQ=255)                                gcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAgc                                                                                                                                       >  5:332014/1‑127 (MQ=255)                                gcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGc                                                                                                                 >  1:184536/1‑149 (MQ=255)                                 ccGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCg                                                                                                                <  6:375034/149‑1 (MQ=255)                                         ccGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTc                                                                                                        >  6:42842/1‑149 (MQ=255)                                               cTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGt                                                                                                                          <  6:996/125‑1 (MQ=255)                                               cTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGt                                                                                                                          >  5:996/1‑125 (MQ=255)                                                  ggAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCATCGCGCAg                                                                                                                                 >  8:331667/1‑115 (MQ=255)                                                   gAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTAt                                                                                              >  3:131471/1‑149 (MQ=255)                                                   gAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTAt                                                                                              >  7:241678/1‑149 (MQ=255)                                                    aaCTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCATGAGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCCGTCGCCGGTGTAtt                                                                                             >  8:167035/1‑149 (MQ=255)                                                         ccccAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGAt                                                                                                                                                   >  6:242379/1‑90 (MQ=255)                                                         ccccAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCaa                                                                                        >  2:19023/1‑149 (MQ=255)                                                         ccccAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCaa                                                                                        <  8:387423/149‑1 (MQ=255)                                                             aaTTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGc                                                                                    >  6:409110/1‑149 (MQ=255)                                                              aactcgcgAT‑CGATCGCATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCt                                                                                   >  1:154276/5‑149 (MQ=255)                                                                  cgcgATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATg                                                                                                                                                  <  3:154400/82‑1 (MQ=255)                                                                  cgcgATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATg                                                                                                                                                  >  4:154400/1‑82 (MQ=255)                                                                  cgcgATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTg                                                                               >  5:171040/1‑149 (MQ=255)                                                                   gcgatgcgatCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGc                                                                                                                 >  6:2771/1‑114 (MQ=255)                                                                   gcgatgcgatCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGc                                                                                                                 <  5:2771/114‑1 (MQ=255)                                                                     gatgcgatCG‑ATCGCGAGTTCTTCCTGCTTTCCGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAg                                                                            <  3:242327/149‑1 (MQ=255)                                                                       tgcgatCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAAt                                                                                                             <  3:62726/114‑1 (MQ=255)                                                                       tgcgatCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAAt                                                                                                             >  4:62726/1‑114 (MQ=255)                                                                          gatcg‑atcGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTc                                                                       >  2:279175/1‑149 (MQ=255)                                                                           atcg‑atcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCg                                                                      >  6:340809/1‑149 (MQ=255)                                                                            tcg‑atcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGCCAGCGCGCAGCGCCAGTGAGCCGGGCGAATCAGTCGCCCGTGTAt                                                                                              >  6:210777/1‑124 (MQ=255)                                                                            tcg‑atcgCGAGTGCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTAt                                                                                              <  5:210777/124‑1 (MQ=255)                                                                                 tcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCTGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGAt                                                                 <  8:241678/149‑1 (MQ=255)                                                                                 tcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGAt                                                                 >  6:278358/1‑149 (MQ=255)                                                                                 tcgCGAGTTCATCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGAt                                                                 <  4:29490/149‑1 (MQ=255)                                                                                     gAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCg                                                              >  6:112593/1‑148 (MQ=255)                                                                                         ccttcCTGCTTTCAGTGTCGCACTACCAGCGCTATATCGCCCCCAGGTTCCTGCGCGATGCCATGGGCAGCGCGGAGGGCCAGTGAGGCGGGCGAATCATTCGCCGGTGTATTCGCAACCGTTGGTGCAGGTTTCGTGGATTCGAACCt                                                         <  1:219745/148‑1 (MQ=255)                                                                                          cttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCt                                                         <  8:325667/148‑1 (MQ=255)                                                                                           ttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                       >  1:350120/1‑149 (MQ=255)                                                                                           ttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                       >  3:292433/1‑149 (MQ=255)                                                                                           ttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCCGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                       >  6:367394/1‑149 (MQ=255)                                                                                           ttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTCCCTGCGCGCTGCCCATGGCATCGCGCAggaag                                                                                                                             >  6:422798/1‑76 (MQ=21)                                                                                           ttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTCCCTGCGCGCTGCCCATGGCATCGCGCAggaag                                                                                                                             <  5:422798/79‑4 (MQ=21)                                                                                             cttGCTTTCAGCGGCGCACTACCAGCGCTATATCGCCCCCCGCTTCCTGCGCTATGCCTTGGGCAGGGCGCAGGGCCAGTGAGCCTGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCtt                                                        <  2:184536/144‑1 (MQ=255)                                                                                             ccTGGTTTCAGAGTCGCACTACCAGCGGTATATCGCCCCGAGCTTCCTGCGCTTAGCCATGGGCAGCGCGAAGGGCCAGTGAGCGGTGCGAATCAGTCGCCGGTGTATTCGCAACCGGTGGTGCAGGTTTCGt                                                                     <  7:232913/133‑1 (MQ=255)                                                                                             ccTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCtt                                                        <  4:131471/146‑1 (MQ=255)                                                                                             ccTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                       <  4:292433/147‑1 (MQ=255)                                                                                             ccTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                       <  1:19023/147‑1 (MQ=255)                                                                                             ccTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTACTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGt                                                                     >  8:232913/1‑133 (MQ=255)                                                                                              cTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGAc                                                    <  4:380811/149‑1 (MQ=255)                                                                                               tGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGa                                                     >  3:349276/1‑147 (MQ=255)                                                                                                 cTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGc                                                 >  4:127750/1‑149 (MQ=255)                                                                                                       gagGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGc                                           <  6:171040/149‑1 (MQ=255)                                                                                                       gagGCGCACTACCAGCGCTATATCG‑CCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCa                                          >  7:387922/1‑149 (MQ=255)                                                                                                        agGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGc                                           <  2:350120/148‑1 (MQ=255)                                                                                                        agGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCa                                          >  6:2643/1‑149 (MQ=255)                                                                                                         ggCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCAt                                         >  7:195906/1‑149 (MQ=255)                                                                                                          gcgcACTACCAGCGATATATCCCCCCCCGCGTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATc                                        <  2:349671/149‑1 (MQ=255)                                                                                                           cgcACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCa                                       <  8:195906/149‑1 (MQ=255)                                                                                                            gcACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAg                                      >  1:222671/1‑149 (MQ=255)                                                                                                                  ccAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGc                                  <  1:279175/147‑1 (MQ=255)                                                                                                                  ccAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCtt                                <  3:127750/149‑1 (MQ=255)                                                                                                                    aGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGa                              >  4:17232/1‑149 (MQ=255)                                                                                                                     gcgcTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACa                                                   >  5:240185/1‑127 (MQ=255)                                                                                                                     gcgcTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGAc                             <  5:340809/149‑1 (MQ=255)                                                                                                                         tataTCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTgg                         >  2:175996/1‑149 (MQ=255)                                                                                                                           tatCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCtt                                >  4:32912/1‑140 (MQ=255)                                                                                                                           tatCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCtt                                <  3:32912/140‑1 (MQ=255)                                                                                                                             tCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTcc                      >  3:212559/1‑148 (MQ=255)                                                                                                                               gCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTcc                      >  2:141569/1‑146 (MQ=255)                                                                                                                                cccccAGCTTCCTGCGCGTTGCCATGGGAAGCGCGCAGGGCCAGTGAGCCCGGCGAATCAGTCCCCGGTGTATTAGCCACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTGCGGCATCAGCGGCTTGACCTGGTcc                      >  5:84316/1‑145 (MQ=255)                                                                                                                                cccccAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGa                  <  5:278358/149‑1 (MQ=255)                                                                                                                                cccccAGCTTCCCGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGa                  <  5:367394/149‑1 (MQ=255)                                                                                                                                 ccccAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGAt                 <  1:175996/149‑1 (MQ=255)                                                                                                                                   ccAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTg                                                                                  <  2:142972/82‑1 (MQ=255)                                                                                                                                   ccAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTg                                                                                  >  1:142972/1‑82 (MQ=255)                                                                                                                                     aGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCAc             <  5:409110/149‑1 (MQ=255)                                                                                                                                         cccTGCGCGCTGCCCATGGCATCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTgg         >  2:61501/2‑149 (MQ=255)                                                                                                                                           cTGCGCGCTGCCCATGGCATCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTc                                                                                            >  2:232139/1‑64 (MQ=21)                                                                                                                                           cTGCGCGCTGCCCATGGCATCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTc                                                                                            <  1:232139/64‑1 (MQ=21)                                                                                                                                           cTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGc        <  4:212559/148‑1 (MQ=255)                                                                                                                                            tGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGa      >  8:46422/1‑149 (MQ=255)                                                                                                                                              cgcgATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGAt     <  5:42842/148‑1 (MQ=255)                                                                                                                                                cgATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCAc  <  7:46422/149‑1 (MQ=255)                                                                                                                                                 gATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATc    >  5:251632/1‑146 (MQ=255)                                                                                                                                                 gATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATc    <  6:251632/146‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCAC  >  NC_002947/2673033‑2673323

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