Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,549,921	A→T	6.8%	H418L (CAC→CTC) ‡	PP_3133 →	oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,549,921	0	A	T	6.8%	173.9 / 3.3	59	H418L (CAC→CTC) ‡	PP_3133	oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (34/21);  new base T (2/2);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

TGCTTAGCGCCGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGTTGCAA  >  NC_002947/3549774‑3550069                                                                                                                                                    |                                                                                                                                                     tGCTTAGCGCCGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCGAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCGGAGAATGGCACCCCATGCCCATca                                                                                                                                                      >  8:173306/1‑148 (MQ=255)   cTTAGCGCCGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATcaca                                                                                                                                                    <  1:13441/148‑1 (MQ=255)           cGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATGTGATGGGcatgg                                                                                                                                           <  2:355521/149‑4 (MQ=255)           cGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGcacca                                                                                                                                           >  3:121678/1‑149 (MQ=255)           cGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGcacca                                                                                                                                           >  7:307151/1‑149 (MQ=255)                   ggCACACCGTCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATg                                                                                                                                   <  4:74612/149‑1 (MQ=255)                    gCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATgg                                                                                                                                  >  1:345211/1‑149 (MQ=255)                     cacaCCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGgc                                                                                                                                 <  2:265288/149‑1 (MQ=255)                      acacCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCAGATGGGCACCACGCGCATGgcc                                                                                                                                >  4:346244/1‑149 (MQ=255)                          cATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCt                                                                                                                            >  7:207861/1‑149 (MQ=255)                              aaGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTcgc                                                                                                                         <  6:203994/148‑1 (MQ=255)                                           gCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGgt                                                                                                              >  2:33660/1‑146 (MQ=255)                                           gCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGgt                                                                                                              <  1:33660/146‑1 (MQ=255)                                              aaCTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGt                                                                                                         >  6:160825/1‑148 (MQ=255)                                                   gCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGc                                                                                                   >  1:38558/1‑149 (MQ=255)                                                   gCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGc                                                                                                   <  7:173306/149‑1 (MQ=255)                                                   gCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGc                                                                                                   <  1:288379/149‑1 (MQ=255)                                                     cAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGcc                                                                                                  >  3:218571/1‑148 (MQ=255)                                                     cAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGc                                                                                                   <  7:202861/147‑1 (MQ=255)                                                     cAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGc                                                                                                   >  8:202861/1‑147 (MQ=255)                                                         cAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATGTGATGGGcatggg                                                                                                                                          <  4:167491/104‑5 (MQ=255)                                                         cAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATGTGATGGGcatggg                                                                                                                                          >  3:167491/1‑100 (MQ=255)                                                           gTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCa                                                                                           <  2:345211/149‑1 (MQ=255)                                                              cGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGCCACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAgg                                                                                        >  6:21391/1‑149 (MQ=255)                                                              cGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAgg                                                                                        >  3:55159/1‑149 (MQ=255)                                                                   aCATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCaa                                                                                                <  5:260890/136‑1 (MQ=255)                                                                   aCATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGCCACCCCATGACCATCACATCGGCACCACCCCCATCCCCGCCTCGACGCAGTTATGTGTGGACGATGACaa                                                                                                >  6:260890/1‑136 (MQ=255)                                                                           cTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCaa                                                                                          <  2:212388/134‑1 (MQ=255)                                                                           cTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCaa                                                                                          >  1:212388/1‑134 (MQ=255)                                                                             ggggTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGaa                                                                          >  3:100863/1‑148 (MQ=255)                                                                               ggTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAAc                                                                       <  5:361759/149‑1 (MQ=255)                                                                                gTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCACTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAAcc                                                                      >  2:279469/1‑149 (MQ=255)                                                                                   cGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTg                                                                    <  1:280057/148‑1 (MQ=255)                                                                                            cAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCg                                                          >  7:325579/1‑149 (MQ=255)                                                                                             accGATTTCGTCTCCGACACCCCGATACCGCTGAGAATGGCACCCCCTGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGcc                                                           <  8:191105/145‑1 (MQ=255)                                                                                               cGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGTGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCgg                                                         >  1:192273/1‑147 (MQ=255)                                                                                                    tCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGgcagcag                                                   >  6:101818/1‑148 (MQ=255)                                                                                                     cGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGgcagcagca                                                 >  4:102816/1‑149 (MQ=255)                                                                                                      gTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCAt                                                >  4:102541/1‑149 (MQ=255)                                                                                                       tttACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAg                                                                                         >  4:275643/1‑107 (MQ=255)                                                                                                       tttACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAg                                                                                         <  3:275643/107‑1 (MQ=255)                                                                                                       tttACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCAtt                                               >  3:358238/1‑149 (MQ=255)                                                                                                          aCGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCAtttt                                             <  8:207861/148‑1 (MQ=255)                                                                                                           cGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTg                                           >  2:349760/1‑149 (MQ=255)                                                                                                                 cTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAgg                                     >  8:145882/1‑149 (MQ=255)                                                                                                                 cTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAgg                                     >  3:361464/1‑149 (MQ=255)                                                                                                                   cGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGggc                                   >  6:86701/1‑149 (MQ=255)                                                                                                                         cGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGcc                             <  1:349760/149‑1 (MQ=255)                                                                                                                           cTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTc                           <  8:325579/149‑1 (MQ=255)                                                                                                                               gaATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTTCCAAGGGCGGCGCCTCCAAc                       <  4:55159/149‑1 (MQ=255)                                                                                                                               gaATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGACAAGGGCGGCGCCTCCAAc                       >  5:184234/1‑149 (MQ=255)                                                                                                                                  tGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAAccc                     >  3:284430/1‑148 (MQ=255)                                                                                                                                    gCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGa                   <  5:86701/148‑1 (MQ=255)                                                                                                                                          ggtGCCCATGTGATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGACAATTGCAAGGTGTACGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCAAAGGGCGGCGCCTCCAACCCGACCATGcc            >  7:341782/3‑149 (MQ=255)                                                                                                                                          cATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCCATTGCAAGGTGTTCGGTACCGAAAACCTGTATCTCGCAGGCAGCAGCATTTTTGCCAACGGCGGCGCCTCCAACCCGACCATGcc            >  7:176488/1‑149 (MQ=255)                                                                                                                                           aTGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCg           <  2:38558/149‑1 (MQ=255)                                                                                                                                           aTGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGAATTTTTGCCAAGGGCGGTGCCTCCAACCCGACCATGCCg           >  3:22469/1‑149 (MQ=255)                                                                                                                                             gCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGAAGCATTTTTGCCAAGGGCGGCGCCTCCCACCCGACCATGCCg           >  8:147651/1‑147 (MQ=255)                                                                                                                                                cATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGtt      >  4:350366/1‑149 (MQ=255)                                                                                                                                                cATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGtt      <  4:218571/149‑1 (MQ=255)                                                                                                                                                  tCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGTTg     >  2:3017/1‑148 (MQ=255)                                                                                                                                                    acaTGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGTTGCaa  >  7:94258/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TGCTTAGCGCCGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGTTGCAA  >  NC_002947/3549774‑3550069

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