Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,843,535	C→T	G420D (GGC→GAC)	gacS ←	sensor protein GacS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,843,535	0	C	T	88.5%	197.4 / 14.5	79	G420D (GGC→GAC)	gacS	sensor protein GacS
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base T (33/36);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGAT  >  NC_002947/1843393‑1843683                                                                                                                                                 |                                                                                                                                                      tGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTcgc                                                                                                                                                 <  7:17531/149‑1 (MQ=255) tGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTcgc                                                                                                                                                 >  5:369239/1‑149 (MQ=255)   aGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCgg                                                                                                                                               <  1:187977/149‑1 (MQ=255)       gTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTg                                                                                                                                            >  4:15994/1‑148 (MQ=255)        tCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGaa                                                                                                                                          <  4:305312/149‑1 (MQ=255)          aGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTGCTGCACGCTGATGCGCAACTCGCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAAc                                                                                                                                         <  1:335927/149‑1 (MQ=255)           gCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACtt                                                                                                                                       <  6:7685/149‑1 (MQ=255)                gCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATg                                                                                                                                   <  7:239979/148‑1 (MQ=255)                 cTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCgg                                                                                                                                               >  1:98230/1‑135 (MQ=255)                 cTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCgg                                                                                                                                               <  2:98230/135‑1 (MQ=255)                  tGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCg                                                                                                                                >  2:244606/1‑149 (MQ=255)                      ggCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACTCTGATGCGCAACT‑GCGCGCATTCTTCGTGCTCGTC‑TT‑CAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTCAACTTGATGGCGTTGCt                                                                                                                           <  4:350127/149‑1 (MQ=255)                      ggCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGc                                                                                                                            >  7:272816/1‑149 (MQ=255)                      ggCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGc                                                                                                                            <  2:114243/149‑1 (MQ=255)                      ggCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TT‑CAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCt                                                                                                                           <  4:356050/149‑1 (MQ=255)                      ggCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TT‑CAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCt                                                                                                                           <  4:350124/149‑1 (MQ=255)                        ccTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCt                                                                                                                           >  7:136074/1‑148 (MQ=255)                         cTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCCTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGa                                                                                                                         <  5:258016/149‑1 (MQ=255)                           ggAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GAGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGAc                                                                                                                        >  2:68632/1‑148 (MQ=255)                            gAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGAc                                                                                                                        >  6:18278/1‑147 (MQ=255)                            gAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGAc                                                                                                                        >  1:34456/1‑147 (MQ=255)                                aGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGtt                                                                                                                  >  2:161425/1‑149 (MQ=255)                                 gtgtGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTg                                                                                                                 <  5:18278/149‑1 (MQ=255)                                 gtgaGCGCACAGCCTGCTGTGACAGGCAGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGGCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGGGTCTTCGCGGGTGAACTTGATGGCGTTGCTGCCCAGGTTg                                                                                                                 <  6:351667/149‑1 (MQ=255)                                  tgtgCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTgg                                                                                                                >  4:194387/1‑149 (MQ=255)                                    tgCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGt                                                                                                               >  5:184240/1‑148 (MQ=255)                                       gcACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAg                                                                                                            >  5:407566/1‑148 (MQ=255)                                       gcACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAg                                                                                                            >  3:354646/1‑148 (MQ=255)                                        cacaTCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GTGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGa                                                                                                          <  1:68632/149‑1 (MQ=255)                                         acaTCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGAt                                                                                                         >  5:156863/1‑149 (MQ=255)                                         acaTCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGAt                                                                                                         <  3:15994/149‑1 (MQ=255)                                              cTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGc                                                                                                     <  2:39123/148‑1 (MQ=255)                                              cTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGc                                                                                                     <  3:393667/148‑1 (MQ=255)                                                 aGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCtt                                                                                                   <  8:272816/147‑1 (MQ=255)                                                    tGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGcc                                                                                              >  5:266698/1‑149 (MQ=255)                                                      aCAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGc                                                                                            <  7:10375/149‑1 (MQ=255)                                                        aGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCaa                                                                                          >  4:17209/1‑149 (MQ=255)                                                        aGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCaa                                                                                          >  1:413996/1‑149 (MQ=255)                                                         ggCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑c                                                                                        <  6:407566/149‑1 (MQ=255)                                                             gATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGt                                                                                    <  6:156863/149‑1 (MQ=255)                                                                  cGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGcc                                                                                >  7:343032/1‑148 (MQ=255)                                                                     tgtCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATAGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTAc                                                                             >  1:324786/1‑148 (MQ=255)                                                                       tCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAg                                                                          <  4:354646/149‑1 (MQ=255)                                                                         cTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGc                                                                         <  4:32683/148‑1 (MQ=255)                                                                         cTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGc                                                                         >  2:401559/1‑148 (MQ=255)                                                                          tGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCg                                                                        >  5:391269/1‑148 (MQ=255)                                                                                 tGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTAcc                                                                            <  8:350289/137‑1 (MQ=255)                                                                                 tGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTAcc                                                                            >  7:344577/1‑137 (MQ=255)                                                                                 tGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTAcc                                                                            >  7:350289/1‑137 (MQ=255)                                                                                 tGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTAcc                                                                            <  8:344577/137‑1 (MQ=255)                                                                                 tGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGgt                                                                >  2:370739/1‑149 (MQ=255)                                                                                  gATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTCTTTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCt                                                                                                    >  1:50545/1‑114 (MQ=255)                                                                                  gATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCt                                                                                                    <  2:50545/113‑1 (MQ=255)                                                                                   aTGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGgt                                                                <  8:343032/147‑1 (MQ=255)                                                                                   aTGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGgt                                                                <  1:370739/147‑1 (MQ=255)                                                                                        cAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTcgc                                                         >  2:331377/1‑149 (MQ=255)                                                                                          aCT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCgg                                                       <  6:369239/149‑1 (MQ=255)                                                                                           cT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCgg                                                       >  3:34751/1‑148 (MQ=255)                                                                                                gcgcAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCtt                                                                                                   <  3:223553/101‑1 (MQ=255)                                                                                                gcgcAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCtt                                                                                                   >  4:223553/1‑101 (MQ=255)                                                                                                  gcAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCa                                                <  3:17209/149‑1 (MQ=255)                                                                                                  gcAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCa                                                <  1:161425/149‑1 (MQ=255)                                                                                                        ttCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCt                                            >  2:309098/1‑147 (MQ=255)                                                                                                          cGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTc                                           <  6:266698/146‑1 (MQ=255)                                                                                                               tCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTc                                   >  7:43428/1‑149 (MQ=255)                                                                                                                 gTC‑TTCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATc                                                                                                        <  2:45804/79‑1 (MQ=255)                                                                                                                 gTC‑TTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCa                                 <  2:34456/149‑1 (MQ=255)                                                                                                                 gTC‑TTCCAGCATGGCGCGCACGACGATGGTGACTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATc                                                                                                        >  1:45804/1‑79 (MQ=255)                                                                                                                     ttCCAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGATCAGCAATTCCAgct                              >  4:48970/1‑149 (MQ=255)                                                                                                                     ttCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAAGTCCAgct                              >  2:137289/1‑149 (MQ=255)                                                                                                                       ccAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACCGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgct                              <  1:331377/148‑1 (MQ=255)                                                                                                                       ccAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgctgc                            >  2:94252/1‑149 (MQ=255)                                                                                                                       ccAGCATGGCGCGCACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgc                               <  1:244606/146‑1 (MQ=255)                                                                                                                        cAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgctgc                            <  6:184240/148‑1 (MQ=255)                                                                                                                                 gcgcACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTggg                  <  4:392230/149‑1 (MQ=255)                                                                                                                                    cACGACGATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGggcgg               >  5:361630/1‑149 (MQ=255)                                                                                                                                     acgacgATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGggcggc              >  6:376363/1‑149 (MQ=255)                                                                                                                                     acgacgATGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGggcggc              >  3:376118/1‑149 (MQ=255)                                                                                                                                           aTGGTGTCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGc        <  1:401559/149‑1 (MQ=255)                                                                                                                                                 ccTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGAt  >  6:66524/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                      TGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACT‑GCGCGCAGTCTTCGTGCTCGTC‑TTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAA‑CGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGAT  >  NC_002947/1843393‑1843683

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