Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,557,034 A→C 6.8% L148R (CTG→CGG)  gntZ ← 6‑phosphogluconate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,557,0340AC 6.8% 213.2 / 4.8 70L148R (CTG→CGG) gntZ6‑phosphogluconate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (33/32);  new base C (2/3);  total (35/35)
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

CTCGATGCCGTTATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATGTCGTCCTT  >  NC_002947/4556889‑4557182
                                                                                                                                                  |                                                                                                                                                    
cTCGATGCCGTTATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCg                                                                                                                                                   <  8:317856/149‑1 (MQ=255)
 tCGATGCCGTTATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCgg                                                                                                                                                  >  6:296201/1‑149 (MQ=255)
   gATGCCGTTATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTGGAACCGctg                                                                                                                                                  <  7:302804/147‑3 (MQ=255)
    aTGCCGTTATGCAACATTTTCACGTTCTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTc                                                                                                                                               <  7:80117/149‑1 (MQ=255)
        cGTTATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGAGCCTTGAACAGAGCTTCCAgg                                                                                                                                           >  1:274036/1‑149 (MQ=255)
         gTTATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTGGAACCGCTGTTCAAGgc                                                                                                                                          <  5:241948/149‑1 (MQ=255)
             tGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTc                                                                                                                                      <  7:103121/149‑1 (MQ=255)
               cACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCgg                                                                                                                                                  >  7:119185/1‑135 (MQ=255)
               cACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCgg                                                                                                                                                  <  8:119185/135‑1 (MQ=255)
                aCCATTTTCACGTAGGGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGaa                                                                                                                                   <  4:248236/149‑1 (MQ=255)
                 ccATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAAGCGCCTTGAACAGCGGTTCCAGGCGCTCGaa                                                                                                                                   <  1:41761/149‑1 (MQ=255)
                   aTTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAacac                                                                                                                                <  1:5707/149‑1 (MQ=255)
                    ttttCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGTCGCC‑AGCGCCTTGAACAGCGGTTcc                                                                                                                                              >  3:270825/1‑134 (MQ=255)
                    ttttCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTcc                                                                                                                                              <  4:270825/134‑1 (MQ=255)
                           gTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCtttt                                                                                                                        >  8:313079/1‑149 (MQ=255)
                              gTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTcg                                                                                                                      <  2:274036/148‑1 (MQ=255)
                                 cccGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTc                                                                                                                                               <  8:28034/120‑1 (MQ=255)
                                 cccGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTc                                                                                                                                               >  7:28034/1‑120 (MQ=255)
                                   cGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGt                                                                                                                              >  5:81642/1‑135 (MQ=255)
                                   cGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGt                                                                                                                              <  6:81642/135‑1 (MQ=255)
                                              ggACCGGCATGGATGTAGCCGTGCTCGGCACGCGGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCa                                                                                                     <  5:274214/149‑1 (MQ=255)
                                                  cGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAg                                                                                                 >  2:369599/1‑149 (MQ=255)
                                                  cGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAg                                                                                                 >  2:90514/1‑149 (MQ=255)
                                                          aTGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTc                                                                                         >  5:245008/1‑149 (MQ=255)
                                                            gTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTc                                                                                         >  8:217574/1‑147 (MQ=255)
                                                             tAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCCGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGAAGCCGCGGTCCAg                                                                                      >  5:20926/1‑149 (MQ=255)
                                                             gaGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCa                                                                                       <  7:291699/147‑1 (MQ=255)
                                                              aGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCCGGCGCTCGAACAGGTCGGTTTCGCCGCCGCTCATCATGCAGTAGCCGCGGTCCAgg                                                                                     >  8:365406/1‑149 (MQ=255)
                                                                ccGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGc                                                                                    <  8:212292/148‑1 (MQ=255)
                                                                       cGGCACGCTGGTGCTCGCCGCTGCGGCCATGCGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACAcg                                                                            >  8:284008/1‑149 (MQ=255)
                                                                        ggCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCgccgcc                                                                                                                 >  6:165036/1‑111 (MQ=255)
                                                                        ggCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCgccgcc                                                                                                                 <  5:165036/111‑1 (MQ=255)
                                                                        ggCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACAcgc                                                                           >  8:110472/1‑149 (MQ=255)
                                                                                 gTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGaa                                                                   <  7:313079/148‑1 (MQ=255)
                                                                                  tGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACa                                                                                                                                                      >  4:143764/1‑64 (MQ=255)
                                                                                  tGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACa                                                                                                                                                      <  3:143764/64‑1 (MQ=255)
                                                                                    cTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTcc                                                                                        >  2:44830/1‑124 (MQ=255)
                                                                                    cTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTcc                                                                                        <  1:44830/124‑1 (MQ=255)
                                                                                     tCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGTCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTccc                                                              >  5:182882/1‑149 (MQ=255)
                                                                                       gccgcTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGa                                                            >  1:168986/1‑149 (MQ=255)
                                                                                           cTGCGGCCATGGGTACGCGGGAGGGGGCCCAGGCCCGGCGCC‑AGCGCGTTCAACAGCGGTTCCAGGCGCTCTAACACGTCCTTTTCGCCGCCTATCATCATGTAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTc                                                        <  1:90514/149‑1 (MQ=255)
                                                                                                gCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCTCGGTCCAGGCCCCACACGCCGCCgg                                                                     <  4:320206/131‑1 (MQ=255)
                                                                                                gCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCgg                                                                     >  3:320206/1‑131 (MQ=255)
                                                                                                   aTGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAAt                                                >  3:209667/1‑149 (MQ=255)
                                                                                                      ggTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGGCGCCGGAAGTCCCGACGTCCAGATAATGc                                              <  7:284008/148‑1 (MQ=255)
                                                                                                       gTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAg                                            >  2:301056/1‑149 (MQ=255)
                                                                                                        tACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAgg                                           <  4:209667/149‑1 (MQ=255)
                                                                                                            cgGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCa                                             >  2:214757/1‑143 (MQ=255)
                                                                                                            cgGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCa                                             <  1:214757/143‑1 (MQ=255)
                                                                                                                aTGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGTCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCtt                                   <  6:182882/149‑1 (MQ=255)
                                                                                                                  gTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCa                                                      >  3:188079/1‑128 (MQ=255)
                                                                                                                  gTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCa                                                      <  4:188079/128‑1 (MQ=255)
                                                                                                                     gCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACAcgc                                                                           >  4:203431/1‑104 (MQ=255)
                                                                                                                     gCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACAcgc                                                                           <  3:203431/104‑1 (MQ=255)
                                                                                                                     gCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGcc                               <  3:215913/148‑1 (MQ=255)
                                                                                                                     gCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGcc                               <  1:143035/148‑1 (MQ=255)
                                                                                                                     gCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGc                                >  3:265067/1‑147 (MQ=255)
                                                                                                                      cccACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACACCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAg                             >  2:166490/1‑149 (MQ=255)
                                                                                                                        cACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACgccgcc                                                                       >  6:37515/1‑105 (MQ=255)
                                                                                                                        cACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACgccgcc                                                                       <  5:37515/105‑1 (MQ=255)
                                                                                                                           gCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTg                        >  7:231084/1‑149 (MQ=255)
                                                                                                                           gCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTg                        <  8:231084/149‑1 (MQ=255)
                                                                                                                             ccGGCGCC‑AGCGCCTTGAACCGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTg                        >  6:244050/1‑147 (MQ=255)
                                                                                                                             ccGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCa                      <  1:369599/149‑1 (MQ=255)
                                                                                                                                     gaGCGCCTGGAACCGCTGTTCAAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCgcggc               >  5:2071/2‑149 (MQ=255)
                                                                                                                                      aGCGCCTGGAACCGCTGTTCAAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCa                      <  3:78256/141‑1 (MQ=255)
                                                                                                                                         gcCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCAt           >  3:193453/1‑149 (MQ=255)
                                                                                                                                           cTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTg                        <  6:352917/134‑1 (MQ=255)
                                                                                                                                           cTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTg                        >  5:352917/1‑134 (MQ=255)
                                                                                                                                             tGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATgtcg       >  2:331086/1‑149 (MQ=255)
                                                                                                                                               aaCAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATgtcgtc     >  6:142602/1‑149 (MQ=255)
                                                                                                                                                aCAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATGTCGTcc    <  4:265067/149‑1 (MQ=255)
                                                                                                                                                  aGCGGTTCCAGGCGCTAGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGAGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCTCGCTAGGACAGTTCTGCAGCGCGGCGCATGTCGTCCtt  >  7:329267/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTCGATGCCGTTATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCC‑AGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATGTCGTCCTT  >  NC_002947/4556889‑4557182

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