Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,883,615	C→G	7.7%	A197G (GCC→GGC) ‡	PP_3428 →	TPR repeat‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,883,615	0	C	G	7.7%	240.6 / 9.9	78	A197G (GCC→GGC) ‡	PP_3428	TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (30/42);  new base G (3/3);  total (33/45)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

CATCGGCCTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGACTCGG  >  NC_002947/3883472‑3883752                                                                                                                                                 |                                                                                                                                          gaTCGGCCGCTATTTCGGACCACTGCAGTCAGTCCATGATTAACGCCATCAAGCCCGGGGGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGAATACACGCACGC‑CGCTGCAGCCCTGCa                                                                                                                                     <  2:959628/149‑1 (MQ=255) cATCGGCGTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGGAAGCCCTGGAGCAACAGCCTGAAGATGCCGGCGCCCAAGTGCACCTGGCGCGGCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCa                                                                                                                                     <  2:979131/150‑1 (MQ=255) cATCGGCCTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCa                                                                                                                                     <  2:401576/150‑1 (MQ=255)  aTCGGCCTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCag                                                                                                                                    >  1:868630/1‑150 (MQ=255)     ggCCTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCagcga                                                                                                                                 <  1:1677434/150‑1 (MQ=255)     ggCCTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCagcga                                                                                                                                 <  2:1067196/150‑1 (MQ=255)     ggCCTCTATTTCGGCCCCCTGCAATGTTATCTAGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCagcga                                                                                                                                 <  1:1433865/150‑1 (MQ=255)       tctctATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCTAGTCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCTCGCTGCAGGGCTGCagcggc                                                                                                                                <  1:1003773/149‑3 (MQ=255)        ctctATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCagcgagcg                                                                                                                              >  2:1374870/1‑150 (MQ=255)        ctctATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCagcgagcg                                                                                                                              >  2:948363/1‑150 (MQ=255)        ctctATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCagcgagcg                                                                                                                              <  1:1111740/150‑1 (MQ=255)           tATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCTAGTCCTGGAGCAACATCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACGTGCACGAAAACGACTACACGCACGCTCGCTGAAGGGCTACagcggcgaga                                                                                                                            >  2:1003773/1‑144 (MQ=255)             tttCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTcgc                                                                                                                         >  2:1305553/1‑150 (MQ=255)                      ccTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGcccc                                                                                                                <  1:1305553/150‑1 (MQ=255)                         gCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCgg                                                                                                             >  2:1912300/1‑150 (MQ=255)                                 cccTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCagcga                                                                                                                                 <  2:648585/122‑1 (MQ=255)                                 cccTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCagcga                                                                                                                                 >  1:648585/1‑122 (MQ=255)                                  ccTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGGGCTGCagcg                                                                                                                                  <  2:1092638/120‑1 (MQ=255)                                  ccTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGGGCTGCagcg                                                                                                                                  >  1:1092638/1‑120 (MQ=255)                                          aaTACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACTTGGCGCGCCTGTACCTGCACGAACACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTgg                                                                                            <  2:1035837/150‑1 (MQ=255)                                          aaTACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTgg                                                                                            <  1:217701/150‑1 (MQ=255)                                           aTACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGc                                                                                           <  1:1900744/150‑1 (MQ=255)                                           aTACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGc                                                                                           <  1:1900801/150‑1 (MQ=255)                                            tACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTg                                                                                             >  2:461365/1‑147 (MQ=255)                                            tACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTg                                                                                             <  1:461365/147‑1 (MQ=255)                                             aCCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGa                                                                                         <  1:393615/150‑1 (MQ=255)                                              ccAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGAc                                                                                        >  1:533650/1‑150 (MQ=255)                                                  cAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAg                                                                                    >  1:839465/1‑150 (MQ=255)                                                   aaGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAgg                                                                                   <  2:1465099/150‑1 (MQ=255)                                                        cTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGc                                                                              <  2:1603436/150‑1 (MQ=255)                                                              cAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCCATACCTGCACGAAAACGACTACATGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTc                                                                        <  1:206009/150‑1 (MQ=255)                                                                   gCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCGGGGGCAACCGCGa                                                                                                    <  2:747393/117‑1 (MQ=255)                                                                      tGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTa                                                                                 >  1:138150/1‑133 (MQ=255)                                                                      tGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTa                                                                                 <  2:138150/133‑1 (MQ=255)                                                                      tGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCa                                                                >  1:1621375/1‑150 (MQ=255)                                                                      tGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCa                                                                >  1:1844063/1‑150 (MQ=255)                                                                        aaGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAgg                                                              >  2:297701/1‑150 (MQ=255)                                                                        aaGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAgg                                                              <  2:1621375/150‑1 (MQ=255)                                                                        aaGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAgg                                                              >  2:1206926/1‑150 (MQ=255)                                                                          gATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGc                                                            <  2:868630/150‑1 (MQ=255)                                                                           aTGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGcc                                                           >  1:1908655/1‑150 (MQ=255)                                                                               cGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCgg                                                                                                                       <  2:1544313/86‑1 (MQ=255)                                                                               cGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCgg                                                                                                                       >  1:1544313/1‑86 (MQ=255)                                                                               cGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAg                                                                                    >  1:388796/1‑121 (MQ=255)                                                                               cGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAg                                                                                    <  2:388796/121‑1 (MQ=255)                                                                                gACGCTCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTc                                                                                                      <  2:59786/102‑1 (MQ=255)                                                                                gACGCTCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTc                                                                                                      >  1:59786/1‑102 (MQ=255)                                                                                   gCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCacgac                                                   <  2:1844063/150‑1 (MQ=255)                                                                                   gCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCacgac                                                   <  1:530534/150‑1 (MQ=255)                                                                                     ccAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGc                                                 >  2:1286227/1‑150 (MQ=255)                                                                                     ccAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGc                                                 >  2:599172/1‑150 (MQ=255)                                                                                     ccAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGc                                                 >  2:184277/1‑150 (MQ=255)                                                                                      cAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGcc                                                >  2:1111165/1‑150 (MQ=255)                                                                                      cAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGcc                                                >  1:1447404/1‑150 (MQ=255)                                                                                      cAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGcc                                                >  1:1282280/1‑150 (MQ=255)                                                                                            cACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGc                                          <  1:1912300/150‑1 (MQ=255)                                                                                                      gcCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTgg                                <  1:1111165/150‑1 (MQ=255)                                                                                                      gcCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTgg                                <  1:948363/150‑1 (MQ=255)                                                                                                           tACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCt                                  <  2:178051/143‑1 (MQ=255)                                                                                                           tACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCt                                  >  1:178051/1‑143 (MQ=255)                                                                                                            cccTGCACGAAAACGACTACACCCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCTCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCGCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAg                          <  2:649583/149‑1 (MQ=255)                                                                                                             ccTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAACGAGCGCTCTCCGTTGTCCTCGGCATCCCCGTAGCCCCGGCGCCTCAGGTAGGCCTTCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGc                         <  1:890364/150‑1 (MQ=255)                                                                                                              cTGCACGAAAACGTCTACCCGCTCGT‑TTCTGCATCTTTGTAGCTATCGTTCTCTGGTGCTCCCGGTATCCTTGATGCCCTGGCGACTCAGGTACGCCTTCTCTAGCGCCATGGTCACGTCGATGTCTCCCGCCATGTGCTGGGCAAATcc                        <  2:1342727/150‑1 (MQ=255)                                                                                                               tGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCTTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGccc                       >  2:820256/1‑150 (MQ=255)                                                                                                               tGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGccc                       <  1:1877460/150‑1 (MQ=255)                                                                                                                      aaaCGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCt                <  2:1709466/150‑1 (MQ=255)                                                                                                                            cTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGccc          <  2:1282280/150‑1 (MQ=255)                                                                                                                            cTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGccc          <  2:839465/150‑1 (MQ=255)                                                                                                                            cTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGccc          <  2:1460297/150‑1 (MQ=255)                                                                                                                            cTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGccc          <  1:184277/150‑1 (MQ=255)                                                                                                                            cTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGccc          >  1:1137253/1‑150 (MQ=255)                                                                                                                                acgcacgc‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGa                                                                      <  1:916400/86‑1 (MQ=255)                                                                                                                                acgcacgc‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGa                                                                      >  2:916400/1‑86 (MQ=255)                                                                                                                                acgcacgc‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGAc      >  2:882505/1‑150 (MQ=255)                                                                                                                                    acgc‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGACTCgg  <  1:1286227/150‑1 (MQ=255)                                                                                                                                    acgc‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGACTCgg  >  1:150339/1‑150 (MQ=255)                                                                                                                                        tcgcTGCAGGGCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTgg                                >  1:1558217/2‑117 (MQ=255)                                                                                                                                        tcgcTGCAGGGCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTgg                                <  2:1558217/116‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                          CATCGGCCTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGC‑CGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGACTCGG  >  NC_002947/3883472‑3883752

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑