Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,715,776	C→T	6.3%	A241V (GCA→GTA)	PP_1510 →	carboxyvinyl‑carboxyphosphonate phosphorylmutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,715,776	0	C	T	6.3%	185.9 / 2.0	63	A241V (GCA→GTA)	PP_1510	carboxyvinyl‑carboxyphosphonate phosphorylmutase
Reads supporting (aligned to +/- strand):  ref base C (30/29);  new base T (2/2);  total (32/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 = 2.89e-01

GTGGAGGAAGTCCGTGCCGTGGTGCAGGCAGTCGCACCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTCCGCCGCTGATGAGGGCATTGCTGGCGCTACTGGCCGGCCTGCTGACGG  >  NC_002947/1715628‑1715924                                                                                                                                                      |                                                                                                                                                     gTGGAGGAAGTCCGTGCCGTGGTGCAGGCAGTCGCACCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGtgctgc                                                                                                                                                      <  5:332181/149‑1 (MQ=255)      ggaAGTCCGTGCCGTGGTGCAGGCAGTCGCACCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCg                                                                                                                                                 >  4:110498/1‑149 (MQ=255)      ggaAGTCCGTGCCGTGGTGCAGGCAGTCGCACCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCg                                                                                                                                                 >  3:301778/1‑149 (MQ=255)      ggaAGTCCGTGCCGTGGTGCAGGCAGTCGCACCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCCCTGG‑CCCGTGCTGCACTCg                                                                                                                                                 >  8:171464/1‑149 (MQ=255)       gaAGTCCGTGCCGTGGTGCAGGCAGTCGCACCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCgg                                                                                                                                                >  6:80269/1‑149 (MQ=255)        aaGTCCGTGCCGTGGTGCAGGCAGTCGCACCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCggg                                                                                                                                               <  5:354741/149‑1 (MQ=255)             cgtgccgtgGTGCAGGCAGTCGCACCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCtt                                                                                                                                          <  8:166625/149‑1 (MQ=255)              gtgccgtgGTGCAGGCAGTCGCACCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTc                                                                                                                                         <  3:124165/149‑1 (MQ=255)               tgccgtgGTGCAGGCAGTCGCACCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCt                                                                                                                                        <  8:242803/149‑1 (MQ=255)               tgccgtgGTGCAGGCAGTCGCACCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCGGCCCTCGGGGCCTTc                                                                                                                                         >  6:71379/1‑148 (MQ=255)                            gcagTCGCACCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCGTTCTCCCCGgcggcg                                                                                                                             >  8:116053/1‑147 (MQ=255)                            gcagTCGCACCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCa                                                                                                                           <  8:277495/149‑1 (MQ=255)                              agTCGCACCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAggcg                                                                                                                                <  8:136250/142‑1 (MQ=255)                              agTCGCACCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAggcg                                                                                                                                >  7:136250/1‑142 (MQ=255)                                  gCACCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGaaa                                                                                                                      <  1:259414/148‑1 (MQ=255)                                   cACCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAAtt                                                                                                                    >  4:193690/1‑149 (MQ=255)                                    aCCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTc                                                                                                                   <  3:110498/149‑1 (MQ=255)                                      cgcgCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTcgc                                                                                                                 >  2:337360/1‑149 (MQ=255)                                           ccTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGa                                                                                                            <  6:319567/149‑1 (MQ=255)                                           ccTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGa                                                                                                            <  1:12041/149‑1 (MQ=255)                                               tgAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAgg                                                                                                          <  4:301778/147‑1 (MQ=255)                                                     tgCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTc                                                                                                  <  1:265630/149‑1 (MQ=255)                                                       cTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTc                                                                                                  >  2:265630/1‑147 (MQ=255)                                                       cTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCg                                                                                                 <  1:190009/148‑1 (MQ=255)                                                        tGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTcgc                                                                                                                 <  8:99276/131‑1 (MQ=255)                                                        tGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTcgc                                                                                                                 >  7:99276/1‑131 (MQ=255)                                                              gCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGGCCCCGTGCTGCACTCGGGGCCTTCTAc                                                                                                                                      <  1:3593/105‑1 (MQ=255)                                                              gCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTAc                                                                                                                                      >  2:3593/1‑104 (MQ=255)                                                              gCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTAGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGt                                                                                         >  6:61246/1‑149 (MQ=255)                                                                 tGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAg                                                                                      >  1:74818/1‑149 (MQ=255)                                                                  ggcgggcgTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAg                                                                                      >  5:218219/1‑148 (MQ=255)                                                                       gcgTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCGACGGTGAGCAGGc                                                                                 >  5:231017/1‑148 (MQ=255)                                                                                  aGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGttt                                                                       <  5:61246/147‑1 (MQ=255)                                                                                     gTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGATGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATCAGGGCACGTTCGGCTACGGTGAGCAGGCAAAGCAGTTTGCCCa                                                                  >  3:311016/1‑149 (MQ=255)                                                                                        aaCCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTc                                                                                                  >  7:318479/1‑114 (MQ=255)                                                                                        aaCCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTc                                                                                                  <  8:318479/114‑1 (MQ=255)                                                                                        aaCCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCt                                                               <  6:218219/149‑1 (MQ=255)                                                                                        aaCCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCt                                                               <  3:193690/149‑1 (MQ=255)                                                                                        aaCCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCt                                                               >  3:77396/1‑149 (MQ=255)                                                                                         aCCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCt                                                               >  4:366211/1‑148 (MQ=255)                                                                                          ccAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTg                                                              >  7:33799/1‑148 (MQ=255)                                                                                                 cAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCt                                                      <  2:74818/149‑1 (MQ=255)                                                                                                  aGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTg                                                     >  1:2415/1‑149 (MQ=255)                                                                                                     aCCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTc                                                  <  7:171464/149‑1 (MQ=255)                                                                                                       cTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGGAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTccg                                                >  2:4916/1‑149 (MQ=255)                                                                                                                cgccgcGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCt                                                               <  5:309009/125‑1 (MQ=255)                                                                                                                cgccgcGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCt                                                               >  6:309009/1‑125 (MQ=255)                                                                                                                 gccgcGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTCCGCCGCtgat                                        >  5:212601/1‑147 (MQ=255)                                                                                                                  ccgcGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTCCGCCGCTGATGAg                                     <  2:2415/149‑1 (MQ=255)                                                                                                                       tCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGttt                                                                       >  2:202695/1‑110 (MQ=255)                                                                                                                       tCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCCGCTACGGTGAGCAGGCAATGCCGttt                                                                       <  1:202695/110‑1 (MQ=255)                                                                                                                        cAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGTACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTCCGCCGCtgatga                                      <  3:361838/142‑1 (MQ=255)                                                                                                                        cAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGTACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTCCGCCGCtgatga                                      >  4:367750/1‑142 (MQ=255)                                                                                                                        cAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGTACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTCCGCCGCtgatga                                      >  4:361838/1‑142 (MQ=255)                                                                                                                        cAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGTACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTCCGCCGCtgatga                                      <  3:367750/142‑1 (MQ=255)                                                                                                                         aGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTCCGCCGCTGATGAGGGCATTg                              <  7:210573/149‑1 (MQ=255)                                                                                                                         aGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTCCGCCGCTGATGAGGGCATTg                              >  1:254970/1‑149 (MQ=255)                                                                                                                             ttGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTCCGCCGCTGATGAGGGCATTGCTgg                          >  3:37117/1‑149 (MQ=255)                                                                                                                                    gtcgCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTAGCCCAGCTGAACGACCTGTTCCGCCGCTGATGAGGGCATTGCTGGCGCTAc                    <  2:303050/148‑1 (MQ=255)                                                                                                                                        cTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTCCGCCGCTGATGAGGGCATTGCTGGCGCTACTggcc               >  7:198342/1‑149 (MQ=255)                                                                                                                                           g‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTCCGCCGCTGATGAGGGCATTGCTGGCGCTACTggccgtc            >  8:156486/1‑147 (MQ=255)                                                                                                                                           g‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTCCGCCGCTGATGAGGGCATTGCTGGCGCTACTggccg              <  6:231017/147‑1 (MQ=255)                                                                                                                                           g‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGAGGAGGGCACGGTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTCCGCCGCTGATGAGGGCATTGCTGGCGCTACTggccggc            <  5:71379/149‑1 (MQ=255)                                                                                                                                             cccGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTCCGCCGCTGATGAGGGCATTGCTGGCGCTACTggccggcc           <  3:366211/149‑1 (MQ=255)                                                                                                                                                      cACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTCCGCCGCTGATGAGGGCATTGCTGGCGCTACTGGCCGGCCTGCTGACgg  >  7:6031/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GTGGAGGAAGTCCGTGCCGTGGTGCAGGCAGTCGCACCGCGCCCTGTGAATGTGCTGATGGGCATGGCGGGCGTGCCGCTGAGCGTCAACCAGTTGCAGGACCTGGGCGTGCGCCGCGTCAGCGTTGGCTCGTCGCTGG‑CCCGTGCTGCACTCGGGGCCTTCTACCAGGCGGCGCAGGAAATTCGCGATGAGGGCACGTTCGGCTACGGTGAGCAGGCAATGCCGTTTGCCCAGCTGAACGACCTGTTCCGCCGCTGATGAGGGCATTGCTGGCGCTACTGGCCGGCCTGCTGACGG  >  NC_002947/1715628‑1715924

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