Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,297,144	T→A	6.0%	D993V (GAC→GTC) ‡	sbcC ←	exonuclease subunit SbcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,297,144	0	T	A	6.0%	241.7 / 4.1	83	D993V (GAC→GTC) ‡	sbcC	exonuclease subunit SbcC
Reads supporting (aligned to +/- strand):  ref base T (36/42);  new base A (3/2);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

TTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCATGGCCAGCAGGCGATCGAG  >  NC_002947/2296996‑2297292                                                                                                                                                     |                                                                                                                                                     ttCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGgtgt                                                                                                                                                      >  4:139006/1‑149 (MQ=255)        tCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACgg                                                                                                                                                <  2:303037/148‑1 (MQ=255)          tGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCTGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCGACGGTCa                                                                                                                                             >  4:24775/1‑149 (MQ=255)            ccgccAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATc                                                                                                                                           <  3:95770/149‑1 (MQ=255)            ccgccAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATc                                                                                                                                           <  2:151449/149‑1 (MQ=255)            ccgccAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGACACCACGGTcc                                                                                                                                             <  5:207551/147‑2 (MQ=255)            ccgccAGCGCCTGGTGGGCCTGCGGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGACACCACGGTccaa                                                                                                                                           <  8:303453/149‑4 (MQ=255)              gccAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCtg                                                                                                                                         >  5:309607/1‑149 (MQ=255)                    gcCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCcgc                                                                                                                                   >  3:1134/1‑149 (MQ=255)                    gcCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCcgc                                                                                                                                   <  7:261230/149‑1 (MQ=255)                    gcCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCcgc                                                                                                                                   <  2:154965/149‑1 (MQ=255)                      cTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTATGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGt                                                                                                                                 <  1:104456/149‑1 (MQ=255)                         gtggGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGAt                                                                                                                              <  2:317000/149‑1 (MQ=255)                           gggCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGAtgc                                                                                                                            >  1:298429/1‑149 (MQ=255)                            ggCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGAGTCCACGGTCATCTGTGCCGCGTGAtgc                                                                                                                            <  6:146680/148‑1 (MQ=255)                             gCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGAtgctg                                                                                                                          >  3:25614/1‑149 (MQ=255)                              ccTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCATGGTCATCTGTGCCGCGTGAtgctgc                                                                                                                         <  4:321751/149‑1 (MQ=255)                              ccTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGAtgctg                                                                                                                          >  8:293799/1‑148 (MQ=255)                              ccTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGAtgctg                                                                                                                          >  3:118712/1‑148 (MQ=255)                                  ctgGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAg                                                                                                                       <  5:53330/147‑1 (MQ=255)                                   tgGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGgcg                                                                                                                    >  4:25940/1‑149 (MQ=255)                                        gcgcTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGtt                                                                                                               <  3:139006/149‑1 (MQ=255)                                        gcgcTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGt                                                                                                                >  7:286187/1‑148 (MQ=255)                                         cgcTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTc                                                                                                              <  4:1134/149‑1 (MQ=255)                                              gtcgtcGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCgcgc                                                                                                          <  3:177191/148‑1 (MQ=255)                                              gtcgtcGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCgcgc                                                                                                          <  5:159239/148‑1 (MQ=255)                                                    ggCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTgc                                                                                                    <  3:25940/148‑1 (MQ=255)                                                        tgttgtGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTcg                                                                                                             >  4:85705/1‑135 (MQ=255)                                                        tgttgtGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTcg                                                                                                             <  3:85705/135‑1 (MQ=255)                                                         gttgtGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGcagcag                                                                                               >  8:227416/1‑148 (MQ=255)                                                          ttgtGCACGCAATTCGGCACATTGCTGCTCGTGGTTTGCCGGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGcagcag                                                                                               <  8:252037/147‑1 (MQ=255)                                                             tgCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGAcgc                                                                                          >  6:43620/1‑149 (MQ=255)                                                               cACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCgg                                                                                        <  3:91724/149‑1 (MQ=255)                                                               cACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCgg                                                                                        <  1:175356/149‑1 (MQ=255)                                                                aCGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGGCCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGc                                                                                       >  3:300155/1‑149 (MQ=255)                                                                aCGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGc                                                                                       <  6:153215/149‑1 (MQ=255)                                                                 cGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGcc                                                                                      >  8:135078/1‑149 (MQ=255)                                                                      ttCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCt                                                                                 >  4:29742/1‑149 (MQ=255)                                                                      ttCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCt                                                                                 >  3:207150/1‑149 (MQ=255)                                                                           cgcATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATg                                                                            <  7:293799/149‑1 (MQ=255)                                                                           cgcATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCACCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATg                                                                            <  1:269932/149‑1 (MQ=255)                                                                              aTTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCt                                                                         >  1:84737/1‑149 (MQ=255)                                                                                  ctgctCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGATGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCttttt                                                                     >  8:204292/1‑149 (MQ=255)                                                                                     ctCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTcc                                                                   <  7:227416/148‑1 (MQ=255)                                                                                      tCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCcgc                                                                 >  8:4292/1‑149 (MQ=255)                                                                                      tCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCcgc                                                                 >  3:194241/1‑149 (MQ=255)                                                                                         tGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGcc                                                              >  3:170111/1‑149 (MQ=255)                                                                                           gTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATTCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCtt                                                            <  5:43620/149‑1 (MQ=255)                                                                                                ccAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGAGCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCaaa                                                       >  1:231838/1‑149 (MQ=255)                                                                                                ccAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTTATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATCGCTTTTTCCGCGCCTTGCaaa                                                       <  4:25614/149‑1 (MQ=255)                                                                                                 cAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAAc                                                      >  1:96338/1‑149 (MQ=255)                                                                                                 cAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAAc                                                      >  2:112731/1‑149 (MQ=255)                                                                                                     cGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTg                                                  <  8:286187/149‑1 (MQ=255)                                                                                                      gTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTgg                                                 <  4:170111/149‑1 (MQ=255)                                                                                                      gTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTgg                                                 >  8:134564/1‑149 (MQ=255)                                                                                                       ttCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGgc                                                <  4:207150/149‑1 (MQ=255)                                                                                                            gcAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAAt                                           >  4:301934/1‑149 (MQ=255)                                                                                                            gcAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAAt                                           >  1:245897/1‑149 (MQ=255)                                                                                                             cAGGTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAAtt                                          <  6:303198/149‑1 (MQ=255)                                                                                                              aGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTc                                         <  2:84737/149‑1 (MQ=255)                                                                                                               gTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCa                                        <  1:112731/149‑1 (MQ=255)                                                                                                               gTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCa                                        <  4:118712/149‑1 (MQ=255)                                                                                                               gTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCa                                        <  5:320304/149‑1 (MQ=255)                                                                                                                  cGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCAtt                                      <  7:134564/148‑1 (MQ=255)                                                                                                                  cGGCCAAGCCTTGCTCCAGGGCTTGGACCGGGGTGTCCACGGTCAGCTGTGCCGCGTGATGCTGCACGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCAtt                                      <  4:57591/148‑1 (MQ=255)                                                                                                                   ggCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGa                                    <  5:329411/149‑1 (MQ=255)                                                                                                                    gCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGAc                                   >  8:126100/1‑149 (MQ=255)                                                                                                                     ccAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACt                                  >  2:151776/1‑149 (MQ=255)                                                                                                                      cAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACtt                                 >  7:251927/1‑149 (MQ=255)                                                                                                                            ttGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGcgcg                            >  6:300551/1‑148 (MQ=255)                                                                                                                               cTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCa                         <  2:96338/148‑1 (MQ=255)                                                                                                                               cTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCAt                        <  2:298429/149‑1 (MQ=255)                                                                                                                               cTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAAGTCAATGAATTgcgccagc                         <  4:194241/148‑5 (MQ=255)                                                                                                                                tCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATc                       >  4:40054/1‑149 (MQ=255)                                                                                                                                tCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATc                       >  7:21020/1‑149 (MQ=255)                                                                                                                                  cAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCa                     >  6:13564/1‑149 (MQ=255)                                                                                                                                  cAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCa                     >  2:165383/1‑149 (MQ=255)                                                                                                                                   aGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCAt                    >  8:11813/1‑149 (MQ=255)                                                                                                                                   aGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCAt                    <  3:40054/149‑1 (MQ=255)                                                                                                                                    gggCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCATg                   <  2:231838/149‑1 (MQ=255)                                                                                                                                        gatgACCGTGGACACCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCATGGc                 >  6:105123/4‑147 (MQ=255)                                                                                                                                           gACCGTGGACACCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCATGGCcag              >  7:267185/1‑147 (MQ=255)                                                                                                                                           gACCGTGGACACCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCATGGCcag              >  2:100818/1‑147 (MQ=255)                                                                                                                                             ccGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCATGGCCAGCAgg          <  1:134937/149‑1 (MQ=255)                                                                                                                                                  gtgtCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCATGGCCAGCAGGCGATc     >  4:227638/1‑149 (MQ=255)                                                                                                                                                     tCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCATGGCCAGCAGGCGATCGa   <  6:131491/148‑1 (MQ=255)                                                                                                                                                     tCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCATGGCCAGCAGGCGATCGAg  <  6:75343/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCATGGCCAGCAGGCGATCGAG  >  NC_002947/2296996‑2297292

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