Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,814,595	C→G	6.3%	intergenic (+95/‑18)	frmC → / → ispF	S‑formylglutathione hydrolase/S‑lactoylglutathione hydrolase/2‑C‑methyl‑D‑erythritol 2,4‑cyclodiphosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,814,595	0	C	G	6.3%	201.2 / 3.7	64	intergenic (+95/‑18)	frmC/ispF	S‑formylglutathione hydrolase/S‑lactoylglutathione hydrolase/2‑C‑methyl‑D‑erythritol 2,4‑cyclodiphosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (26/34);  new base G (2/2);  total (28/36)
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.12e-01

CGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACGGTGACGTGCTGCTGCACG  >  NC_002947/1814455‑1814739                                                                                                                                             |                                                                                                                                                 cGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAggg                                                                                                                                          >  5:58850/1‑149 (MQ=255) cGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAggg                                                                                                                                          >  7:192022/1‑149 (MQ=255)  gaggagCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGc                                                                                                                                         <  4:289386/149‑1 (MQ=255)    ggagCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGAATGAGTCAGGGCGt                                                                                                                                       <  5:82907/149‑1 (MQ=255)      agCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGttt                                                                                                                                     <  2:280898/149‑1 (MQ=255)       gCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGtttt                                                                                                                                    <  3:3661/149‑1 (MQ=255)        gaCCTGCGCCATGATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCATCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTc                                                                                                                                   <  7:23472/148‑1 (MQ=255)         aCCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGAATGAGTCAGGGCGtgattc                                                                                                                                  <  4:42325/149‑6 (MQ=255)          ccTGCGCCATCATGCGGTGGCACTGGGGCGGGTCTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCt                                                                                                                                  >  5:79839/1‑148 (MQ=255)            tGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTAt                                                                                                                               >  6:164028/1‑149 (MQ=255)              cgcCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTAt                                                                                                                               <  1:251847/147‑1 (MQ=255)              agcCATGATGCGGTGGCACTGGGGCGGGTTTAACGCCTCGTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTAt                                                                                                                               <  8:187796/146‑1 (MQ=255)                   tcatGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTAtt                                                                                                                        >  3:269343/1‑149 (MQ=255)                               aCTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTAc                                                                                                            >  2:14145/1‑149 (MQ=255)                               aCTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTAc                                                                                                            >  2:31847/1‑149 (MQ=255)                               aCTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTAc                                                                                                            >  3:100665/1‑149 (MQ=255)                               aCTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTAc                                                                                                            >  7:30947/1‑149 (MQ=255)                                 tGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGa                                                                                                          >  1:33615/1‑149 (MQ=255)                                         gTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCAcc                                                                                                  <  2:188107/149‑1 (MQ=255)                                          tttAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGAGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCg                                                                                                 >  1:221444/1‑149 (MQ=255)                                              aCGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTc                                                                                             >  3:248736/1‑149 (MQ=255)                                                gCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTg                                                                                           <  6:117526/149‑1 (MQ=255)                                                 cctcctGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGc                                                                                          <  2:221444/149‑1 (MQ=255)                                                 cctcctGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGc                                                                                          <  4:61699/149‑1 (MQ=255)                                                     ctGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACg                                                                                      <  6:79839/149‑1 (MQ=255)                                                        tgGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGt                                                                                    <  7:260000/148‑1 (MQ=255)                                                        tgGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTg                                                                                   <  7:126939/149‑1 (MQ=255)                                                          ggAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGAt                                                                                 >  2:275926/1‑149 (MQ=255)                                                           gAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGAtt                                                                                <  4:100665/149‑1 (MQ=255)                                                             gCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTc                                                                              <  5:164028/149‑1 (MQ=255)                                                              cGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCa                                                                             <  7:294041/149‑1 (MQ=255)                                                               gggCGTGCCCGCGAAGGGCTTCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTc                                                                              <  6:225075/147‑1 (MQ=255)                                                                ggCGTGCCCGCGAATGGCTGCAAGGCAGACCCAATCCCACAAATCCACTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCAtt                                                                           <  2:276842/149‑1 (MQ=255)                                                                      cccGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCt                                                                      <  8:30947/148‑1 (MQ=255)                                                                       ccGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTgg                                                                    <  3:310265/149‑1 (MQ=255)                                                                       ccGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTgg                                                                    <  3:68103/149‑1 (MQ=255)                                                                       ccGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTgg                                                                    >  3:70015/1‑149 (MQ=255)                                                                       ccGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTgg                                                                    >  3:176944/1‑149 (MQ=255)                                                                            aaGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCggg                                                               <  4:176944/149‑1 (MQ=255)                                                                             aGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCgggg                                                              <  5:205218/149‑1 (MQ=255)                                                                             aGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCgggg                                                              <  3:64406/149‑1 (MQ=255)                                                                                gCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCGTTACCCTGGGCGGGGTGc                                                           <  2:190247/149‑1 (MQ=255)                                                                                gCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGc                                                           >  2:125300/1‑149 (MQ=255)                                                                                        gcagcCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATcccc                                                   <  8:192022/149‑1 (MQ=255)                                                                                        gcagcCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATcccc                                                   <  2:33615/149‑1 (MQ=255)                                                                                             cccAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGCCTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGAGTATCCCCcaca                                               >  2:26033/1‑148 (MQ=255)                                                                                              ccAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCAcc                                                                                                  <  8:289958/96‑1 (MQ=255)                                                                                              ccAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCAcc                                                                                                  >  7:289958/1‑96 (MQ=255)                                                                                              ccAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACaaa                                             >  4:275286/1‑149 (MQ=255)                                                                                                  tCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACg                                         <  1:31847/149‑1 (MQ=255)                                                                                                        aaaTCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCctgc                                   <  3:275286/149‑1 (MQ=255)                                                                                                           tCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTgg                                <  6:58850/149‑1 (MQ=255)                                                                                                           tCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTgg                                >  8:106464/1‑149 (MQ=255)                                                                                                             cccTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGcc                              >  5:2233/1‑149 (MQ=255)                                                                                                              ccTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGcc                              >  4:197082/1‑148 (MQ=255)                                                                                                                      gCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCct                                     >  8:154050/1‑133 (MQ=255)                                                                                                                       cAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCCGGCCCACTCCGACg                    <  4:248736/149‑1 (MQ=255)                                                                                                                        aGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACgg                   >  5:317436/1‑149 (MQ=255)                                                                                                                              aaTCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATcccc                                                   >  5:306869/1‑111 (MQ=255)                                                                                                                              aaTCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATcccc                                                   <  6:306869/111‑1 (MQ=255)                                                                                                                                 aaCGCCCTGAATGAGTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACGGTGACGtgct          >  7:243054/2‑149 (MQ=255)                                                                                                                                  aCGCCCTGAATGAGTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACGGTGACGtgctg         >  6:171225/1‑149 (MQ=255)                                                                                                                                     cccTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACGGTGACGtgctgctg      <  1:14145/149‑1 (MQ=255)                                                                                                                                         gACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACGGTGACGTGCTGCTGCACg  >  7:302030/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                 CGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACGGTGACGTGCTGCTGCACG  >  NC_002947/1814455‑1814739

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