Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,066,366	G→C	5.4%	V153L (GTG→CTG)	PP_2678 →	hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,066,366	0	G	C	5.4%	300.8 / 6.0	93	V153L (GTG→CTG)	PP_2678	hydrolase
Reads supporting (aligned to +/- strand):  ref base G (45/43);  new base C (2/3);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

CCAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCT  >  NC_002947/3066224‑3066511                                                                                                                                               |                                                                                                                                                  ccAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCa                                                                                                                                            >  2:716962/1‑150 (MQ=255)  cAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAg                                                                                                                                           >  1:2353401/1‑150 (MQ=255)   aGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGc                                                                                                                                          >  1:756745/1‑150 (MQ=255)   aGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGc                                                                                                                                          >  2:594370/1‑150 (MQ=255)   aGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGc                                                                                                                                          >  2:2073264/1‑150 (MQ=255)    ggCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGcc                                                                                                                                         >  1:607784/1‑150 (MQ=255)    ggCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGcc                                                                                                                                         <  2:1319052/150‑1 (MQ=255)        tttGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGc                                                                                                                                     <  2:787066/150‑1 (MQ=255)        tttGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGc                                                                                                                                     <  1:925382/150‑1 (MQ=255)            ccGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGc                                                                                                                                 <  1:232114/150‑1 (MQ=255)              gATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCACGCCCGGCTGCAGCACCCgggtcggcag                                                                                                                               <  2:454218/150‑8 (MQ=255)                 gtgCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAg                                                                                                                            <  1:780560/150‑1 (MQ=255)                  tgCCGATCGGTGCCCTGGCCGGCACCCGCGACATGCTGCGCGAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAgg                                                                                                                           <  2:529019/150‑1 (MQ=255)                      gATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtggg                                                                                                                       >  1:342908/1‑150 (MQ=255)                        tCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggtg                                                                                                                     <  2:607784/150‑1 (MQ=255)                         cGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGGTGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGgtgggtgg                                                                                                                    <  1:1983073/150‑1 (MQ=255)                         cGGTGCCCTGGCCGGCACCCGCCACCGGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGgtgggtgg                                                                                                                    <  1:733086/150‑1 (MQ=255)                           gTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGCTGCGCGGTACCGAGGTGGTGCTGCCGACCCGGCTGCAGCACCCgg                                                                                                                                       <  2:2203150/129‑1 (MQ=255)                           gTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGCTGCAGCACCCgg                                                                                                                                       >  1:2203150/1‑129 (MQ=255)                               ccTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCg                                                                                                              >  2:2309873/1‑150 (MQ=255)                                 tGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggt                                                                                                                      >  2:1678501/1‑140 (MQ=255)                                 tGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggt                                                                                                                      <  1:1678501/140‑1 (MQ=255)                                 tGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAg                                                                                                                            >  2:139797/1‑134 (MQ=255)                                 tGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAg                                                                                                                            <  1:139797/134‑1 (MQ=255)                                   gCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGCTGCAGCACCCgggtc                                                                                                                                    >  2:593231/1‑122 (MQ=255)                                   gCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGCTGCAGCACCCgggtc                                                                                                                                    <  1:593231/124‑3 (MQ=255)                                     aagCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAg                                                                                                        >  2:1902570/3‑150 (MQ=255)                                        cACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATTCTCGGCACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGGGGTGGGTGGGCACCTGCTGAGGtt                                                                                                     >  1:1489512/1‑150 (MQ=255)                                         aCCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTg                                                                                                    >  1:2389573/1‑150 (MQ=255)                                           ccGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGACGGGGACTGGATGCTCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTgct                                                                                                  >  1:1499851/1‑150 (MQ=255)                                            cgcgACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtg                                                                                                                         >  1:2075293/1‑126 (MQ=255)                                            cgcgACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtg                                                                                                                         <  2:2075293/126‑1 (MQ=255)                                             gcgACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTgctgg                                                                                                <  2:1155675/150‑1 (MQ=255)                                             gcgACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTgctgg                                                                                                <  2:952946/150‑1 (MQ=255)                                               gACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTgctgggg                                                                                              >  2:1115427/1‑150 (MQ=255)                                                  ctgctgCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTCGGTGGGCACCGCCTGAGGTTgctggggctg                                                                                           >  2:663418/1‑150 (MQ=255)                                                   tgctgcGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGGCTGGTTGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGGGGGCACCGCCTGAGGTTgctggggctgg                                                                                          >  1:1623594/1‑150 (MQ=255)                                                   tgctgcGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTgctggggctgg                                                                                          >  1:755742/1‑150 (MQ=255)                                                     ctgcGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTgggc                                                                                        <  1:2309873/150‑1 (MQ=255)                                                      tgcGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTgggcg                                                                                       <  2:1884751/150‑1 (MQ=255)                                                       gcgcCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTgggcgg                                                                                      <  1:716962/150‑1 (MQ=255)                                                          ccAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCa                                                                                   >  2:789985/1‑150 (MQ=255)                                                             gCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCAt                                                                                                                               <  2:2318363/103‑1 (MQ=255)                                                             gCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCAt                                                                                                                               >  1:2318363/1‑103 (MQ=255)                                                             gCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATAc                                                                                >  2:2185495/1‑150 (MQ=255)                                                                 ggCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGt                                                                            <  2:2389573/150‑1 (MQ=255)                                                                    gATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggtg                                                                                                                     <  1:1760891/106‑1 (MQ=255)                                                                    gATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggtg                                                                                                                     >  2:1760891/1‑106 (MQ=255)                                                                    gATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGcc                                                                         >  2:2228099/1‑150 (MQ=255)                                                                    gATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGAAGCTGCCGACCCGGGAGCTGCAGCCGGGCGTGAAGGAGGTGGGTGGGGACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATCGGGGTAcc                                                                         >  2:1166293/1‑150 (MQ=255)                                                                      tGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGa                                                                       <  2:1443213/150‑1 (MQ=255)                                                                        ccATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCa                                                                                                                                            >  2:53758/1‑79 (MQ=255)                                                                        ccATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCa                                                                                                                                            <  1:53758/79‑1 (MQ=255)                                                                          aTGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTg                                                                   >  1:315153/1‑150 (MQ=255)                                                                           tGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTgg                                                                  >  1:339742/1‑150 (MQ=255)                                                                            ggCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGc                                                                 <  2:315153/150‑1 (MQ=255)                                                                               cGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCa                                                                                                                                <  1:1666049/84‑1 (MQ=255)                                                                               cGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCa                                                                                                                                >  2:1666049/1‑84 (MQ=255)                                                                               cGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGAt                                                              >  2:1412448/1‑150 (MQ=255)                                                                                   aaCATGTATCGCCTGGTCGGCGACTGGATGTGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTc                                                          <  2:2198294/150‑1 (MQ=255)                                                                                      aTGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGAt                                                       <  1:1412448/150‑1 (MQ=255)                                                                                      aTGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGAt                                                       >  1:230515/1‑150 (MQ=255)                                                                                      aTGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGAt                                                       >  2:847321/1‑150 (MQ=255)                                                                                        gTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGa                                                     >  1:1751072/1‑150 (MQ=255)                                                                                          aTCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGaaa                                                   <  2:1903475/150‑1 (MQ=255)                                                                                           tCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggt                                                                                                                      >  2:1436933/1‑82 (MQ=255)                                                                                           tCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggt                                                                                                                      <  1:1436933/82‑1 (MQ=255)                                                                                           tCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTg                                                                                                          >  1:1014651/1‑94 (MQ=255)                                                                                           tCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTg                                                                                                          <  2:1014651/94‑1 (MQ=255)                                                                                             gCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGAGCTGGCGATAATCGAGTAAAAGa                                                >  2:1268737/1‑150 (MQ=255)                                                                                                tGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTTCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTCAGACGTGTCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCg                                             <  1:1268737/150‑1 (MQ=255)                                                                                                       cGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGAGCGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCt                                      <  2:342908/150‑1 (MQ=255)                                                                                                        gACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGa                                                                                                         >  2:1726434/1‑82 (MQ=255)                                                                                                        gACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGa                                                                                                         <  1:1726434/82‑1 (MQ=255)                                                                                                         aCTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGa                                                        <  1:650200/130‑1 (MQ=255)                                                                                                         aCTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGa                                                        >  2:650200/1‑130 (MQ=255)                                                                                                           tGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTc                                  <  2:1751072/150‑1 (MQ=255)                                                                                                              aTGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGcc                               <  1:594370/150‑1 (MQ=255)                                                                                                               tGCGCGGTACCGAGGTGGTGCTGGCGACCCGGGTGCTGCAGGCGGGCGTGCATGAGGTGGGTGGGCACCGCCGGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCg                              >  2:1853667/1‑150 (MQ=255)                                                                                                                gcgcGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCAt                                                                                                                               <  1:1732908/52‑1 (MQ=255)                                                                                                                gcgcGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCAt                                                                                                                               >  2:1732908/1‑52 (MQ=255)                                                                                                                       aCCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTg                      >  2:586242/1‑150 (MQ=255)                                                                                                                               ggtgCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGACGGGGCTGGGCGGGCTGACGGGGTCCGTCCTGGCGATTTTCTATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTa              <  2:2078868/150‑1 (MQ=255)                                                                                                                               ggtcCTGCCGACCCGGGTGCTGCAGGCGGGCGTGCATGTGTGGGGTGGGCACCGCCTGAGGTTGCTGGGGCTAGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTa              <  1:1737510/150‑1 (MQ=255)                                                                                                                                gtgCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCTCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTAc             <  1:1341134/150‑1 (MQ=255)                                                                                                                                gtgCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTAc             <  1:1115427/150‑1 (MQ=255)                                                                                                                                gtgCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTAc             <  1:586242/150‑1 (MQ=255)                                                                                                                                   ctgcCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGCCGGTTTTCGATGAAAAGACCGGGGTGCTTTGCGCCGGCGACCTCGTGCTCTACCaa          >  2:542769/1‑150 (MQ=255)                                                                                                                                    tgcCGACCCGGGTGCTGCAGGCGGGCGTGCATGACGTGGGTGCGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGAGCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAAc         <  2:1082686/150‑1 (MQ=255)                                                                                                                                    tgcCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAAc         >  2:2030359/1‑150 (MQ=255)                                                                                                                                           ccGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGCCTGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCt  >  1:1039658/1‑150 (MQ=255)                                                                                                                                           ccGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCt  >  2:2274310/1‑150 (MQ=255)                                                                                                                                           ccGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCt  <  1:2030359/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                  CCAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCT  >  NC_002947/3066224‑3066511

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