Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,610,644	C→T	6.8%	S20L (TCG→TTG) ‡	emrE →	small multidrug resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,610,644	0	C	T	6.8%	261.2 / 6.0	88	S20L (TCG→TTG) ‡	emrE	small multidrug resistance protein
Reads supporting (aligned to +/- strand):  ref base C (38/44);  new base T (3/3);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

TCAACCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATT  >  NC_002947/5610499‑5610789                                                                                                                                                  |                                                                                                                                                     tCAACCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGAt                                                                                                                                                   <  3:235690/149‑1 (MQ=255)      cTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑gg                                                                                                                                             >  7:91179/1‑149 (MQ=255)      cTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑gg                                                                                                                                             >  2:19596/1‑149 (MQ=255)      cTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑gg                                                                                                                                             <  6:184476/149‑1 (MQ=255)      cTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑gg                                                                                                                                             <  5:109586/149‑1 (MQ=255)       tATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCAATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGc                                                                                                                                            >  2:79167/1‑149 (MQ=255)            cccATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTTCATCGA‑GGCGGTgg                                                                                                                                       <  7:292676/149‑2 (MQ=255)            cccATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTTCATCGA‑GGCGGTgg                                                                                                                                       <  4:10182/149‑2 (MQ=255)             ccATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCgg                                                                                                                                          <  4:238674/145‑1 (MQ=255)                tatCGATACGTAATGACTTTACCCCATCAAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAgg                                                                                                                                   <  4:277942/149‑1 (MQ=255)                  tCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTTCATCGA‑GGCGGTggcgatga                                                                                                                                 <  5:363693/149‑8 (MQ=255)                   cGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCt                                                                                                                                >  2:264807/1‑149 (MQ=255)                    gATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCt                                                                                                                                <  3:195008/148‑1 (MQ=255)                     aTAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGa                                                                                                                              <  5:45552/149‑1 (MQ=255)                      tAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGa                                                                                                                              <  7:260714/148‑1 (MQ=255)                        ggTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGc                                                                                                                            <  5:165774/148‑1 (MQ=255)                         gTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCAc                                                                                                                          <  8:55028/149‑1 (MQ=255)                          tAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCt                                                                                                                                >  2:169902/1‑142 (MQ=255)                          tAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCt                                                                                                                                <  1:169902/142‑1 (MQ=255)                               aCTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGtt                                                                                                                    >  5:324922/1‑149 (MQ=255)                                 tttACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑Gc                                                                                                                 >  6:324196/1‑149 (MQ=255)                                 tttACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑Gc                                                                                                                 >  3:24082/1‑149 (MQ=255)                                 tttACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATCAA‑GGCGGTGAAAGGCCTGAGCACGCCCTT‑Gc                                                                                                                 >  8:223507/1‑149 (MQ=255)                                    aCCCCATCCAGGCAAGGTCTAGAATGCCCGCTGAGTAACGGAGGATGCTGCCTTGAATGCCTATACCTACCTCGGCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCa                                                                                                               <  7:183028/148‑1 (MQ=255)                                    aCCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCa                                                                                                               <  7:188261/148‑1 (MQ=255)                                    aCCCCATCCAGCCAACGTCTAGAATGCCCTCTCAGTAACGGAGGATCCTGCCATGAATGCCTCTACCTACCTCGCCAATGCCATCTGCGCCGAAGTCAGCGCCACCGCCTCTATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCa                                                                                                               <  6:315797/148‑1 (MQ=255)                                            cAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCActgctgct                                                                                                       <  1:340457/148‑1 (MQ=255)                                             aGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGa                                                                                                     >  1:243955/1‑149 (MQ=255)                                               gCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGAtg                                                                                                   <  4:24082/149‑1 (MQ=255)                                                 aaCGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGAtggt                                                                                                 >  1:203505/1‑149 (MQ=255)                                                 aaCGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGAtggt                                                                                                 >  8:178525/1‑149 (MQ=255)                                                  aCGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAACGGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGAtggt                                                                                                 <  1:264807/149‑1 (MQ=255)                                                   cGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGAtggtgg                                                                                               >  8:143793/1‑149 (MQ=255)                                                     tCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGt‑t                                                                                            <  7:178525/149‑1 (MQ=255)                                                       tAGAAAGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGAT‑CCCACTGCTGCTGATGGTGGT‑Tgg                                                                                          >  8:92455/1‑149 (MQ=255)                                                        aGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑Tgg                                                                                          <  7:223507/148‑1 (MQ=255)                                                         gAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑Tg                                                                                           <  7:143793/146‑1 (MQ=255)                                                         gAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGc                                                                                         >  2:368041/1‑148 (MQ=255)                                                              cccGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATgg                                                                                    <  1:19596/148‑1 (MQ=255)                                                               ccGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCa                                                                                  >  2:127504/1‑149 (MQ=255)                                                                     aGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCg                                                                             <  8:91179/148‑1 (MQ=255)                                                                       tAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGNTNNTGGT‑TGGCTATGGCATTGCGTTc                                                                          >  8:279971/1‑149 (MQ=255)                                                                       tAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTc                                                                          >  8:275030/1‑149 (MQ=255)                                                                         aCGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCt                                                                         <  2:243955/148‑1 (MQ=255)                                                                            gaggaTCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGAt                                                                     <  5:263911/149‑1 (MQ=255)                                                                            gaggaTCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGAt                                                                     >  3:221448/1‑149 (MQ=255)                                                                                 tCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCAGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGa                                                                <  4:282008/149‑1 (MQ=255)                                                                                  ccTGCCATGAATGCCTATACCTACCTCGCCATTGTCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGAc                                                               >  5:82214/1‑149 (MQ=255)                                                                                  ccTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGAc                                                               >  2:76707/1‑149 (MQ=255)                                                                                    tGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGc                                                             >  7:252023/1‑149 (MQ=255)                                                                                    tGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGc                                                             >  6:286002/1‑149 (MQ=255)                                                                                     gCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCt                                                            >  1:301143/1‑149 (MQ=255)                                                                                      ccATGAATGCCTATACCTACCTCGCCATTGTCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGc                                                             <  6:82214/147‑1 (MQ=255)                                                                                       cATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCtgg                                                          <  8:252023/149‑1 (MQ=255)                                                                                        aTGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCtggt                                                         <  6:324922/149‑1 (MQ=255)                                                                                         tGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCtggtg                                                        <  5:343783/149‑1 (MQ=255)                                                                                         tGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCtggtg                                                        <  1:79167/149‑1 (MQ=255)                                                                                          gAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCtggtgg                                                       >  3:112055/1‑149 (MQ=255)                                                                                           caTGCCTATCCCTACCTCGCCATTCCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTggtggt                                                      <  4:45060/148‑1 (MQ=255)                                                                                                cTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgca                                                 <  2:131246/149‑1 (MQ=255)                                                                                                cTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCg                                                   >  8:81439/1‑147 (MQ=255)                                                                                                   tacctaccTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGTCTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgcag                                                <  1:76707/148‑1 (MQ=255)                                                                                                      ctacctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTc                                           <  1:368041/149‑1 (MQ=255)                                                                                                      ctacctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTc                                           >  3:182425/1‑149 (MQ=255)                                                                                                         cctCGCCATTGCCCTCTGCGCCTAAGTGCTCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCg                                         <  6:18374/148‑1 (MQ=255)                                                                                                            cGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTgg                                     <  7:298742/149‑1 (MQ=255)                                                                                                             gCCATTGCCATCTTCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTgg                                     <  5:324196/148‑1 (MQ=255)                                                                                                             gCCATTGCCATCTGCGCCGAAGTCCTCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTggg                                    <  7:279971/149‑1 (MQ=255)                                                                                                             gCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTggg                                    <  7:275030/149‑1 (MQ=255)                                                                                                                            gcCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTTAGCCAATGCTGCTGAGGGTGGT‑TGGCTATGGAATTGCGTGCTGGATGATGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCa                      >  1:289618/1‑149 (MQ=255)                                                                                                                             ccGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATc                    >  1:246179/1‑149 (MQ=255)                                                                                                                                  gTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTc               >  3:151244/1‑149 (MQ=255)                                                                                                                                     aTCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCggg            >  1:197263/1‑149 (MQ=255)                                                                                                                                     aTCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCggg            <  4:112055/149‑1 (MQ=255)                                                                                                                                     aTCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCggg            >  7:311022/1‑149 (MQ=255)                                                                                                                                      tCGCCACCGGCTCGATGAA‑GGCGGTGAAAGGCCTGCGCACGCCGTT‑GCCACTGCTGGTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGGATTCCGGTGGGGATTGCCTACGCCATCTGGTCgggg           <  7:92455/149‑1 (MQ=255)                                                                                                                                      tCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTggg                                    >  6:50357/1‑124 (MQ=255)                                                                                                                                      ctttcACCGCCTTCATCGA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCgggg           >  2:6062/5‑149 (MQ=255)                                                                                                                                       tttcACCGCCTTCATCGA‑GGCGGTGGAAGGCATGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGc          >  1:190789/4‑149 (MQ=255)                                                                                                                                       ttccACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGc          >  7:225840/3‑149 (MQ=255)                                                                                                                                       cGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCt                            >  4:186610/1‑131 (MQ=255)                                                                                                                                       cGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCt                            <  3:186610/131‑1 (MQ=255)                                                                                                                                         tcACCGCCTTCATCGA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTg        >  7:14390/2‑149 (MQ=255)                                                                                                                                         ccACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTg        >  6:264271/1‑149 (MQ=255)                                                                                                                                         ccACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTg        >  2:69735/1‑149 (MQ=255)                                                                                                                                            ccGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTggg      >  5:253992/1‑148 (MQ=255)                                                                                                                                               ccTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGAtt  <  7:84716/149‑1 (MQ=255)                                                                                                                                               ccTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGAtt  <  1:127504/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     TCAACCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAA‑GGCGGTGAAAGGCCTGAGCACGCCGTT‑GCCACTGCTGCTGATGGTGGT‑TGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATT  >  NC_002947/5610499‑5610789

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