Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,453,977	A→C	5.9%	D204A (GAT→GCT)	PP_1272 →	multidrug MFS transporter membrane fusion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,453,977	0	A	C	5.9%	212.7 / 4.2	68	D204A (GAT→GCT)	PP_1272	multidrug MFS transporter membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base A (34/29);  new base C (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

CCGAGCAGGGCGCGGGTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCGGGCGCTCGCCTGCTGTCG  >  NC_002947/1453830‑1454125                                                                                                                                                     |                                                                                                                                                     ccGAGCAGGGCGCGGTTACCGTGCAGGACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑At                                                                                                                                                     >  1:5541/1‑149 (MQ=255) ccGAGCAGGGCGCGGGTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑At                                                                                                                                                     <  1:282984/149‑1 (MQ=255)   gAGCAGGGCGCGGGTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑a                                                                                                                                                      >  4:26869/1‑146 (MQ=255)   gAGCAGGGCGCGGGTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑a                                                                                                                                                      <  3:26869/146‑1 (MQ=255)       aGGGCGCGGGTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAg                                                                                                                                               >  4:215633/1‑149 (MQ=255)          gcgcGGGTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTc                                                                                                                                            <  8:280950/149‑1 (MQ=255)            gcgGGTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGACCATCGGCG‑CTGGCCacttgag                                                                                                                                          <  5:40597/149‑7 (MQ=255)              gggTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCATCGGCG‑CTGGCCactt                                                                                                                                             >  2:351320/1‑141 (MQ=255)              gggTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCATCGGCG‑CTGGCCactt                                                                                                                                             <  1:351320/144‑4 (MQ=255)                      tGCAGAACACCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCg                                                                                                                                  >  8:137367/1‑147 (MQ=255)                       gCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCt                                                                                                                                             <  6:273539/135‑1 (MQ=255)                       gCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCt                                                                                                                                             >  5:273539/1‑135 (MQ=255)                        cAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGc                                                                                                                              <  4:76923/149‑1 (MQ=255)                          gAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAg                                                                                                                            >  5:30406/1‑149 (MQ=255)                               cccAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTg                                                                                                                       <  6:297236/149‑1 (MQ=255)                                   gcagGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTg                                                                                                                       <  2:365846/145‑1 (MQ=255)                                   gcagGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTg                                                                                                                       >  1:365846/1‑145 (MQ=255)                                    cagGAGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGaa                                                                                                                    >  6:42096/1‑147 (MQ=255)                                     agGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAg                                                                                                                 <  6:359258/149‑1 (MQ=255)                                      ggCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAgg                                                                                                                <  8:77706/149‑1 (MQ=255)                                        cgcgcAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGcc                                                                                                              >  1:55426/1‑149 (MQ=255)                                         gcgcAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGccc                                                                                                             <  6:275400/149‑1 (MQ=255)                                                 cTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGAt                                                                                                     <  7:63788/149‑1 (MQ=255)                                                      gATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTc                                                                                                >  4:232051/1‑149 (MQ=255)                                                         cAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGAt                                                                                                     >  4:4318/1‑141 (MQ=255)                                                         cAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGAt                                                                                                     <  3:4318/141‑1 (MQ=255)                                                           aGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTaca                                                                                           >  7:45033/1‑149 (MQ=255)                                                               cGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCg                                                                                       <  2:82691/149‑1 (MQ=255)                                                                 tGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCa                                                                                     <  8:88952/149‑1 (MQ=255)                                                                      gCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATca                                                                                  <  7:273112/147‑1 (MQ=255)                                                                      gCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTCTCCTACACCCGCATcaca                                                                                <  1:361180/149‑1 (MQ=255)                                                                              aaCGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGt                                                                        >  8:285219/1‑149 (MQ=255)                                                                              aaCGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGt                                                                        <  5:42096/149‑1 (MQ=255)                                                                              aaCGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCAAGCTCGATCTGTCATACACCCGCATACCAGCGCCGGt                                                                        >  3:235633/1‑149 (MQ=255)                                                                                 gCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGa                                                                     <  7:285219/149‑1 (MQ=255)                                                                                   tCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCaa                                                                                                                                          >  6:390147/1‑78 (MQ=255)                                                                                     aaGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGTCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACgg                                                                  >  2:237754/1‑148 (MQ=255)                                                                                     aaGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACgg                                                                  >  1:26658/1‑148 (MQ=255)                                                                                             cgGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCgg                                                         <  3:232051/149‑1 (MQ=255)                                                                                                  ggcggcGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCGAATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGaa                                                     >  1:284579/1‑149 (MQ=255)                                                                                                  ggcggcGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑TTGGCCAGCTGAAGCGTGCCGAAGCAGGCCTGGAATCGGCCCATCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTTGACGGTATGGTCGGCGAAc                                                    <  8:152413/149‑1 (MQ=255)                                                                                                   gcggcgCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAAcg                                                   >  2:58203/1‑149 (MQ=255)                                                                                                           cAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGt                                                                        <  8:282599/120‑1 (MQ=255)                                                                                                           cAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGt                                                                        >  7:282599/1‑120 (MQ=255)                                                                                                              gCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGaa                                                                                                                    <  5:327374/73‑1 (MQ=255)                                                                                                              gCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGaa                                                                                                                    >  6:327374/1‑73 (MQ=255)                                                                                                                  gTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTggg                                    <  2:55426/149‑1 (MQ=255)                                                                                                                     gACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGc                                 >  5:268214/1‑149 (MQ=255)                                                                                                                        aTCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTa                              >  6:218010/1‑149 (MQ=255)                                                                                                                         tCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTAc                             >  8:36824/1‑149 (MQ=255)                                                                                                                          ccTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTc                                                           >  7:175003/1‑118 (MQ=255)                                                                                                                          ccTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTc                                                           <  8:175003/118‑1 (MQ=255)                                                                                                                            tCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGt                           >  1:99670/1‑148 (MQ=255)                                                                                                                            tCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGt                           >  3:140255/1‑148 (MQ=255)                                                                                                                                ggCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAAcc                      >  4:88467/1‑149 (MQ=255)                                                                                                                                ggCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAAcc                      >  4:1711/1‑149 (MQ=255)                                                                                                                                   tCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCgg                   >  5:109815/1‑149 (MQ=255)                                                                                                                                    gagctGGCCATCGGCG‑CTGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCggg                  >  1:221274/5‑149 (MQ=255)                                                                                                                                     aaGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCggg                  <  7:180179/148‑1 (MQ=255)                                                                                                                                     aaGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCgg                   <  1:206517/147‑1 (MQ=255)                                                                                                                                     aaGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCgg                   <  4:140255/147‑1 (MQ=255)                                                                                                                                     aaGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCGGgc                 >  2:89236/1‑149 (MQ=255)                                                                                                                                     aaGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCGGgc                 <  2:284579/149‑1 (MQ=255)                                                                                                                                        tGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCGGgcgc               >  4:87198/1‑148 (MQ=255)                                                                                                                                        tGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCGGgcgc               <  3:87198/148‑1 (MQ=255)                                                                                                                                        tGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCGGGCGCt              >  8:189815/1‑149 (MQ=255)                                                                                                                                          gCCAGCGCCG‑ATGGTCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCGGGCGCTCg            <  1:237754/149‑1 (MQ=255)                                                                                                                                          gCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCCACCCGGGCGCTCg            >  8:64214/1‑149 (MQ=255)                                                                                                                                           ccAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCGGGCGCTCGc           <  2:308132/149‑1 (MQ=255)                                                                                                                                                   g‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCGGGCGCTCGCCTGCTGTc   >  8:328350/1‑149 (MQ=255)                                                                                                                                                     aTGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGGCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCGGGCGCTCGCCTGCTGTCg  <  7:137367/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CCGAGCAGGGCGCGGGTACCGTGCAGAACGCCCAGCAGGCGCGCAGCGCTGTCGATCAAGCCCGTGCTCGCCTGGCCAACGCTCAAGCTGCGCGGGTGGCGGCGCGCAAGCAGGTGGACATCCTCACGGCTCAAGTGGCCAGCGCCG‑ATGGCCAGCTCAAGCGTGCCGAAGCAGGCCTGGAAAAGGCCCAGCTCGATCTGTCCTACACCCGCATCACAGCGCCGGTGGACGGTATGGTCGGCGAACGCGCCCTGCGCGTGGGGGCCTACGTCAACCCGGGCGCTCGCCTGCTGTCG  >  NC_002947/1453830‑1454125

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑