Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,502,237	T→G	7.7%	intergenic (+29/‑135)	yccA → / → xanP	FtsH‑mediated proteolysis substrate/modulator/xanthine permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,502,237	0	T	G	7.7%	190.4 / 3.3	62	intergenic (+29/‑135)	yccA/xanP	FtsH‑mediated proteolysis substrate/modulator/xanthine permease
Reads supporting (aligned to +/- strand):  ref base T (29/28);  new base G (2/3);  total (31/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 = 1.00e+00

ATTCATGGTGGCGAACGTAACTACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACAC  >  NC_002947/4502095‑4502383                                                                                                                                                |                                                                                                                                                   attcatGGTGGCGAACGTAACTACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTcggcgg                                                                                                                                              <  5:328504/149‑1 (MQ=255)     atGGTGGCGAACGTAACTACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCtt                                                                                                                                          >  2:48862/1‑149 (MQ=255)        gtggCGAACGTAACTACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCtttt                                                                                                                                        <  3:284032/148‑1 (MQ=255)        gtggCGAACGTAACTACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGGTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGAATTGCTGTACCGCAAGAAGCCCGCCGAAGCGGGCttc                                                                                                                                         <  1:62073/148‑2 (MQ=255)            cGAACGTAACTACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTtgtg                                                                                                                                   <  2:61097/149‑1 (MQ=255)            cGAACGTAACTACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTtgtg                                                                                                                                   <  2:111134/149‑1 (MQ=255)                      tACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCCGAAGCGGGCttc                                                                                                                                         >  5:285498/1‑132 (MQ=255)                      tACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCCGAAGCGGGCttc                                                                                                                                         <  6:285498/133‑2 (MQ=255)                          tcatGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGATGTACCGCAAGAACCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATgg                                                                                                                     >  3:13851/1‑149 (MQ=255)                            atGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATggg                                                                                                                    <  3:328993/148‑1 (MQ=255)                             tGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGtt                                                                                                                  >  1:173226/1‑149 (MQ=255)                             tGGCAACCATCAGCCTGTATGTGGCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGAGTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGtt                                                                                                                  <  1:165165/149‑1 (MQ=255)                               gCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAAAAGCCCGCCGAAGCGGGcttcttg                                                                                                                                      >  5:247669/1‑126 (MQ=255)                               gCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAAAAGCCCGCCGAAGCGGGcttcttg                                                                                                                                      <  6:247669/127‑2 (MQ=255)                                        aGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAgg                                                                                                        <  7:173526/148‑1 (MQ=255)                                          ccTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTtgt                                                                                                                                    >  5:326829/1‑118 (MQ=255)                                          ccTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTtgt                                                                                                                                    <  6:326829/118‑1 (MQ=255)                                          ccTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTtgt                                                                                                                                    <  6:331725/118‑1 (MQ=255)                                          ccTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTtgt                                                                                                                                    >  5:331725/1‑118 (MQ=255)                                                    tCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGt                                                                                           >  2:193818/1‑149 (MQ=255)                                                         cTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATcc                                                                                      >  2:227980/1‑149 (MQ=255)                                                           aCAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCg                                                                                    >  4:102446/1‑149 (MQ=255)                                                           aCAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCg                                                                                    >  2:248432/1‑149 (MQ=255)                                                                cTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGGAGGATGCTCGTGCCGTCATCCCGTTTtc                                                                               <  1:193818/149‑1 (MQ=255)                                                                   tttGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTtctc                                                                             >  5:2672/1‑148 (MQ=255)                                                                           ccTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCTCTa                                                                     >  2:334904/1‑148 (MQ=255)                                                                               ctgcAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGc                                                                <  2:173226/149‑1 (MQ=255)                                                                                tgcAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGc                                                                >  2:178381/1‑148 (MQ=255)                                                                                  cAATTGTTTGGCATCTTGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGcc                                                               <  2:196534/147‑1 (MQ=255)                                                                                  cAATTGTTTGGCATCTTGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGcc                                                               >  1:196534/1‑147 (MQ=255)                                                                                  cAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCa                                                              >  1:96422/1‑148 (MQ=255)                                                                                  cAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCAt                                                             <  6:2672/149‑1 (MQ=255)                                                                                    aTTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCa                                                           <  7:226141/149‑1 (MQ=255)                                                                                      tGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATgg                                                                                                                     <  2:166982/89‑1 (MQ=255)                                                                                      tGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATgg                                                                                                                     >  1:166982/1‑89 (MQ=255)                                                                                       gTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACAGCTATCCGCCATCAccc                                                        >  4:295842/1‑149 (MQ=255)                                                                                        tttGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCAcccc                                                       >  1:13336/1‑149 (MQ=255)                                                                                           ggCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGt                                                     <  1:248432/148‑1 (MQ=255)                                                                                            gCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCg                                                      >  8:23544/1‑146 (MQ=255)                                                                                            gCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTa                                                    <  1:48862/148‑1 (MQ=255)                                                                                             catcatGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTtc                                                                               <  8:90446/120‑1 (MQ=255)                                                                                             catcatGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTtc                                                                               >  7:90446/1‑120 (MQ=255)                                                                                              atcatGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGaa                                                 >  3:59928/1‑149 (MQ=255)                                                                                                     gTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTc                                          >  1:101861/1‑149 (MQ=255)                                                                                                         tGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCttt                                      >  8:61191/1‑149 (MQ=255)                                                                                                            cgacTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCtttttt                                   <  1:334904/149‑1 (MQ=255)                                                                                                              acTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTttct                                 >  8:203752/1‑149 (MQ=255)                                                                                                                        ttgCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGggc                         <  4:59928/148‑1 (MQ=255)                                                                                                                          gCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGggcgg                       <  2:101861/148‑1 (MQ=255)                                                                                                                            tGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGggc                         <  8:240557/144‑1 (MQ=255)                                                                                                                            tGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGggc                         >  7:240557/1‑144 (MQ=255)                                                                                                                              tACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCttt                   <  7:61191/148‑1 (MQ=255)                                                                                                                                ccGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAg                <  2:96422/149‑1 (MQ=255)                                                                                                                                   cAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCaa             >  3:256874/1‑149 (MQ=255)                                                                                                                                    aagaagCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAAt            >  3:163080/1‑149 (MQ=255)                                                                                                                                      gaagCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGa          <  1:227980/149‑1 (MQ=255)                                                                                                                                        agCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGc        >  5:352725/1‑149 (MQ=255)                                                                                                                                           ccGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGt                                                                                                <  1:342843/58‑1 (MQ=255)                                                                                                                                           ccGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGt                                                                                                >  2:342843/1‑58 (MQ=255)                                                                                                                                            cGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATcc                                                                  <  1:98533/87‑1 (MQ=255)                                                                                                                                            cGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATcc                                                                  >  2:98533/1‑87 (MQ=255)                                                                                                                                              cTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCacac  <  7:203752/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   ATTCATGGTGGCGAACGTAACTACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTG‑ATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCACAC  >  NC_002947/4502095‑4502383

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