Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,743,598	G→A	5.3%	A48T (GCT→ACT)	mdtK →	multidrug efflux system transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,743,598	0	G	A	5.3%	238.6 / 2.2	76	A48T (GCT→ACT)	mdtK	multidrug efflux system transporter
Reads supporting (aligned to +/- strand):  ref base G (39/32);  new base A (2/2);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

CAGAAGTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTCTGGCTGGCAGGTTT  >  NC_000913/1743460‑1743746                                                                                                                                           |                                                                                                                                                     cAGAAGTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTa                                                                                                                                             <  5:222875/148‑1 (MQ=255) cAGAAGTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTAc                                                                                                                                            <  7:125946/149‑1 (MQ=255) cAGAAGTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTAc                                                                                                                                            >  7:321468/1‑149 (MQ=255)  aGAAGTATATCAGTGCAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATGCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGACGATACCGTGATGGCGGGCGGCTATAGTGCAACCGACATGGCGGCGGGCGCTATCGGTACt                                                                                                                                           >  5:315611/1‑149 (MQ=255)        atatCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATc                                                                                                                                     >  1:255994/1‑149 (MQ=255)             aGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGTGTTTTGTCGATAAAGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATGGGAACTTCTATCTGGCt                                                                                                                                >  2:355396/1‑149 (MQ=255)             aGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTTTCGGTACTTCTATCTGGCt                                                                                                                                >  8:196094/1‑149 (MQ=255)             aGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCt                                                                                                                                >  4:215744/1‑149 (MQ=255)              gTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCtt                                                                                                                               <  6:153765/149‑1 (MQ=255)                 aaGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCg                                                                                                                            >  5:381774/1‑149 (MQ=255)                            ttattaGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTTTCGGTACTTCTATCTGGCTTCCGGCGATCctct                                                                                                                  <  7:196094/148‑1 (MQ=255)                               ttagcattagcaATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGAGGACGGCTATAGTGCCACCGACATGGCGGCTGTCGCTATCGGTAGTTCTATCTGGCTTCCGGCGATCCTCTTTgg                                                                                                              >  5:119626/1‑149 (MQ=255)                               ttagcattagcaATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTGTTGTCGATAGCGTGATGGCGGGCGGCTATATTGCCACCGTCATGGCTGCGGTCGCTTTGGTTACTTCTATCTGGCTTCCGGCGATCCTCTTTg                                                                                                               <  3:128072/148‑1 (MQ=255)                                    attagcaATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCAc                                                                                                          <  4:14696/148‑1 (MQ=255)                                    attagcaATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCAc                                                                                                          >  3:14696/1‑148 (MQ=255)                                    attagcaATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACg                                                                                                         <  8:374544/149‑1 (MQ=255)                                     ttagcaATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACg                                                                                                         >  1:404104/1‑148 (MQ=255)                                        ggaaTCCCGTTGATTCTCGCTCAAATCGCCCAAACTGCGTTGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGGCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTg                                                                                                               <  6:340364/137‑1 (MQ=255)                                        gcaATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTg                                                                                                               >  5:340364/1‑139 (MQ=255)                                           aTCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATTGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGTTCCTCTTTGGTCACGGActgct                                                                                                  <  1:185054/149‑1 (MQ=255)                                            tCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGActgctg                                                                                                 >  2:409801/1‑149 (MQ=255)                                              ccGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATc                                                                                                                                                 >  2:34811/1‑99 (MQ=255)                                              ccGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATc                                                                                                                                                 <  1:34811/99‑1 (MQ=255)                                                ggTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTgg                                                                                             <  2:404104/149‑1 (MQ=255)                                                ggTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTgg                                                                                             <  8:321468/149‑1 (MQ=255)                                                 gTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTgg                                                                                             >  7:227596/1‑148 (MQ=255)                                                          gcgcAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAAc                                                                                      <  4:190814/146‑1 (MQ=255)                                                          gcgcAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCACTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGcc                                                                                   >  2:94824/1‑149 (MQ=255)                                                             cAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGt                                                                                <  8:413858/149‑1 (MQ=255)                                                               aaTCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTa                                                                              <  3:152907/149‑1 (MQ=255)                                                                     ccAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATcgcg                                                                         >  2:287482/1‑148 (MQ=255)                                                                     ccAAACTGCGATGGGTTTTGTCGATACCGGGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCCCGGACTGCTGCTGGCATTAACGCCGGTTATcgcg                                                                         >  2:282505/1‑148 (MQ=255)                                                                           tGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTa                                                                   >  4:22656/1‑148 (MQ=255)                                                                             cGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCACTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTaaa                                                                 <  1:94824/148‑1 (MQ=255)                                                                              gATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCgg                                                                                 <  7:75226/131‑1 (MQ=255)                                                                              gATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCgg                                                                                 >  8:75226/1‑131 (MQ=255)                                                                               aTGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTg                                                                                                               <  5:118033/100‑1 (MQ=255)                                                                               aTGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTg                                                                                                               >  6:118033/1‑100 (MQ=255)                                                                                tGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGt                                                             >  8:272761/1‑149 (MQ=255)                                                                                  ggTGTTGTCGTTACCGTGATGGAGGGCGGCTATATTGCCAGCGACTTGGCGGCGGTCGCTATCGGTTCTTCTATCTGGCTTCCGGCGATACTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGtt                                                            <  1:355396/148‑1 (MQ=255)                                                                                   gTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCg                                                                                                                            >  6:137729/1‑83 (MQ=255)                                                                                   gTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCg                                                                                                                            <  5:137729/83‑1 (MQ=255)                                                                                    ttttGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCCCGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTaaa                                                                 >  2:278298/1‑141 (MQ=255)                                                                                    ttttGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTaaa                                                                 <  1:278298/141‑1 (MQ=255)                                                                                           gATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGAc                                                   <  3:215744/148‑1 (MQ=255)                                                                                           gATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACg                                                  >  2:75359/1‑149 (MQ=255)                                                                                           gATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGG‑TATCGCGCAATTAAATGGTTCCGGTCGACGt                                                 >  1:8324/1‑149 (MQ=255)                                                                                              aCCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAAt                                                                <  8:196474/132‑1 (MQ=255)                                                                                              aCCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAAt                                                                >  7:196474/1‑132 (MQ=255)                                                                                               ccGTGATGGCGGGCGGCTTTAGTGCCACCGACATGGCGGCGGTCTCTATCGGTACTTCTATCTGGCTTCCTGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGa                                               <  3:358148/148‑1 (MQ=255)                                                                                               ccGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGa                                               <  5:200777/148‑1 (MQ=255)                                                                                               ccGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAg                                              >  2:280428/1‑149 (MQ=255)                                                                                                 gTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTaaa                                                                 >  6:268842/1‑128 (MQ=255)                                                                                                 gTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTaaa                                                                 <  5:268842/128‑1 (MQ=255)                                                                                                    aTGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCAt                                         >  5:230416/1‑149 (MQ=255)                                                                                                        cgggcggCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTgcg                                     <  6:230416/149‑1 (MQ=255)                                                                                                         gggcggCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTACGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTgag                                     >  6:12820/1‑146 (MQ=255)                                                                                                             ggcTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAGCGCCGGGTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCAt                                  >  5:251194/1‑147 (MQ=255)                                                                                                              gcTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAg                               >  1:254672/1‑149 (MQ=255)                                                                                                               cTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATAAgg                              >  5:137928/1‑149 (MQ=255)                                                                                                                tataGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTgc                                      <  8:3991/140‑1 (MQ=255)                                                                                                                tataGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTgc                                      >  7:3991/1‑140 (MQ=255)                                                                                                                 ataGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTg                            >  6:277548/1‑149 (MQ=255)                                                                                                                  taGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCTCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGc                           >  3:400999/1‑149 (MQ=255)                                                                                                                  taGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGt                             <  1:409801/147‑1 (MQ=255)                                                                                                                     tGCCACCGACATGGCGGCGGTCACTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCaa                                                                      >  4:7615/1‑103 (MQ=255)                                                                                                                     tGCCACCGACATGGCGGCGGTCACTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCaa                                                                      <  3:7615/103‑1 (MQ=255)                                                                                                                        cACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAg                     <  3:22656/149‑1 (MQ=255)                                                                                                                         aCCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATACTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAgg                    >  5:231654/1‑149 (MQ=255)                                                                                                                             aCATGGCGGCGGTCGCTATCGGTACTTCAATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTc                <  6:98510/149‑1 (MQ=255)                                                                                                                             aCATGGCGGCGGTCGCAATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGtt                                                            <  5:318623/105‑1 (MQ=255)                                                                                                                             aCATGGCGGCGGTCGCAATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGtt                                                            >  6:318623/1‑105 (MQ=255)                                                                                                                              cATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTct               <  2:255994/149‑1 (MQ=255)                                                                                                                                tGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTctgg             >  3:352348/1‑149 (MQ=255)                                                                                                                                  gcggcgGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGAAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTctggct           <  1:287482/149‑1 (MQ=255)                                                                                                                                        gTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTCTGGCTGGCAgg     >  7:29233/1‑149 (MQ=255)                                                                                                                                           gCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGACGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTCTGGCTGGCAGGttt  >  1:84188/1‑149 (MQ=255)                                                                                                                                           |                                                                                                                                                     CAGAAGTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTCTGGCTGGCAGGTTT  >  NC_000913/1743460‑1743746

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