Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,212,988	T→A	10.0%	L19M (TTG→ATG)	mexC →	multidrug RND transporter membrane fusion protein MexC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,212,988	0	T	A	10.0%	187.2 / 11.4	70	L19M (TTG→ATG)	mexC	multidrug RND transporter membrane fusion protein MexC
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base A (5/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

TAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGAGGTGC  >  NC_002947/3212843‑3213132                                                                                                                                                  |                                                                                                                                                  tAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTgg                                                                                                                                                >  3:309187/1‑149 (MQ=255) tAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGGGTGGGATGGGCAAGACAGGTTCGAGAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTgg                                                                                                                                                >  8:115249/1‑149 (MQ=255)  aaGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGTTAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTgg                                                                                                                                                <  1:259376/148‑1 (MQ=255)   agagTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGc                                                                                                                                               <  1:61723/148‑1 (MQ=255)       tCTTTTATTATTTAAATGAGACAATTTTGACTCTTGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCTATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGgt                                                                                                                                          <  4:232912/149‑1 (MQ=255)       tCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGGTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGgt                                                                                                                                          <  7:206562/149‑1 (MQ=255)         ttttATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGgtgt                                                                                                                                        >  5:344846/1‑149 (MQ=255)         ttttATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGgtgt                                                                                                                                        >  4:45187/1‑149 (MQ=255)            tctGATTTTAATGCCACAATCTAGACTCATGGGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGCCCATGGCGATCGCCTTGGCCGGGTGTGa                                                                                                                                      <  4:309187/146‑1 (MQ=255)            tATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGa                                                                                                                                      <  5:40784/148‑1 (MQ=255)                aTTTTAATGAGACACTATTGACTCATGTGGGTTTGAAAATGAGCATCTGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGATCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGccc                                                                                                                                  <  4:350578/148‑1 (MQ=255)                 ttttAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCgc                                                                                                                                <  1:72632/149‑1 (MQ=255)                 ttttAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCg                                                                                                                                 <  4:295661/148‑1 (MQ=255)                 ttttAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCg                                                                                                                                 >  2:86445/1‑148 (MQ=255)                 ttttAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCg                                                                                                                                 >  3:256225/1‑148 (MQ=255)                 ttttAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGAGAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCg                                                                                                                                 >  4:113565/1‑148 (MQ=255)                  tttAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCg                                                                                                                                 >  3:349620/1‑147 (MQ=255)                        gagaCAATATTGACTCATGTCGGTTTGAAAAGTAGCATCCGCCCCGATTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTCTGAGCCCtcagca                                                                                                                            >  6:307245/1‑146 (MQ=255)                         agaCAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagag                                                                                                                         >  8:12250/1‑148 (MQ=255)                          gaCAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCACTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagc                                                                                                                        <  3:102800/148‑1 (MQ=255)                               tatTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACg                                                                                                                  >  4:18802/1‑149 (MQ=255)                                        aTGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑c                                                                                                          <  4:349620/147‑1 (MQ=255)                                              gTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGc                                                                                                    >  3:121262/1‑147 (MQ=255)                                              gTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGa                                                                                                  <  3:113565/149‑1 (MQ=255)                                              gTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGa                                                                                                  <  1:363384/149‑1 (MQ=255)                                               tttGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAg                                                                                                 >  6:331454/1‑149 (MQ=255)                                                   aaaaGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGGGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTAc                                                                                             <  5:76010/149‑1 (MQ=255)                                                        tAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTAcc                                                                                            >  3:56606/1‑145 (MQ=255)                                                        tAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTAcc                                                                                            <  4:56606/145‑1 (MQ=255)                                                        tAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGt                                                                                        >  4:22993/1‑149 (MQ=255)                                                        tAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGt                                                                                        >  6:218823/1‑149 (MQ=255)                                                            aTCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagcag                                                                                                                      <  5:101294/116‑1 (MQ=255)                                                            aTCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagcag                                                                                                                      >  6:101294/1‑116 (MQ=255)                                                                gCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCg                                                                                >  5:287591/1‑149 (MQ=255)                                                                   ccGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGa                                                                             <  4:121262/149‑1 (MQ=255)                                                                     gTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCGTGGCGATGGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGGTACGCGGTGGAAATCGTGACa                                                                           >  8:60718/1‑149 (MQ=255)                                                                         tcaatcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTg                                                                        <  6:344846/148‑1 (MQ=255)                                                                          caatcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGTTAAGGGTGGGCTCATTGGCCATGGCGATCGCCTTGGCAGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTgg                                                                       <  5:307245/148‑1 (MQ=255)                                                                            atcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGgcgc                                                                    <  5:331454/149‑1 (MQ=255)                                                                            atcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGgcgc                                                                    >  1:390223/1‑149 (MQ=255)                                                                               aaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGccc                                                                  >  8:226926/1‑148 (MQ=255)                                                                                aGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACa                                                                           <  7:198324/138‑1 (MQ=255)                                                                                aGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACa                                                                           >  8:198324/1‑138 (MQ=255)                                                                                aGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGccc                                                                  >  6:282242/1‑147 (MQ=255)                                                                                       ggTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTc                                                         >  5:339952/1‑149 (MQ=255)                                                                                          gtgtgATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGCGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGcc                                                      <  3:297303/149‑1 (MQ=255)                                                                                              gATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTgg                                                  >  1:287559/1‑149 (MQ=255)                                                                                               aTGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTggc                                                 <  7:12250/149‑1 (MQ=255)                                                                                               aTGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTggc                                                 <  3:307424/149‑1 (MQ=255)                                                                                               aTGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTggc                                                 >  5:144473/1‑149 (MQ=255)                                                                                                  ggCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTggcggc                                              <  4:256225/149‑1 (MQ=255)                                                                                                   gCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCt                                             <  2:390223/149‑1 (MQ=255)                                                                                                   gCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCt                                             <  3:18802/149‑1 (MQ=255)                                                                                                          aGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgc                                      >  1:78592/1‑149 (MQ=255)                                                                                                            gTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgccg                                    <  6:287591/149‑1 (MQ=255)                                                                                                                gATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgccgggc                                 <  3:22993/148‑1 (MQ=255)                                                                                                                 aTAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgccgggccg                               >  2:53663/1‑149 (MQ=255)                                                                                                                              acccggCCAAGGCGATCGCCATGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATg                  >  2:376108/5‑149 (MQ=255)                                                                                                                               aTTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCACAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGCCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATgc                 <  6:64862/149‑1 (MQ=255)                                                                                                                                ccggCCAAGGCGATCGCCATGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATgcg                >  8:346626/3‑149 (MQ=255)                                                                                                                                ccggCCAAGGCGATCGCCATGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATgcg                >  7:52211/3‑149 (MQ=255)                                                                                                                                  ggCCAAGGCGATCGCCATGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGccgggccggg                             >  4:237856/1‑134 (MQ=255)                                                                                                                                  ggCCAAGGCGATCGCCATGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGccgggccggg                             <  3:237856/134‑1 (MQ=255)                                                                                                                                   gCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGt             <  6:339952/149‑1 (MQ=255)                                                                                                                                   gCCAAGGCGATCGCCATGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGATGAACTGCCGGGCCGGGTCGAGCCGATGCGCGt             >  4:204555/1‑149 (MQ=255)                                                                                                                                    ccATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGgc                                                                      >  2:265551/1‑91 (MQ=255)                                                                                                                                    ccATGGCGATCGCCTTGGCCGGGTGTGACCCCGCATCAGAGCAGGACGGTGCAGACCAATGCGAGCTACCCGGTGGAAATCGTGACATTGgc                                                                      <  1:265551/92‑1 (MQ=255)                                                                                                                                         gCGATCGCCATGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGa       <  7:346626/149‑1 (MQ=255)                                                                                                                                           gATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGAgg     >  6:150677/1‑149 (MQ=255)                                                                                                                                               gCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGAGGTGc  <  7:226926/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  TAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAG‑CCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGAGGTGC  >  NC_002947/3212843‑3213132

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