Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,717,958	A→C	5.8%	I175L (ATC→CTC)	fabB →	3‑oxoacyl‑ACP synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,717,958	0	A	C	5.8%	283.9 / 4.0	92	I175L (ATC→CTC)	fabB	3‑oxoacyl‑ACP synthase
Reads supporting (aligned to +/- strand):  ref base A (37/49);  new base C (4/2);  total (41/51)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCGTCGGCCCATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTCGCGCGCCTACGACGCTGAC  >  NC_002947/4717810‑4718104                                                                                                                                                     |                                                                                                                                                   gcgcGTCGGCCCATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGa                                                                                                                                                    <  5:102578/149‑1 (MQ=255)        ggCCCATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCTATCTCGTCGGCCTGCGCCACCTCTGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGt                                                                                                                                              <  2:286846/148‑1 (MQ=255)        ggCCCATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTg                                                                                                                                             <  2:296216/149‑1 (MQ=255)        ggCCCATACCGCGTGACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTg                                                                                                                                             >  7:370635/1‑149 (MQ=255)         gCCCATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTgg                                                                                                                                            >  1:249814/1‑149 (MQ=255)         gCCCATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTgg                                                                                                                                            <  8:370635/149‑1 (MQ=255)          cccATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCAGCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTggg                                                                                                                                           <  1:357466/149‑1 (MQ=255)           ccATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTgggg                                                                                                                                          >  5:117693/1‑149 (MQ=255)           ccATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTgggg                                                                                                                                          >  1:436660/1‑149 (MQ=255)             aTACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGc                                                                                                                                         >  1:128891/1‑148 (MQ=255)             aTACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGc                                                                                                                                         >  7:269142/1‑148 (MQ=255)             aTACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACATCGGCTCACTGCAGCGGTACCCCACTGGAGCAGATCCAGTGGGGc                                                                                                                                         >  2:431607/1‑148 (MQ=255)               aCCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTggg                                                                                                                                           <  2:194288/144‑1 (MQ=255)               aCCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTggg                                                                                                                                           >  1:194288/1‑144 (MQ=255)                ccGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAg                                                                                                                                      <  2:225944/148‑1 (MQ=255)                  gcgTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAg                                                                                                                                   <  4:271565/149‑1 (MQ=255)                    gTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGa                                                                                                                                 <  4:83878/149‑1 (MQ=255)                        cccGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATc                                                                                                                             <  6:72187/149‑1 (MQ=255)                         ccGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCGACCGCGTTCATGATGTAGGGCAGCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCGGACTGCAGCGGTACCGCCCTGGAGGAGATCCCGTGGGGCAAGCAGGCCATCg                                                                                                                            >  4:169986/1‑149 (MQ=255)                          cGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGATCTGCTCCAGTgcgg                                                                                                                                          >  8:292421/1‑134 (MQ=255)                          cGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGATCTGCTCCAGTgcgg                                                                                                                                          <  7:292421/134‑1 (MQ=255)                          cGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGt                                                                                                                           <  7:104479/149‑1 (MQ=255)                           gCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTc                                                                                                                          >  5:175405/1‑149 (MQ=255)                            cACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATTGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCt                                                                                                                         >  2:86972/1‑149 (MQ=255)                               cATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCg                                                                                                                      >  3:136400/1‑149 (MQ=255)                                   ggCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCg                                                                                                                   <  3:422108/148‑1 (MQ=255)                                     cagcaCCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTg                                                                                                                <  2:436660/149‑1 (MQ=255)                                        gacggtgtCGGCATGCCGGGCCACCCCGTTCAAGATCAACGGCATCATCTTCTCGATCTCGTCGGCGTGCGCCACCTCGGCTCACTGCATCGGTACCGCCCTGCAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCTGTGGCg                                                                                                             <  8:412129/145‑1 (MQ=255)                                        caCCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTAGTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTAGCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCg                                                                                                             <  5:318315/149‑1 (MQ=255)                                        caCCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCg                                                                                                             <  6:375400/149‑1 (MQ=255)                                        caCCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCg                                                                                                             <  8:108872/149‑1 (MQ=255)                                          ccGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGt                                                                                                           <  5:254657/149‑1 (MQ=255)                                           cGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTg                                                                                                          <  3:136969/149‑1 (MQ=255)                                           cGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTg                                                                                                          <  2:249814/149‑1 (MQ=255)                                             tgtCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGTCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTgaa                                                                                                        <  6:111341/149‑1 (MQ=255)                                                 ggCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGaagaag                                                                                                    <  1:339425/149‑1 (MQ=255)                                                 ggCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGaagaag                                                                                                    <  1:174246/149‑1 (MQ=255)                                                  gCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAagaaga                                                                                                   <  5:125409/149‑1 (MQ=255)                                                   cATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAagag                                                                                                  <  6:175405/149‑1 (MQ=255)                                                    aTGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGc                                                                                                 >  8:334429/1‑149 (MQ=255)                                                       ccTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGGGGCGGTGAAGAAGAGCAc                                                                                               <  8:237479/148‑1 (MQ=255)                                                         tcgCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGACTGCGCCACCTCGCCTCACTCCATCGGTACCGCAATGGAGAACAACAACAGGGGCAACCAGGACATCGTCTTCCCCCGAGGCCCTGAAGAAcaccacag                                                                                             >  8:94063/3‑146 (MQ=255)                                                           gCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAg                                                                                          >  2:87008/1‑149 (MQ=255)                                                            ccACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCTCTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGc                                                                                         <  6:285047/149‑1 (MQ=255)                                                            ccACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGc                                                                                         <  4:336835/149‑1 (MQ=255)                                                                    ttCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGTACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGtt                                                                                 <  3:169986/149‑1 (MQ=255)                                                                     tCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTc                                                                                >  3:408854/1‑149 (MQ=255)                                                                      cAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTc                                                                                <  2:128891/148‑1 (MQ=255)                                                                      cAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTc                                                                                <  1:87008/148‑1 (MQ=255)                                                                       aaGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCAc                                                                                               <  4:383896/132‑1 (MQ=255)                                                                       aaGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCAc                                                                                               >  3:383896/1‑132 (MQ=255)                                                                       aaGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCt                                                                              >  5:276257/1‑149 (MQ=255)                                                                       aaGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCt                                                                              <  8:292589/149‑1 (MQ=255)                                                                       aaGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTAACGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCt                                                                              >  3:287148/1‑149 (MQ=255)                                                                         gATCAAGGGCATCAACTACTCGATCTTGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGCCATCGTCTTAGCCGGTGGCGGTGAAGAAGAGCACTGTAGCCAGACGTTCCTGt                                                                            >  6:285611/1‑149 (MQ=255)                                                                         gATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGt                                                                            >  5:26442/1‑149 (MQ=255)                                                                             aaGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCgt                                                                        >  7:409404/1‑148 (MQ=255)                                                                             aaGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGa                                                                        <  1:431607/149‑1 (MQ=255)                                                                             aaGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACAGCACTGGAGCAGATCCAGTGGGGAAAGCAGGACAACGTCTTCGCCGGTGGCGGTGAAGAAGAGCACAGCAGCCAGTCGTTACTGttaca                                                                        >  7:142661/1‑146 (MQ=255)                                                                                   aTCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCAt                                                                  >  5:245331/1‑149 (MQ=255)                                                                                    tCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATg                                                                 >  7:96/1‑149 (MQ=255)                                                                                    tCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATg                                                                 >  1:288734/1‑149 (MQ=255)                                                                                      aaCTACTCGATCTTGTCGGCCTGCGACAACTCGGCTCACTTCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATgg                                                                <  5:285611/148‑1 (MQ=255)                                                                                       actactCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGATCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGgc                                                              <  7:94063/149‑1 (MQ=255)                                                                                       actactCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGgc                                                              >  6:194856/1‑149 (MQ=255)                                                                                           ctCGATCTCGTCGGCCTGCGCTACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGccc                                                          >  1:255533/1‑149 (MQ=255)                                                                                           ctCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATTGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAGGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGgcgc                                                            <  1:86972/147‑1 (MQ=255)                                                                                           ctCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGcc                                                           <  7:271760/148‑1 (MQ=255)                                                                                            tCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCt                                                         >  7:237623/1‑149 (MQ=255)                                                                                             cGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTg                                                        >  7:280443/1‑149 (MQ=255)                                                                                              gATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCt                                                         <  1:43872/147‑1 (MQ=255)                                                                                                tctcGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCg                                                     >  5:319021/1‑149 (MQ=255)                                                                                                 ctcGTCGGCCTGCCCCACCTAGGCTACCTCCCACGCTACCGCAGTCGACCACCCCCAGGGGGGCAAGCAGGACATCGACTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGa                                                    <  8:379597/149‑1 (MQ=255)                                                                                                    gtcgGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTc                                                      <  6:245331/144‑1 (MQ=255)                                                                                                    gtcgGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCg                                                     <  4:136400/145‑1 (MQ=255)                                                                                                         ccTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAgcgc                                            <  7:110247/149‑1 (MQ=255)                                                                                                             cgcCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCaa                                          <  2:288734/147‑1 (MQ=255)                                                                                                                cACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGaaa                                     >  3:109395/1‑149 (MQ=255)                                                                                                                 aCCTCGGCTCACTGCATCGGTACTGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAAc                                    >  6:13954/1‑149 (MQ=255)                                                                                                                      ggCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCg                                <  7:334429/148‑1 (MQ=255)                                                                                                                           aCTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACt                           <  5:194856/148‑1 (MQ=255)                                                                                                                             tGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGc                        >  3:413220/1‑149 (MQ=255)                                                                                                                               cATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGc                        <  8:237623/147‑1 (MQ=255)                                                                                                                               cATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGGCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGc                        <  6:117693/147‑1 (MQ=255)                                                                                                                                aTCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTc                     <  7:159505/149‑1 (MQ=255)                                                                                                                                    gTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAACCCCCGGAACTGGCCTcgcgc                 >  8:140383/1‑149 (MQ=255)                                                                                                                                       ccccACTGGATCTGCTCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTCgcgcgc               >  8:260630/4‑148 (MQ=255)                                                                                                                                         ccACTGGATCTGCTCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAAGCCCGGAACTGGCCTCGCGCGCCTa            >  6:379890/2‑149 (MQ=255)                                                                                                                                         ccACTGGATCTGCTCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTCGCGCGCCTa            >  7:62931/2‑149 (MQ=255)                                                                                                                                               gAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTCGCGCGCCTacgacg       <  8:280443/148‑1 (MQ=255)                                                                                                                                                 gCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTCGCGCGCCTACGACGCTg    <  6:276257/149‑1 (MQ=255)                                                                                                                                                 gCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGAAAAGAAACCCCGCAACTGGCCTCGCGCGCCTACTACGCt     >  8:383952/1‑148 (MQ=255)                                                                                                                                                   aGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAACAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTCGCCCGCCTACGACGCTGAc  >  1:338698/1‑149 (MQ=255)                                                                                                                                                     aTCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTcgcgcg                >  8:202377/1‑133 (MQ=255)                                                                                                                                                     aTCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTcgcgcg                <  7:202377/133‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GCGCGTCGGCCCATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTCGCGCGCCTACGACGCTGAC  >  NC_002947/4717810‑4718104

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