Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,812,659	A→G	6.7%	D328G (GAC→GGC)	PP_3368 →	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,812,659	0	A	G	6.7%	173.5 / 2.8	60	D328G (GAC→GGC)	PP_3368	MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (27/29);  new base G (2/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

GCGCCGTGGGCTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGCGCCGCGACTG  >  NC_002947/3812512‑3812805                                                                                                                                                    |                                                                                                                                                   gcgcCGTGGGCTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGa                                                                                                                                                    >  8:350075/1‑148 (MQ=255) gcgcCGTGGGCTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGa                                                                                                                                                    >  4:329134/1‑148 (MQ=255)          gCTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCAt                                                                                                                                          <  7:32744/149‑1 (MQ=255)              ggaaCTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAgg                                                                                                                                      <  6:244535/146‑1 (MQ=255)               gTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGGCAATGTCCATGAgg                                                                                                                                      <  1:259763/148‑1 (MQ=255)                     aTGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAg                                                                                                                               <  5:351538/149‑1 (MQ=255)                              gtcgCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAgc                                                                                                                      >  6:239663/1‑149 (MQ=255)                               tcgCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAgcg                                                                                                                     <  7:350075/149‑1 (MQ=255)                                  ccGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGACGCCTAGCCTGGCAgcggcg                                                                                                                  >  7:19413/1‑149 (MQ=255)                                     tgctgcTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTcc                                                                                                               <  3:71673/149‑1 (MQ=255)                                        tgctgTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAAc                                                                                                            <  7:308126/149‑1 (MQ=255)                                         gctgTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAAc                                                                                                            >  6:26927/1‑148 (MQ=255)                                          ctgTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAAc                                                                                                            <  6:160887/147‑1 (MQ=255)                                           tgTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCAt                                                                                                         >  6:391235/1‑149 (MQ=255)                                               cTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATgcg                                                                                                      <  3:323216/148‑1 (MQ=255)                                                tttACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATgcggc                                                                                                    >  7:130391/1‑149 (MQ=255)                                                  tACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATgcggc                                                                                                    <  3:329134/147‑1 (MQ=255)                                                     cTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCCGCGGCGTCCAACCATGCGGCGTTc                                                                                                >  3:376321/1‑148 (MQ=255)                                                      ttCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCaa                                                                                              >  3:318743/1‑149 (MQ=255)                                                         gccgccCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAAcc                                                                                            <  2:131253/148‑1 (MQ=255)                                                          ccgccCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTg                                                                                          >  6:81275/1‑149 (MQ=255)                                                             cccAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGcc                                                                                       <  5:239663/149‑1 (MQ=255)                                                              ccAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTTGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGcc                                                                                       >  8:284231/1‑148 (MQ=255)                                                               cAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCaa                                                                                     >  4:238415/1‑149 (MQ=255)                                                                aGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAAt                                                                                    <  8:130391/149‑1 (MQ=255)                                                                        tgGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGGCTTGGAGGCAGCGCTACAGATGCGCCTGATGGACTTTGCCCATGAGGCGCCTTGCCTGGCAGCGGGGTCCAACCATGCGGCGTTCAAGCTGGCCAATGCCCttg                                                                             <  6:314538/148‑3 (MQ=255)                                                                        tgGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCCCCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCATTGCCCTggg                                                                            <  5:391235/149‑1 (MQ=255)                                                                        tgGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTggg                                                                            <  5:26927/149‑1 (MQ=255)                                                                        tgGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTg                                                                              <  6:101799/147‑1 (MQ=255)                                                                          gAGCNTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGc                                                                          >  5:109889/1‑149 (MQ=255)                                                                                gctgctgGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGc                                                                    >  5:262551/1‑149 (MQ=255)                                                                                 ctgctgGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCt                                                                   <  4:262577/149‑1 (MQ=255)                                                                                  tgctgGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTggg                                                                            <  1:151236/139‑1 (MQ=255)                                                                                  tgctgGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTggg                                                                            >  2:151236/1‑139 (MQ=255)                                                                                        gCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGc                                                            <  2:214237/149‑1 (MQ=255)                                                                                                   gtggGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCTGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCAc                                                 >  1:123110/1‑149 (MQ=255)                                                                                                     gggCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAg                                                              >  3:218792/1‑134 (MQ=255)                                                                                                     gggCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAg                                                              <  4:218792/134‑1 (MQ=255)                                                                                                        cACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCg                                            >  1:151431/1‑149 (MQ=255)                                                                                                         aCCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGGCTAGCCTGTCAGCGGCGTCCAACCTTGCGGCGTTACAACTGGCCAATCGCCTGCGGCGGTGGCTGGGAGGCATGGCGATCACTGCCgg                                           >  5:327141/1‑149 (MQ=255)                                                                                                         aCCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGGCGTGGCTGGGAGGCAAGGCGATCACTGCCgg                                           >  5:175397/1‑149 (MQ=255)                                                                                                         aCCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCgg                                           >  5:171228/1‑149 (MQ=255)                                                                                                         aCCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCgg                                           >  4:222790/1‑149 (MQ=255)                                                                                                          ccATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCgg                                           >  2:200/1‑148 (MQ=255)                                                                                                           cATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCa                                         >  7:293893/1‑149 (MQ=255)                                                                                                               aTCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATgg                                      <  1:346585/148‑1 (MQ=255)                                                                                                                       ggcggcACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGAc                               <  2:151431/147‑1 (MQ=255)                                                                                                                        gcggcACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCa                             <  3:3930/148‑1 (MQ=255)                                                                                                                          ggcACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCa                          >  7:69504/1‑149 (MQ=255)                                                                                                                          ggcACAGCTACAGCTCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGAGGCGTACAACCTGGCCAAAGAACTGGGGCCGTCGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCa                          >  8:162511/1‑149 (MQ=255)                                                                                                                               cGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGAcc                              >  4:135361/1‑140 (MQ=255)                                                                                                                               cGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGAcc                              <  3:135361/140‑1 (MQ=255)                                                                                                                               cGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGa                     <  6:109889/149‑1 (MQ=255)                                                                                                                               cGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGa                     <  5:81275/149‑1 (MQ=255)                                                                                                                                gCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGAt                    <  4:376321/149‑1 (MQ=255)                                                                                                                                gCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGAt                    <  4:318743/149‑1 (MQ=255)                                                                                                                                      gATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCg              <  3:138315/149‑1 (MQ=255)                                                                                                                                           gCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGcgcc         >  5:181767/1‑149 (MQ=255)                                                                                                                                            ccTCATGGGCAATGTCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGcgccg        >  4:307045/1‑149 (MQ=255)                                                                                                                                             cTCATGGGCAATGTCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGcgccgc       <  3:307045/149‑1 (MQ=255)                                                                                                                                              tCATGGGCAATGTCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGcgccgc       >  3:9970/1‑148 (MQ=255)                                                                                                                                                  ggACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGCGCCGCGACTg  >  4:37456/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCGCCGTGGGCTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGCGCCGCGACTG  >  NC_002947/3812512‑3812805

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