Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,894,867	G→A	84.9%	S64S (AGC→AGT)	yegH ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,894,867	0	G	A	84.9%	170.0 / 21.5	73	S64S (AGC→AGT)	yegH	putative transporter
Reads supporting (aligned to +/- strand):  ref base G (4/7);  new base A (31/31);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

GTGACCTTCCAGGCGTTCGTGCAGCTCCATGGTGGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATGCTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAG  >  NC_002947/3894723‑3895013                                                                                                                                                 |                                                                                                                                                   gTGACCTTCCAGGCGTTCGTGCAGCTCCATGGTGGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTgg                                                                                                                                                >  8:201594/1‑149 (MQ=255) gTGACCTTCCAGGCGTTCGTGCAGCTCCATGGTGGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACCAAATACTgg                                                                                                                                                >  7:68257/1‑149 (MQ=255)  tGACCTTCCAGGCGTTCGTGCAGCTCCATGGTGGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATGCTGGc                                                                                                                                               <  6:313380/149‑1 (MQ=255)  agACCTTCCAGGCGTTCGTGCAGCTCCATGGTGGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGc                                                                                                                                               >  1:248699/2‑149 (MQ=255)      cTTCCAGGCGTTCGTGCAGCTCCATGGTGGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCaa                                                                                                                                            >  3:146682/1‑148 (MQ=255)        tCCAGGCGTTCGTGCAGCTCCATGGTGGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACa                                                                                                                                          <  5:293769/148‑1 (MQ=255)          cAGGCGTTCGTGCAGCTCCATGGTGGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGc                                                                                                                                       <  2:248699/149‑1 (MQ=255)          cAGGCGTTCGTGCAGCTCCATGGTGGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGc                                                                                                                                       <  6:190495/149‑1 (MQ=255)           aGGCGTTCGTGCAGCTCCATGGTGGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAgg                                                                                                                                        >  5:9577/1‑147 (MQ=255)           aGGCGTTCGTGCAGCTCCATGGTGGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAgg                                                                                                                                        <  6:9577/147‑1 (MQ=255)             gCGTTCGTGCAGCTCCATGGTGGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGACGAACTCTTCGAACAGAGGTGCGGTCAGGGTGACCATCCACGAAATACTGGACAACAGGcc                                                                                                                                      >  6:191757/1‑147 (MQ=255)              cGTTCGTGCAGCTCCATGGTCGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAg                                                                                                                                   <  4:248570/149‑1 (MQ=255)               gTTCGTGCAGCTCCATGGTGGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATGCTGGCCAACAGGCCCAgg                                                                                                                                  >  3:191651/1‑149 (MQ=255)                    tGCAGCTCCATGGTGGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCat                                                                                                                             >  7:170526/1‑149 (MQ=255)                     gCAGCTCCATGGTGGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATGCTGGCCAACAGGCCCAGGCGCat                                                                                                                             >  8:297550/1‑148 (MQ=255)                      cAGCTCCATGGTGGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCatga                                                                                                                           <  6:79333/149‑1 (MQ=255)                       aGCTCCATGGTGGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCatgat                                                                                                                          <  3:218667/149‑1 (MQ=255)                                gtggCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCa                                                                                                                  <  3:156466/148‑1 (MQ=255)                                 tggCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAgg                                                                                                                >  8:105704/1‑149 (MQ=255)                                   gCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCatgat                                                                                                                          <  3:274847/137‑1 (MQ=255)                                   gCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCatgat                                                                                                                          >  4:274847/1‑137 (MQ=255)                                     cTTGAACAACAGGAACACACCACCGAACAGCATGATCATGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATGCTGGCCAACATGCCCAGGCGCATGATGAGCGCCAGGCTCa                                                                                                            <  2:100416/149‑1 (MQ=255)                                     cTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATTATGAGCGCCAGGCt                                                                                                              >  8:252616/1‑147 (MQ=255)                                     cTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACACCGGTGCGGTCAGCGTGACCATCCACGAAATACTGGCCAACATGCCCAGGCGCATGATGAGCGCAAGGCt                                                                                                              >  3:208535/1‑147 (MQ=255)                                        gAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAg                                                                                                           <  7:252616/147‑1 (MQ=255)                                         aacaacAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGc                                                                                                         <  7:201594/148‑1 (MQ=255)                                          acaacaGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGcc                                                                                                        >  6:24535/1‑148 (MQ=255)                                            aacaGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATGCTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAAt                                                                                                     <  3:185734/149‑1 (MQ=255)                                             acaGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAAt                                                                                                     >  2:14907/1‑148 (MQ=255)                                                   acacacCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCcgc                                                                                              >  8:46086/1‑149 (MQ=255)                                                               aCAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCgcgc                                                                                  <  8:68257/149‑1 (MQ=255)                                                                 aGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCt                                                                                 >  4:314641/1‑148 (MQ=255)                                                                 aGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCt                                                                                 >  6:49814/1‑148 (MQ=255)                                                                  gCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGa                                                                               <  4:146682/149‑1 (MQ=255)                                                                   cATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGAt                                                                              <  8:170526/149‑1 (MQ=255)                                                                     tGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATGCTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATgg                                                                            >  6:333630/1‑149 (MQ=255)                                                                     tGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATgg                                                                            >  7:263460/1‑149 (MQ=255)                                                                       atgaGGTCTCGGCCAGTGAAGATCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGTTGACCATACAATAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATggg                                                                           <  4:208535/145‑1 (MQ=255)                                                                        tCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGAt                                                                              <  7:259756/144‑1 (MQ=255)                                                                        tCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGAt                                                                              >  8:259756/1‑144 (MQ=255)                                                                        tCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGggc                                                                         >  5:237796/1‑149 (MQ=255)                                                                            gTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATgggg                                                                          <  3:130456/144‑1 (MQ=255)                                                                            gTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATgggg                                                                          >  4:130456/1‑144 (MQ=255)                                                                               aCGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACa                                                                                                                                          >  4:193326/1‑77 (MQ=255)                                                                               aCGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACa                                                                                                                                          <  3:193326/77‑1 (MQ=255)                                                                               aCGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACt                                                                  >  1:96989/1‑149 (MQ=255)                                                                                    cAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCg                                                             <  6:237796/149‑1 (MQ=255)                                                                                    cAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCg                                                             <  7:105704/149‑1 (MQ=255)                                                                                      gagaAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATGCTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGc                                                           <  5:333630/149‑1 (MQ=255)                                                                                       agaAGCTCTTGCAGAACACTTCGAACAGAGGTGCGGTCAGGGTGACAATCCACGAAATACTGGCCAACAGTCCAAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTAGGc                                                           >  4:94116/1‑148 (MQ=255)                                                                                        gaAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATGCTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCaa                                                         <  4:191651/149‑1 (MQ=255)                                                                                         aaGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAg                                                        <  8:263460/149‑1 (MQ=255)                                                                                               ttGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCAAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGc                                                  >  7:105008/1‑149 (MQ=255)                                                                                                 gCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGa                                                >  7:106279/1‑149 (MQ=255)                                                                                                 gCCGAACACTTCGAACAGAGGTTCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGTGGCGGCAACTTGTCTGCAAGGATgt                                                   <  3:94116/146‑2 (MQ=255)                                                                                                  ccGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGc                                                  >  1:299347/1‑146 (MQ=255)                                                                                                  ccGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGc                                                  <  2:299347/146‑1 (MQ=255)                                                                                                  ccGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGa                                                <  2:96989/148‑1 (MQ=255)                                                                                                       cacTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAaca                                          <  5:49814/149‑1 (MQ=255)                                                                                                          ttCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAATTTGTCGGCAAGGATGGCGATGAacac                                         >  3:19407/1‑147 (MQ=255)                                                                                                          ttCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAacac                                         >  8:201464/1‑147 (MQ=255)                                                                                                          ttCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAacac                                         >  6:172988/1‑147 (MQ=255)                                                                                                               aCAGCGGTGCGGTCAGGGTGACCATCCACGAAATGCTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTg                                  >  8:295796/1‑149 (MQ=255)                                                                                                                 aGCTGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTc                                <  5:191757/149‑1 (MQ=255)                                                                                                                 aGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTc                                <  7:46086/149‑1 (MQ=255)                                                                                                                 aGCGGTGCGGTCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTc                                <  7:201464/149‑1 (MQ=255)                                                                                                                     gtgcggtCAGGGTGACCATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATa                            <  1:14907/149‑1 (MQ=255)                                                                                                                                     cATCCACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAATTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGa            <  4:19407/149‑1 (MQ=255)                                                                                                                                      aTCCACGAAATGCTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAg           <  7:297550/149‑1 (MQ=255)                                                                                                                                         cACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGa         <  8:106279/148‑1 (MQ=255)                                                                                                                                         cACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGAc        >  7:174486/1‑149 (MQ=255)                                                                                                                                         cACGAAATACTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGAc        >  4:31361/1‑149 (MQ=255)                                                                                                                                               aTGCTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAg  <  7:295796/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   GTGACCTTCCAGGCGTTCGTGCAGCTCCATGGTGGCCTTGAACAACAGGAACACACCACCGAACAGCATGATCAGGTCACGGCCAGAGAAGCTCTTGCCGAACACTTCGAACAGCGGTGCGGTCAGGGTGACCATCCACGAAATGCTGGCCAACAGGCCCAGGCGCATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAG  >  NC_002947/3894723‑3895013

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