Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,879,748	T→C	6.6%	F528S (TTC→TCC)	mexF →	multidrug RND transporter MexF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,879,748	0	T	C	6.6%	174.1 / 2.4	61	F528S (TTC→TCC)	mexF	multidrug RND transporter MexF
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base C (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

ACCGTCATTTCGGCGTTCAACTCGCTGACCCTGTCGCCAGCCTTGGCCGCGGTGCTGCTCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAG  >  NC_002947/3879603‑3879893                                                                                                                                                  |                                                                                                                                                   aCCGTCATTTCGGCGTTCAACTCGCTGACCCTGTCGCCAGCCTTGGCCGCGGTGCTGCTCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGttctt                                                                                                                                                 <  1:146165/149‑1 (MQ=255)  ccGTCATTTCGGCGTTCAACTCGCTGACCCTGTCGCCAGCCTTGGCCGCGGTGCTGCTCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGttcttc                                                                                                                                                <  1:241864/149‑1 (MQ=255)               gTTCAACTCGCTGACCCTGTCGCCAGCCTTGGCCGCGGTGCTGCTCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGgccagc                                                                                                                                    >  7:312324/1‑148 (MQ=255)                   aaCTCGCTGACCCTGTCGCCAGCCTTGGCCGCGGTGCTGCTCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCgg                                                                                                                                          <  7:218946/138‑1 (MQ=255)                   aaCTCGCTGACCCTGTCGCCAGCCTTGGCCGCGGTGCTGCTCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCgg                                                                                                                                          >  8:218946/1‑138 (MQ=255)                   aaCTCGCTGACCCTGTCGCCAGCCTTGGCCGCGGTGCTGCTCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCcaca                                                                                                                                <  4:352488/148‑1 (MQ=255)                   aaCTCGCTGACCCTGTCGCCAGCCTTGGCCGCGGTGCTGCTCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCcaca                                                                                                                                <  2:117073/148‑1 (MQ=255)                            aCCCTGTCGCCAGCCTTGGCCGCGGTGCTGCTCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGttcttc                                                                                                                                                <  4:342906/123‑1 (MQ=255)                            aCCCTGTCGCCAGCCTTGGCCGCGGTGCTGCTCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGttcttc                                                                                                                                                >  3:342906/1‑123 (MQ=255)                             cccTGTCGCCAGCCTTGGCCGCGGTGCTGCTCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGt                                                                                                                     <  2:82461/149‑1 (MQ=255)                              ccTGTCGCCAGCCTTGGCCGCGGTGCTGCTCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTg                                                                                                                    >  1:82506/1‑149 (MQ=255)                                          cTTGGCCGCGGTGCTGCTCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGTTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCggg                                                                                                        >  8:274254/1‑149 (MQ=255)                                                 gcgGTGCTGCTCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCg                                                                                                 >  5:50512/1‑149 (MQ=255)                                                  cgGTGCTGCTCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCg                                                                                                 >  5:323965/1‑148 (MQ=255)                                                        tgctCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGc                                                                                         >  8:248037/1‑149 (MQ=255)                                                         gctCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCg                                                                                        >  4:290590/1‑149 (MQ=255)                                                          ctCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCggg                                                                                                                                         <  2:322728/100‑1 (MQ=255)                                                          ctCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCggg                                                                                                                                         >  1:322728/1‑100 (MQ=255)                                                          ctCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCgg                                                                                       >  6:20936/1‑149 (MQ=255)                                                                gATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGc                                                                                 >  3:104339/1‑149 (MQ=255)                                                                  tcatcaCGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGccc                                                                               >  3:296422/1‑149 (MQ=255)                                                                       aCGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTc                                                                          >  3:255211/1‑149 (MQ=255)                                                                       aCGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTc                                                                          <  4:104339/149‑1 (MQ=255)                                                                       aCGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTc                                                                          <  6:354136/149‑1 (MQ=255)                                                                       aCGCACCCAAGGACCGTTTCTCGCGGTTCCACGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTc                                                                          >  1:29028/1‑149 (MQ=255)                                                                             ccAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTAt                                                                    >  8:28053/1‑149 (MQ=255)                                                                             ccAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTAt                                                                    >  2:201124/1‑149 (MQ=255)                                                                               aaGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGc                                                                  >  4:235784/1‑149 (MQ=255)                                                                               aaGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACACCTACGTCGGTGGCGTTCGCCCGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTAttc                                                                  >  4:261350/1‑147 (MQ=255)                                                                                 ggACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCg                                                                <  7:28053/149‑1 (MQ=255)                                                                                   aCCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGc                                                              <  5:20936/149‑1 (MQ=255)                                                                                       tttCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCt                                                                                                                             <  6:19857/83‑1 (MQ=255)                                                                                       tttCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCt                                                                                                                             >  5:19857/1‑83 (MQ=255)                                                                                         tgtcGCGGTTCCTCGACAAGCTGCTTGGCAGCTGGCTGTTCTCCCCTTTCCACCGGTGCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGAt                                                         <  3:261350/146‑1 (MQ=255)                                                                                             gcgGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGt                                                     <  7:274254/148‑1 (MQ=255)                                                                                              cgGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTc                                                    >  1:133428/1‑148 (MQ=255)                                                                                                  tCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGacc                                               >  5:180619/1‑149 (MQ=255)                                                                                                    cTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGacct                                              <  2:82506/148‑1 (MQ=255)                                                                                                    cTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATCCCGGCCTGATGGGTCTGaccta                                             >  6:258183/1‑149 (MQ=255)                                                                                                           aGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGc                                      <  6:323965/149‑1 (MQ=255)                                                                                                            gctgctGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCt                                     >  5:36427/1‑149 (MQ=255)                                                                                                                ctgGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTtctc                                 <  8:312324/149‑1 (MQ=255)                                                                                                                ctgGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTtctc                                 <  2:29028/149‑1 (MQ=255)                                                                                                                ctgGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTtctc                                 <  6:180619/149‑1 (MQ=255)                                                                                                                 tgGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCACCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATgg                                                       <  6:75772/126‑1 (MQ=255)                                                                                                                   ggCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGt                               <  6:36427/148‑1 (MQ=255)                                                                                                                        ctggctgTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTtctc                                 >  5:54156/1‑141 (MQ=255)                                                                                                                        ctggctgTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTtctc                                 <  6:54156/141‑1 (MQ=255)                                                                                                                          ggctgTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTa                                                                     <  6:331744/103‑1 (MQ=255)                                                                                                                          ggctgTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTa                                                                     >  5:331744/1‑103 (MQ=255)                                                                                                                               ttCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCg                   <  4:296422/148‑1 (MQ=255)                                                                                                                               ttCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCg                   <  1:201124/148‑1 (MQ=255)                                                                                                                                  tcCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTc                                                                          <  2:38897/90‑1 (MQ=255)                                                                                                                                  tcCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTc                                                                          >  1:38897/1‑90 (MQ=255)                                                                                                                                  tcCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTc                                         <  2:213874/123‑1 (MQ=255)                                                                                                                                  tcCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTc                                         >  1:213874/1‑123 (MQ=255)                                                                                                                                    cccTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTc              <  7:248037/148‑1 (MQ=255)                                                                                                                                    cccTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCACTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGtt                <  6:50512/146‑1 (MQ=255)                                                                                                                                    cccTTTCAACCGGTCCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGAcc                    >  6:33329/1‑142 (MQ=255)                                                                                                                                    cccTTTCAACCGGTCCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGAcc                    <  5:27849/142‑1 (MQ=255)                                                                                                                                    cccTTTCAACCGGTCCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGAcc                    >  6:27849/1‑142 (MQ=255)                                                                                                                                    cccTTTCAACCGGTCCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGAcc                    <  5:33329/142‑1 (MQ=255)                                                                                                                                     ccTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGt            >  4:241480/1‑149 (MQ=255)                                                                                                                                             ccGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCCAATCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGcc     >  1:111255/1‑149 (MQ=255)                                                                                                                                               ggTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAg  <  2:133428/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   ACCGTCATTTCGGCGTTCAACTCGCTGACCCTGTCGCCAGCCTTGGCCGCGGTGCTGCTCAAGGATCATCACGCACCCAAGGACCGTTTCTCGCGGTTCCTCGACAAGCTGCTGGGCAGCTGGCTGTTCTCCCCTTTCAACCGGTTCTTCGACCGGGCCAGCCACAGCTACGTCGGTGGCGTGCGCCGGGTCATCCGCTCC‑AGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAG  >  NC_002947/3879603‑3879893

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