Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,144,912	T→A	6.9%	V107D (GTC→GAC)	PP_2760 →	ribose ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,144,912	0	T	A	6.9%	165.7 / 2.9	58	V107D (GTC→GAC)	PP_2760	ribose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (26/28);  new base A (2/2);  total (28/30)
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.34e-01

GCAGCCTTGCCTTCGGCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGACCGAAAACA  >  NC_002947/3144765‑3145058                                                                                                                                                    |                                                                                                                                                    gCAGCCTTGCCTTCGGCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTc                                                                                                                                                    <  5:10632/149‑1 (MQ=255) gCAGCCTTGCCTTCGGCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTc                                                                                                                                                    <  2:376805/149‑1 (MQ=255)       ttgccttCGGCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATc                                                                                                                                              >  2:1442/1‑149 (MQ=255)          ccttCGGCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGCTGACAACATcgc                                                                                                                                            <  5:214991/148‑1 (MQ=255)              cGGCCTGCCCACCGCGATCATCGGCGGCGGCTCGAACGTGGGGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCGGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCCTCGCcgcg                                                                                                                                        <  4:345273/148‑1 (MQ=255)               ggCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGAc                                                                                                                                      >  6:62696/1‑149 (MQ=255)                gCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGCTGACAACATcgccga                                                                                                                                         <  5:269400/145‑2 (MQ=255)                gCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGCTGACAACATcgccga                                                                                                                                         >  6:269400/1‑144 (MQ=255)                       cACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCa                                                                                                                              >  7:48140/1‑149 (MQ=255)                          ggTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCActg                                                                                                                           <  1:232384/149‑1 (MQ=255)                                            cTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGc                                                                                                        >  1:406152/1‑149 (MQ=255)                                                       tgtCCGGTGGGCTCGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCtg                                                                                             <  6:394773/149‑1 (MQ=255)                                                         tCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCtggt                                                                                           <  5:211455/149‑1 (MQ=255)                                                            ggTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCtggtggt                                                                                        <  5:140525/149‑1 (MQ=255)                                                                        cTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGttt                                                                               <  7:277791/146‑1 (MQ=255)                                                                        cTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGttt                                                                               >  8:277791/1‑146 (MQ=255)                                                                         tGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGt                                                                           >  7:325190/1‑149 (MQ=255)                                                                           tCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTcc                                                                             <  7:15660/145‑1 (MQ=255)                                                                           tCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGttt                                                                         <  5:62696/149‑1 (MQ=255)                                                                               tGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTgccg                                                                     <  2:406152/149‑1 (MQ=255)                                                                                   ggcCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCcgccgc                                                                 >  7:312731/1‑149 (MQ=255)                                                                                      ctaCCTCGGGTTGTGCGTTGCGGTCTATGCCAGCCTGAACAACGGCGGTATCGGGGATGTTGTCAGCATCGCCGCGACGCTTGGCTCTGGCCTGGCC‑ATCGGCACGCTGATTGGGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCgctgc                                                              <  3:282845/147‑1 (MQ=255)                                                                                        aCCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGctgctg                                                            <  6:405204/149‑1 (MQ=255)                                                                                        aCCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGctgctg                                                            <  6:47473/149‑1 (MQ=255)                                                                                         ccTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTgg                                                           <  1:1442/149‑1 (MQ=255)                                                                                         ccTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTgg                                                           >  6:216041/1‑149 (MQ=255)                                                                                           tcaggTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTAGTCGGTTTCCGTTTGCCGCCGCTGCTGGCg                                                         >  4:23694/4‑149 (MQ=255)                                                                                                 tGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTg                                                   >  4:182841/1‑149 (MQ=255)                                                                                                  gTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTg                                                   <  4:280555/148‑1 (MQ=255)                                                                                                   tCCGCTGCGGTCTATGCCTGGCTGAACCACGCCGGTTTCTGCGATGTTGTCAGCATCGCCGCGCCCCTTGGCAAGGGCCTGGCCTATCGGCACGCGGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTg                                                   <  1:324984/148‑1 (MQ=255)                                                                                                   tCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGc                                                 >  7:43351/1‑149 (MQ=255)                                                                                                   tCCGCTGCGGTCTATGCCAACCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGc                                                 >  2:33179/1‑149 (MQ=255)                                                                                                     cGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCa                                               >  8:174996/1‑149 (MQ=255)                                                                                                       cTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCAc                                              >  2:403884/1‑148 (MQ=255)                                                                                                          cGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAACATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGAt                                           >  2:134773/1‑148 (MQ=255)                                                                                                                 tGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGcc                                                >  7:125351/1‑136 (MQ=255)                                                                                                                 tGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGcc                                                <  8:125351/136‑1 (MQ=255)                                                                                                                     aGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGc                               <  8:24155/149‑1 (MQ=255)                                                                                                                      gCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGcc                              >  5:380237/1‑149 (MQ=255)                                                                                                                      gCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGcc                              >  5:3887/1‑149 (MQ=255)                                                                                                                       ccTGAACAACGCCGGTTTCGGCGATGCTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCg                             >  1:73715/1‑149 (MQ=255)                                                                                                                         tGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCggg                           >  7:14980/1‑149 (MQ=255)                                                                                                                          gAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGa                                            >  6:203809/1‑131 (MQ=255)                                                                                                                          gAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGa                                            <  5:203809/131‑1 (MQ=255)                                                                                                                               acGCCGGTTTCGGCGATGCTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGaa                     <  2:73715/149‑1 (MQ=255)                                                                                                                                cGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCAGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAAc                    >  3:297812/1‑149 (MQ=255)                                                                                                                                 gCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACt                   >  8:62966/1‑149 (MQ=255)                                                                                                                                  ccGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTg                  >  3:36204/1‑149 (MQ=255)                                                                                                                                  ccGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTg                  <  8:325190/149‑1 (MQ=255)                                                                                                                                     gTTTCGGCGATCTGCTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTg               <  1:344177/149‑1 (MQ=255)                                                                                                                                          ggCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGAc          <  8:312731/149‑1 (MQ=255)                                                                                                                                          ggCGATGTTGTCAACATCGCCGCGACGCTTGGCACTGGCATGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGa           <  1:134773/148‑1 (MQ=255)                                                                                                                                          ggCGATGCTGACAACATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGAc          >  1:65707/1‑149 (MQ=255)                                                                                                                                           gCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCtggtg                                 <  2:164975/125‑1 (MQ=255)                                                                                                                                           gCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCtggtg                                 >  1:164975/1‑125 (MQ=255)                                                                                                                                           gCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCggg                           <  3:391798/131‑1 (MQ=255)                                                                                                                                           gCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCggg                           >  4:391798/1‑131 (MQ=255)                                                                                                                                              aTGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGACCGaa      <  1:403884/149‑1 (MQ=255)                                                                                                                                                gttgtCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGAc          <  3:324665/143‑1 (MQ=255)                                                                                                                                                gttgtCAGCATCGCCGCGACGCTTGGAACTCGCCTGACC‑ATCGGCACGCTGAATGCGCTAGTCGTGGTCGGTATCAGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGAAGAACCTGGTGGCCGCGATGGAACTGGTGCTGAc          >  4:324665/1‑143 (MQ=255)                                                                                                                                                  tgtCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGACCGAAAaca  >  3:229826/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCAGCCTTGCCTTCGGCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGACCGAAAACA  >  NC_002947/3144765‑3145058

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