Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,144,915 G→A 9.8% S108N (AGC→AAC)  PP_2760 → ribose ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,144,9150GA 9.8% 168.3 / 6.7 61S108N (AGC→AAC) PP_2760ribose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (27/28);  new base A (3/3);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.75e-01

TTGCCTTCGGCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGACCGAAAACACCGT  >  NC_002947/3144771‑3145062
                                                                                                                                                |                                                                                                                                                    
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                                                                                                                                             >  6:269400/1‑144 (MQ=255)
         gCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGCTGACAACATcgccga                                                                                                                                             <  5:269400/145‑2 (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                                                <  5:203809/131‑1 (MQ=255)
                                                                                                                   gAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGa                                                >  6:203809/1‑131 (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                               >  4:391798/1‑131 (MQ=255)
                                                                                                                                    gCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCggg                               <  3:391798/131‑1 (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)
                                                                                                                                               aGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCCGGCTGGAACTGGTGCTGACCGAAAACACCg   <  1:33179/148‑1 (MQ=255)
                                                                                                                                               aGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCAGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGACCGAAAACACCGt  <  4:297812/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TTGCCTTCGGCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCCAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCC‑ATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGACCGAAAACACCGT  >  NC_002947/3144771‑3145062

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