Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,893,962	T→A	6.6%	Q366L (CAG→CTG)	yegH ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,893,962	0	T	A	6.6%	235.7 / 3.5	81	Q366L (CAG→CTG)	yegH	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (35/40);  new base A (4/2);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTGCGCCAG  >  NC_002947/3893816‑3894107                                                                                                                                                   |                                                                                                                                                  tGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctgg                                                                                                                                                 <  3:294483/149‑1 (MQ=255) tGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGCTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctgg                                                                                                                                                 <  1:316481/149‑1 (MQ=255)    tCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGgctggctg                                                                                                                                              >  8:217216/1‑149 (MQ=255)     cAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGAGTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGCAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGc                                                                                                                                             >  4:235694/1‑149 (MQ=255)     cAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGc                                                                                                                                             >  6:158559/1‑149 (MQ=255)      aaCAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGcc                                                                                                                                            >  8:216629/1‑149 (MQ=255)       aCAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGccc                                                                                                                                           >  8:96112/1‑149 (MQ=255)          gCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCAGCCAGCCCgac                                                                                                                                        <  3:360952/149‑2 (MQ=255)           cccGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCgaga                                                                                                                                       <  4:306536/149‑1 (MQ=255)           cccGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCga                                                                                                                                         >  7:254928/1‑147 (MQ=255)            ccGGTGAAGTCCCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGagag                                                                                                                                      <  8:200477/149‑1 (MQ=255)             cGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGTGTCGACTGACCGCGAATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGTTCGGGCTGGATGCCCGagag                                                                                                                                      >  6:359058/1‑148 (MQ=255)             cGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagc                                                                                                                                     <  7:216629/149‑1 (MQ=255)                tGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagt                                                                                                                                  <  7:96112/149‑1 (MQ=255)                   aGTCACCGAACTCGTTGACCACGAAGGCAATATGGGGCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcc                                                                                                                               <  1:286862/149‑1 (MQ=255)                    gTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcct                                                                                                                              >  2:381432/1‑149 (MQ=255)                    gTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCAGCCAGCCCgacctgga                                                                                                                                   >  6:391697/1‑138 (MQ=255)                    gTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCAGCCAGCCCgacctgga                                                                                                                                   <  5:391697/144‑7 (MQ=255)                       aCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttca                                                                                                                           <  6:339183/149‑1 (MQ=255)                       aCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttc                                                                                                                            <  4:345114/148‑1 (MQ=255)                       aCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcctt                                                                                                                             <  8:254928/147‑1 (MQ=255)                        ccGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcag                                                                                                                          <  5:370813/149‑1 (MQ=255)                         cGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCgag                                                                                                                                        >  1:143983/1‑134 (MQ=255)                         cGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCgag                                                                                                                                        <  2:143983/134‑1 (MQ=255)                              tCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcagcagttc                                                                                                                    >  1:140403/1‑149 (MQ=255)                                gTTGACCACGAAGTCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGATGATGATGAAGCTCTCCAGCAGGTTAAAAGGCGCGCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcagcagttcc                                                                                                                   <  5:359058/148‑1 (MQ=255)                                       aCGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTtgtg                                                                                                           >  1:187806/1‑149 (MQ=255)                                         gAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGca                                                                                                         <  3:70215/149‑1 (MQ=255)                                            ggCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGcaca                                                                                                       >  1:222175/1‑148 (MQ=255)                                             gCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACaaa                                                                                                     <  6:215377/149‑1 (MQ=255)                                               aaTATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGc                                                                                                   >  3:192099/1‑149 (MQ=255)                                               aaTATGGGTAGAATGACAGCGCATTTGATCCACGGCATTCAGGATGCTGAAGCTCTCCAGCAGGTTCAAAGGCGCCCGTGCCATGCTTTCCAGGTCTGGCAGGCTGCCCGAGAGCAGTTCCTACAGCAGTTCCATCTTTAGCACaaaca                                                                                                   >  8:319976/1‑147 (MQ=255)                                                  atGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGCGCCATGCTTTCCACGACGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGAACAAAGccc                                                                                                 >  1:385688/1‑148 (MQ=255)                                                      gTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAggg                                                                                            >  4:627/1‑149 (MQ=255)                                                                 gcgcATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCAcgc                                                                                 >  6:260608/1‑149 (MQ=255)                                                                    cATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGAc                                                                              <  2:222175/149‑1 (MQ=255)                                                                     aTTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGAcc                                                                             <  3:291289/149‑1 (MQ=255)                                                                     aTTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGAc                                                                              >  8:122916/1‑148 (MQ=255)                                                                           tCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCACCGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGGCCACGCGACCATcgcg                                                                       >  8:305219/1‑149 (MQ=255)                                                                           tCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATcgcg                                                                       >  5:107815/1‑149 (MQ=255)                                                                             cAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATcgcg                                                                       <  5:158559/147‑1 (MQ=255)                                                                             cAGGGATTTGAGGATGCTGAAGCTCTCAAGCAGGTTCAACGGCGCCCCGGCCATGCTTTCCTGGTCGGGCTGGCTGCCGGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATcgc                                                                        <  7:319976/146‑1 (MQ=255)                                                                              aGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGAt                                                                    >  8:61460/1‑149 (MQ=255)                                                                               gggCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGAt                                                                    >  6:2984/1‑148 (MQ=255)                                                                                ggCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCa                                                                  <  3:19923/149‑1 (MQ=255)                                                                                 gCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAg                                                                 <  2:187806/149‑1 (MQ=255)                                                                                         gATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGgc                                                         >  3:197672/1‑149 (MQ=255)                                                                                           tGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAGAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGgcgc                                                       <  7:217216/149‑1 (MQ=255)                                                                                            gCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGgcgc                                                       >  4:161618/1‑148 (MQ=255)                                                                                             cTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGgcgc                                                       >  3:69929/1‑147 (MQ=255)                                                                                             cTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGa                                                     <  3:235694/149‑1 (MQ=255)                                                                                              tGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAg                                                    >  8:64233/1‑149 (MQ=255)                                                                                              tGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAg                                                    >  7:241244/1‑149 (MQ=255)                                                                                               gAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGa                                                     <  3:354395/147‑1 (MQ=255)                                                                                               gAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGa                                                     >  4:354395/1‑147 (MQ=255)                                                                                                aaGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAg                                                    <  7:64233/147‑1 (MQ=255)                                                                                                aaGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGt                                                   <  7:324560/148‑1 (MQ=255)                                                                                                 aGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTAc                                                 <  5:2984/149‑1 (MQ=255)                                                                                                    tctcCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGc                                              >  2:47718/1‑149 (MQ=255)                                                                                                      tcCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATcgcg                                                                       <  3:241338/122‑1 (MQ=255)                                                                                                      tcCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATcgcg                                                                       >  4:241338/1‑122 (MQ=255)                                                                                                          gcagGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTg                                        <  2:140403/149‑1 (MQ=255)                                                                                                          gcagGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTg                                        <  4:192099/149‑1 (MQ=255)                                                                                                           cagGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGGTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTgg                                       <  7:305219/149‑1 (MQ=255)                                                                                                              gTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCa                                    <  7:122916/149‑1 (MQ=255)                                                                                                                 cAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAggg                                 >  7:255167/1‑149 (MQ=255)                                                                                                                 cAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAgg                                  <  2:385688/148‑1 (MQ=255)                                                                                                                  aaCGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGc                                <  4:191911/149‑1 (MQ=255)                                                                                                                     ggCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTg                             <  8:255167/149‑1 (MQ=255)                                                                                                                      gcgcCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTgg                            <  4:197672/149‑1 (MQ=255)                                                                                                                         cccGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGc                          >  5:348273/1‑148 (MQ=255)                                                                                                                         cccGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGc                          >  7:100842/1‑148 (MQ=255)                                                                                                                          ccGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCATTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGTCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCaa                        <  8:302816/149‑1 (MQ=255)                                                                                                                          ccGGGCCATGCTTTCCAGGTCGGGCTGGCTGCACGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCaa                        >  7:341286/1‑149 (MQ=255)                                                                                                                              gCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGCCCAGGCCCTGGCCAAttcg                    <  6:107815/149‑4 (MQ=255)                                                                                                                                cATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGt                  >  1:344640/1‑149 (MQ=255)                                                                                                                                    cTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCATCGGc              <  8:241244/149‑1 (MQ=255)                                                                                                                                     tttCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGcc             >  2:10473/1‑149 (MQ=255)                                                                                                                                        ccAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCgg           <  3:161618/148‑1 (MQ=255)                                                                                                                                          aGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTg        >  3:50841/1‑149 (MQ=255)                                                                                                                                           tctCGGGCAGCCAGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgc       >  5:377436/3‑149 (MQ=255)                                                                                                                                           tctCGGGCAGCCAGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgc       >  7:151415/3‑149 (MQ=255)                                                                                                                                             tCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgcgc     <  1:381432/149‑1 (MQ=255)                                                                                                                                             tCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgcgc     <  3:627/149‑1 (MQ=255)                                                                                                                                                ggctggctGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTGCGCCAg  >  3:39792/1‑149 (MQ=255)                                                                                                                                                ggctggctGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTGCGCCAg  <  5:260608/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTGCGCCAG  >  NC_002947/3893816‑3894107

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