Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,644,840	G→C	5.2%	P825R (CCG→CGG) ‡	PP_2316 ←	ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,644,840	0	G	C	5.2%	238.9 / 2.5	76	P825R (CCG→CGG) ‡	PP_2316	ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (33/39);  new base C (2/2);  total (35/41)
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.64e-01

GGTGAAGCGACCTGGTATGTGAGCAGCATCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACTGCAAATCGAAGGCATAGCGGTACAGCACCCAGGTAATCAGTTCGCA  >  NC_002947/2644692‑2644988                                                                                                                                                     |                                                                                                                                                     ggTGAAGCGACCTGGTATGTGAGCAGCATCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCg                                                                                                                                                      <  1:374563/149‑1 (MQ=255)  gTGAAGCGACCTGGTATGTGAGCAGCATCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCgg                                                                                                                                                     >  4:278777/1‑149 (MQ=255)       gCGACCTGGTATGTGAGCAGCATCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTgg                                                                                                                                                >  5:226399/1‑149 (MQ=255)       gCGACCTGGTATGTGAGCAGCATCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCCCGCTgg                                                                                                                                                >  3:322786/1‑149 (MQ=255)          aCCTGGTATGTGAGCAGCATCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCAt                                                                                                                                             >  3:43406/1‑149 (MQ=255)          aCCTGGTATGTGAGCAGCATCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCCCGCTGGcgg                                                                                                                                             <  4:322786/149‑3 (MQ=255)           ccTGGTATGTGAGCAGCATCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCAtt                                                                                                                                            <  8:319019/149‑1 (MQ=255)            cTGGTATGTGAGCAGCATCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCAtt                                                                                                                                            >  8:45373/1‑148 (MQ=255)               gTATGTGAGCAGCATCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCaa                                                                                                                                         <  8:126999/148‑1 (MQ=255)                  tgtgAGCAGCATCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAAcgc                                                                                                                                      <  7:45373/148‑1 (MQ=255)                    tgAGCAGCATCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAAcgcccg                                                                                                                                   >  8:72869/1‑149 (MQ=255)                     gAGCAGCATCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAAcgcccgc                                                                                                                                  >  5:207759/1‑149 (MQ=255)                     gAGCAGCATCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAAcgcccgc                                                                                                                                  >  3:52438/1‑149 (MQ=255)                         agcaTCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCgtg                                                                                                                              <  4:43406/149‑1 (MQ=255)                          gcaTCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCgtgt                                                                                                                             <  1:332883/149‑1 (MQ=255)                            aTCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTc                                                                                                                           >  8:57525/1‑149 (MQ=255)                             tCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGc                                                                                                                                                   >  5:256462/1‑124 (MQ=255)                             tCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGc                                                                                                                                                   <  6:256462/124‑1 (MQ=255)                             tCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTcc                                                                                                                          >  7:26900/1‑149 (MQ=255)                             tCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTcc                                                                                                                          >  3:24955/1‑149 (MQ=255)                             tCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTcc                                                                                                                          >  4:161342/1‑149 (MQ=255)                              ccGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAAc                                                                                                                                        <  5:370722/134‑1 (MQ=255)                              ccGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAAc                                                                                                                                        >  6:370722/1‑134 (MQ=255)                               cGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTccca                                                                                                                        <  4:295925/149‑1 (MQ=255)                               cGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTccca                                                                                                                        <  6:226399/149‑1 (MQ=255)                               cGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTccca                                                                                                                        <  7:113462/149‑1 (MQ=255)                               cGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTccc                                                                                                                         <  3:278777/148‑1 (MQ=255)                                gTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTcccag                                                                                                                       <  4:293962/149‑1 (MQ=255)                                    cGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTcccagcac                                                                                                                    <  8:26900/148‑1 (MQ=255)                                     ggCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCCCGCTGGcgg                                                                                                                                             >  6:182210/1‑120 (MQ=255)                                     ggCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCCCGCTGGcgg                                                                                                                                             <  5:182210/122‑3 (MQ=255)                                      gCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTcccagcaccc                                                                                                                  >  5:9228/1‑148 (MQ=255)                                       ccATGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTcccagcacccca                                                                                                                <  3:124118/149‑1 (MQ=255)                                       ccAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTcccagcaccc                                                                                                                  >  1:393257/1‑147 (MQ=255)                                            gCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGTGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTcc                                                                                                                          <  6:394016/134‑1 (MQ=255)                                            gCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGCGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTcc                                                                                                                          >  5:394016/1‑134 (MQ=255)                                               tctcGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCaccacc                                                                                                        <  7:72869/149‑1 (MQ=255)                                                 tcGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCCGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCCGCACCCCAGCAccaccca                                                                                                      >  5:28448/1‑149 (MQ=255)                                                   gACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCAccacccac                                                                                                     >  6:400910/1‑148 (MQ=255)                                                   gACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCAccacccac                                                                                                     <  6:11040/148‑1 (MQ=255)                                                    aCGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACccacca                                                                                                   <  4:70994/149‑1 (MQ=255)                                                      ggACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACcag                                                                                                  >  2:1733/1‑148 (MQ=255)                                                       gactgactGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACcagca                                                                                                <  3:203354/149‑1 (MQ=255)                                                          tgactgGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTg                                                                                             >  1:62420/1‑149 (MQ=255)                                                            actgGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGcc                                                                                           >  7:18283/1‑149 (MQ=255)                                                            actgGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGc                                                                                            <  5:400910/148‑1 (MQ=255)                                                                gCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGccccc                                                                                        <  4:52438/148‑1 (MQ=255)                                                                     tCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGc                                                                                  <  1:57955/149‑1 (MQ=255)                                                                      cTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAg                                                                                   >  5:383625/1‑147 (MQ=255)                                                                      cTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAg                                                                                   <  6:383625/147‑1 (MQ=255)                                                                         gCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGc                                                                               >  2:64908/1‑148 (MQ=255)                                                                              aGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAacac                                                                          <  6:9228/148‑1 (MQ=255)                                                                                   gAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTccca                                                                                                                        <  2:347761/97‑1 (MQ=255)                                                                                   gAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTccca                                                                                                                        >  1:347761/1‑97 (MQ=255)                                                                                       cATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCAc                                                                <  4:24955/149‑1 (MQ=255)                                                                                           tgCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGt                                                            <  1:1733/149‑1 (MQ=255)                                                                                           tgCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGt                                                            <  3:161342/149‑1 (MQ=255)                                                                                              tGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTg                                                                                                                                                 >  5:262711/1‑61 (MQ=255)                                                                                              tGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTg                                                                                                                                                 <  6:262711/61‑1 (MQ=255)                                                                                                  ccGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCt                                                     >  6:385456/1‑149 (MQ=255)                                                                                                   cGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACcagcag                                                                    <  3:267325/133‑1 (MQ=255)                                                                                                   cGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCAACACCACTGCCCCCCTCAGCGGCAACACcagcag                                                                    >  4:267325/1‑133 (MQ=255)                                                                                                       aGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACt                                                >  8:305005/1‑149 (MQ=255)                                                                                                          aGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACTGc                                              <  6:207759/148‑1 (MQ=255)                                                                                                                gACTGCCAGGCCTCAAGCCGCGCTCAACACCGCCCGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACTGCAAATCGa                                       <  8:305162/149‑1 (MQ=255)                                                                                                                   tGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACTGCAAATCGAAgg                                    >  6:59653/1‑149 (MQ=255)                                                                                                                       aGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACTGCAAATCGAAGGCa                                  >  8:99092/1‑147 (MQ=255)                                                                                                                        ggCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACTGCAAATCGAAGGCATAg                               <  7:305005/149‑1 (MQ=255)                                                                                                                         gCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGAGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACTGCAAATCGAAGGCATAGc                              <  5:61073/149‑1 (MQ=255)                                                                                                                           cTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACTGCAAATCGAAGGCATAGAgg                            >  4:274546/1‑149 (MQ=255)                                                                                                                            tCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACTGCAAATCGAAGGCATAGCgg                            >  7:314489/1‑148 (MQ=255)                                                                                                                            tCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACTGCAAATCGAAGGCATAGCgg                            >  6:198639/1‑148 (MQ=255)                                                                                                                            tCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACTGCAAATCGAAGGCATAGCg                             >  4:352609/1‑147 (MQ=255)                                                                                                                            tCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACTGCAAATCGAAGGCATAGCg                             <  3:352609/147‑1 (MQ=255)                                                                                                                                       gcAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACTGCAAATCGAAGGCATAGCGGTACAGCACCCAg                <  7:57525/149‑1 (MQ=255)                                                                                                                                       gcAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACTGCAAATCGAAGGCATAGCGGTACAGCACCCAg                <  2:62420/149‑1 (MQ=255)                                                                                                                                       gcAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACTGCAAATCGAAGGCATAGCGGTACAGCACCCAg                <  2:393257/149‑1 (MQ=255)                                                                                                                                                 cAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACTGCAAATCGAAGGCATAGCGGTACAGCACCCAGGTAATCAGtt      <  5:59653/149‑1 (MQ=255)                                                                                                                                                     gggCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACTGCAAATCGAAGGCATAGCGGTACAGCACCCAGGTAATCAGTTCGCa  <  1:64908/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GGTGAAGCGACCTGGTATGTGAGCAGCATCCGTTACGGCCAGGGCGTCTCGACGGACTGACTGGCCCGTCTTGCGTCAGCTTGAGGCATGTGCTGTCCCGGCAGAAGCGCTGACTGCCAGGCCTCAAGCCTCGCGCAACACCGCCAGCGGGCTGGCATTCAACGCCCGCCGTGTTCCCAGCACCCCAGCACCACCCACCAGCAGTGCCCCCGTCAGCGGCAACACCAGCAGCCACGGGTGCGGGCTCCACTGCAAATCGAAGGCATAGCGGTACAGCACCCAGGTAATCAGTTCGCA  >  NC_002947/2644692‑2644988

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