Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,549,920	C→G	6.8%	H418D (CAC→GAC) ‡	PP_3133 →	oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,549,920	0	C	G	6.8%	182.8 / 3.9	59	H418D (CAC→GAC) ‡	PP_3133	oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (34/21);  new base G (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 = 1.89e-01

GGTGCTTAGCGCCGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGTTG  >  NC_002947/3549772‑3550066                                                                                                                                                     |                                                                                                                                                   ggTGCTTAGCGCCGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATc                                                                                                                                                    <  4:158899/149‑1 (MQ=255) ggTGCTTAGCGCCGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATc                                                                                                                                                    <  4:154380/149‑1 (MQ=255)   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:288379/149‑1 (MQ=255)                                                     gCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGc                                                                                                <  7:173306/149‑1 (MQ=255)                                                     gCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGc                                                                                                >  1:38558/1‑149 (MQ=255)                                                       cAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGcc                                                                                               >  3:218571/1‑148 (MQ=255)                                                       cAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGc                                                                                                >  8:202861/1‑147 (MQ=255)                                                       cAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGc                                                                                                <  7:202861/147‑1 (MQ=255)                                                           cAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATGTGATGGGcatggg                                                                                                                                       >  3:167491/1‑100 (MQ=255)                                                           cAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATGTGATGGGcatggg                                                                                                                                       <  4:167491/104‑5 (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                                                                                       >  1:212388/1‑134 (MQ=255)                                                                             cTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCaa                                                                                       <  2:212388/134‑1 (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                                  >  3:361464/1‑149 (MQ=255)                                                                                                                   cTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAgg                                  >  8:145882/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:218571/149‑1 (MQ=255)                                                                                                                                                  cATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGtt   >  4:350366/1‑149 (MQ=255)                                                                                                                                                    tCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGTTg  >  2:3017/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GGTGCTTAGCGCCGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGTTG  >  NC_002947/3549772‑3550066

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