Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	508,954:1	+C	7.5%	coding (57/1890 nt)	estP ←	esterase EstP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	508,954	1	.	C	7.5%	258.3 / 14.0	67	coding (57/1890 nt)	estP	esterase EstP
Reads supporting (aligned to +/- strand):  ref base . (35/27);  new base C (3/2);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

CGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTCGTTTCC  >  NC_002947/508810‑509099                                                                                                                                                  |                                                                                                                                                   cGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTg                                                                                                                                                 >  4:67310/1‑148 (MQ=255)   gggCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAg                                                                                                                                              <  3:124022/149‑1 (MQ=255)      cGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAggcc                                                                                                                                           <  3:342002/149‑1 (MQ=255)          gTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAggccaggg                                                                                                                                       >  1:356516/1‑149 (MQ=255)          gTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAggccaggg                                                                                                                                       <  3:53042/149‑1 (MQ=255)           ttGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGGCCTGCa‑gccaggcg                                                                                                                                      <  6:188024/149‑3 (MQ=255)              ccGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGGCCTGCa‑gccaggcgtt                                                                                                                                    <  3:181084/148‑5 (MQ=255)              ccGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAggccagggccag                                                                                                                                   <  7:144648/149‑1 (MQ=255)                gTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGt                                                                                                                                  >  8:19975/1‑148 (MQ=255)                   ggCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGaa                                                                                                                               <  2:77764/148‑1 (MQ=255)                      gTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGc                                                                                                                             >  6:242034/1‑147 (MQ=255)                       tCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAg                                                                                                                          <  5:273083/149‑1 (MQ=255)                       tCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAg                                                                                                                          <  6:279821/149‑1 (MQ=255)                         aCGGGTGGTGTAACGCGCGGCTGGGGTACCGCCAACAAGATCGGGAAAGTGCCCTGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAggg                                                                                                                        <  2:99253/149‑1 (MQ=255)                                      cgcgcgCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATaa                                                                                                          >  6:209928/1‑149 (MQ=255)                                                 ggTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGGGGGTAAAGCGC‑AATAACGGGGCTTTTc                                                                                               <  1:354658/149‑1 (MQ=255)                                                 ggTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTc                                                                                               <  2:356516/149‑1 (MQ=255)                                                   tACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGGTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGc                                                                                             <  6:241505/149‑1 (MQ=255)                                                        ccAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCgg                                                                                        >  7:173651/1‑149 (MQ=255)                                                          aacaaGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCggg                                                                                       >  3:169421/1‑148 (MQ=255)                                                           acaaGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGa                                                                                     <  8:173651/149‑1 (MQ=255)                                                                 tCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACCCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGTATCGGCGAACATcc                                                                               >  4:334525/1‑149 (MQ=255)                                                                    ggAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCttt                                                                            >  5:365366/1‑149 (MQ=255)                                                                     gAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCtttt                                                                           >  4:95290/1‑149 (MQ=255)                                                                      aaaCTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACCTCCTTTTc                                                                          >  7:79095/1‑149 (MQ=255)                                                                       aaCTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCt                                                                         >  2:268730/1‑149 (MQ=255)                                                                        aCTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCtt                                                                        >  3:331207/1‑149 (MQ=255)                                                                         cTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCttt                                                                       <  6:365366/149‑1 (MQ=255)                                                                          tGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTc                                                                          <  8:211917/145‑1 (MQ=255)                                                                          tGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTc                                                                          >  7:211917/1‑145 (MQ=255)                                                                             ccGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGttt                                                                   >  5:193269/1‑149 (MQ=255)                                                                                     ggCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAgg                                                           <  8:241083/149‑1 (MQ=255)                                                                                     ggCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAgg                                                           >  7:241083/1‑149 (MQ=255)                                                                                     ggCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAg                                                            <  5:242034/148‑1 (MQ=255)                                                                                             tGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGCGCAAGcc                                                    >  7:318825/1‑148 (MQ=255)                                                                                                cGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATc                                                >  6:56416/1‑149 (MQ=255)                                                                                                   cAAACACGAGCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGAc                                             <  7:122072/149‑1 (MQ=255)                                                                                                        acGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAg                                        >  4:215205/1‑149 (MQ=255)                                                                                                        acGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAg                                        <  3:313831/149‑1 (MQ=255)                                                                                                        acGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAg                                        >  7:139467/1‑149 (MQ=255)                                                                                                         cGATCAGGGTTGAATAGGGCGAAGGGCCGTCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGc                                       <  3:67310/149‑1 (MQ=255)                                                                                                              aGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAAtc                                  <  6:193269/149‑1 (MQ=255)                                                                                                                 gTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAAtcgtc                               <  5:361536/149‑1 (MQ=255)                                                                                                                 gTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAAtcgtc                               >  6:54691/1‑149 (MQ=255)                                                                                                                  ttGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGGACAAATAAGGGGGGTTTAGGATGCGGGGACCTGCATTTGGTTTTTTTTTTTATGGGC‑AGCCCATGCACCTTTGCAAAttcgag                              >  1:311191/1‑144 (MQ=255)                                                                                                                      aTAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCt                          >  5:185533/1‑149 (MQ=255)                                                                                                                        aGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTc                        >  3:332416/1‑149 (MQ=255)                                                                                                                          ggCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGc                      >  4:17586/1‑149 (MQ=255)                                                                                                                          ggCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTCTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGc                      >  3:380045/1‑149 (MQ=255)                                                                                                                           gCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCa                                                  <  7:22948/120‑1 (MQ=255)                                                                                                                             gAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTcc                   <  7:19975/149‑1 (MQ=255)                                                                                                                             gAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTcc                   >  4:113051/1‑149 (MQ=255)                                                                                                                              aaGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTccc                  <  4:332416/149‑1 (MQ=255)                                                                                                                                gggCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTccc                  >  5:153603/1‑147 (MQ=255)                                                                                                                                 ggccggccAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAg               >  1:128616/1‑149 (MQ=255)                                                                                                                                  gccggccAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCgttgtt                                                                    >  2:361565/1‑95 (MQ=255)                                                                                                                                   ccggccAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTc             >  4:194384/1‑149 (MQ=255)                                                                                                                                       actggcccTGGCCTGCA‑GCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCAtt         >  1:189043/7‑149 (MQ=255)                                                                                                                                        cAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTc        <  4:169421/149‑1 (MQ=255)                                                                                                                                         tggcccTGGCCTGCA‑GCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTCg       >  2:20897/5‑149 (MQ=255)                                                                                                                                         aaCGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGAAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTCg       >  2:346013/1‑149 (MQ=255)                                                                                                                                          atgCCGGG‑CTGCAGGCCAGGGCCAGTGAATCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTCGt      <  7:250665/147‑1 (MQ=255)                                                                                                                                          aCGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTCGt      >  7:356922/1‑149 (MQ=255)                                                                                                                                           gcccTGGCCTGCA‑GCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTCGtt     >  1:133187/3‑149 (MQ=255)                                                                                                                                           cGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTCGtt     <  3:215205/149‑1 (MQ=255)                                                                                                                                             ccTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTCGTTTc   >  7:159404/1‑149 (MQ=255)                                                                                                                                              cTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTCGTTTcc  >  6:264885/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGC‑AATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTCGTTTCC  >  NC_002947/508810‑509099

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