Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,825,963 A→T 5.1% V28D (GTC→GAC)  ptxS ← 2‑ketogluconate utilization repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,825,9630AT 5.1% 253.0 / 2.7 79V28D (GTC→GAC) ptxS2‑ketogluconate utilization repressor
Reads supporting (aligned to +/- strand):  ref base A (35/40);  new base T (2/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

GCATGCCGATCAGTCGCGTCTGGCCGCGCTTTAGGCCGCGAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAACCGGTTTCAATACAGCAAAAATGCGCGG  >  NC_002947/3825820‑3826108
                                                                                                                                               |                                                                                                                                                 
gCATGCCGATCAGTCGCGTCTGGCCGCGCTTTAGGCCGCGAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTg                                                                                                                                              <  4:145855/149‑1 (MQ=255)
  aTGCCGATCAGTCGCGTCTGGCCGCGCTTTAGGCCGCGAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGc                                                                                                                                            >  1:166953/1‑149 (MQ=255)
     ccGATCAGTCGCGTCTGGCCGCGCTTTAGGCCGCGAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCttt                                                                                                                                         <  1:223508/149‑1 (MQ=255)
      cGATCAGTCGCGTCTGGCCGCGCTTTAGGCCGCGAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTg                                                                                                                                        <  5:88235/149‑1 (MQ=255)
           aGTCGCGTCTGGCCGCGCTTTAGGCCGCGAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAAc                                                                                                                                   <  2:166953/149‑1 (MQ=255)
            gTCGCGTCTGGCCGCGCTTTAGGCCGCGAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACg                                                                                                                                  >  3:226197/1‑149 (MQ=255)
             tCGCGTCTGGCCGCGCTTTAGGCCGCGAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGc                                                                                                                                 >  2:109607/1‑149 (MQ=255)
               gcgTCTGGCCGCGCTTTAGGCCGCGAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCATGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCa                                                                                                                               <  6:113488/149‑1 (MQ=255)
                    tGGCCGCGCTTTAGGCCGCGAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCaa                                                                                                                           >  6:45577/1‑148 (MQ=255)
                        cgcgCTTTAGGCCGCGAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACAcgc                                                                                                                      <  4:226197/149‑1 (MQ=255)
                        cgcgCTTTAGGCCGCGAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACAcgc                                                                                                                      >  5:82340/1‑149 (MQ=255)
                                  gCCGCGAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCAc                                                                                                            >  7:129967/1‑149 (MQ=255)
                                     gcgAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGa                                                                                                         >  3:385634/1‑149 (MQ=255)
                                       gAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATg                                                                                                       <  3:209707/149‑1 (MQ=255)
                                         gCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACACCAGCAACACGCGCCACTTCACTGATGGt                                                                                                     >  3:309411/1‑149 (MQ=255)
                                           cATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGa                                                                                                   <  8:202907/149‑1 (MQ=255)
                                           cATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGa                                                                                                   <  8:42037/149‑1 (MQ=255)
                                             tCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGTTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGAcc                                                                                                 >  5:198534/1‑149 (MQ=255)
                                             tCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGAcc                                                                                                 >  3:142992/1‑149 (MQ=255)
                                             tCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGAcc                                                                                                 >  3:43279/1‑149 (MQ=255)
                                                ggTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACg                                                                                                                                                 >  6:160016/1‑98 (MQ=255)
                                                ggTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGt                                                                                              >  8:332012/1‑149 (MQ=255)
                                                 gTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGtt                                                                                             >  6:177307/1‑149 (MQ=255)
                                                     gggCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGc                                                                                                                                 >  3:66073/1‑109 (MQ=255)
                                                     gggCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGc                                                                                                                                 <  4:66073/109‑1 (MQ=255)
                                                     gggCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTcgcg                                                                                         <  6:82340/149‑1 (MQ=255)
                                                     gggCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTcgc                                                                                          <  4:385634/148‑1 (MQ=255)
                                                        cGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGt                                                                                      <  5:177307/149‑1 (MQ=255)
                                                             gCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATgg                                                                                  <  1:109607/148‑1 (MQ=255)
                                                              ccGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATggg                                                                                 >  8:30410/1‑148 (MQ=255)
                                                                gAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGTCGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAg                                                                              >  4:181616/1‑149 (MQ=255)
                                                                gAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGTCGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAg                                                                              >  4:187881/1‑149 (MQ=255)
                                                                 aGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAgg                                                                             >  3:185495/1‑149 (MQ=255)
                                                                     gTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGcgcg                                                                         >  2:86123/1‑149 (MQ=255)
                                                                      ttCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCATCAATTGGCGGTCGCCACCGATGTAACGGGAGTCGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGt                                                                        <  3:187881/149‑1 (MQ=255)
                                                                      ttCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGTCGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGt                                                                        <  3:181616/149‑1 (MQ=255)
                                                                             aCCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTTGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCAt                                                                 <  1:135565/149‑1 (MQ=255)
                                                                              ccTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATg                                                                >  5:128475/1‑149 (MQ=255)
                                                                              ccTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATg                                                                <  2:305429/149‑1 (MQ=255)
                                                                              ccTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATg                                                                >  6:403350/1‑149 (MQ=255)
                                                                               ctctTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATgg                                                               >  5:64257/1‑149 (MQ=255)
                                                                               ctctTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATg                                                                <  7:305835/148‑1 (MQ=255)
                                                                               ctctTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATg                                                                >  8:305835/1‑148 (MQ=255)
                                                                                  ttCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGtt                                                            >  2:303417/1‑149 (MQ=255)
                                                                                     cAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTc                                                         <  1:86123/149‑1 (MQ=255)
                                                                                     cAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTc                                                         <  4:185495/149‑1 (MQ=255)
                                                                                     cAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACACCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGtttt                                                          <  4:309411/148‑1 (MQ=255)
                                                                                        ccGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCtt                                                      <  4:43279/149‑1 (MQ=255)
                                                                                          gCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATgg                                                                                                      >  4:99331/1‑99 (MQ=255)
                                                                                          gCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATgg                                                                                                      <  3:99331/99‑1 (MQ=255)
                                                                                          gCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTAt                                                    <  8:129967/149‑1 (MQ=255)
                                                                                          gCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTAt                                                    <  7:30410/149‑1 (MQ=255)
                                                                                          gCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTAt                                                    <  7:332012/149‑1 (MQ=255)
                                                                                                ccGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGt                                              >  3:66131/1‑149 (MQ=255)
                                                                                                     gCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTaa                                         <  5:404301/149‑1 (MQ=255)
                                                                                                       ttCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTa                                                     >  6:95984/1‑135 (MQ=255)
                                                                                                       ttCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTa                                                     <  5:95984/135‑1 (MQ=255)
                                                                                                       ttCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTaa                                         >  2:20818/1‑147 (MQ=255)
                                                                                                           gCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGgc                                                                            <  8:360394/108‑1 (MQ=255)
                                                                                                           gCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGgc                                                                            >  7:360394/1‑108 (MQ=255)
                                                                                                            ccAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGtgat                                  <  4:142992/149‑1 (MQ=255)
                                                                                                              aGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGtga                                   <  7:371862/146‑1 (MQ=255)
                                                                                                              aGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGtga                                   >  8:371862/1‑146 (MQ=255)
                                                                                                                 aaTTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAAc                             <  4:66131/149‑1 (MQ=255)
                                                                                                                  aTTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAAcc                            >  6:351138/1‑149 (MQ=255)
                                                                                                                         gTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAACCGGTTTc                      >  6:30790/1‑148 (MQ=255)
                                                                                                                          tCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAACCGGTTTCaa                    <  6:198534/149‑1 (MQ=255)
                                                                                                                          tCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAACCGGTTTCaa                    <  6:128475/149‑1 (MQ=255)
                                                                                                                            gCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAACCGGTTTCAAt                   >  2:207475/1‑148 (MQ=255)
                                                                                                                            gCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAACCGGTTTCAATa                  <  5:30790/149‑1 (MQ=255)
                                                                                                                             ccACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAACCGGTTTCAATAc                 <  1:207475/149‑1 (MQ=255)
                                                                                                                              cACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAAcc                            <  7:103122/137‑1 (MQ=255)
                                                                                                                              cACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAAcc                            >  8:103122/1‑137 (MQ=255)
                                                                                                                              cACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAACCGGTTTCAATACa                <  6:64257/149‑1 (MQ=255)
                                                                                                                                       aaCGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAACCGGTTTCAATACAGCAAAAATg       >  1:299309/1‑149 (MQ=255)
                                                                                                                                       aaCGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAACCGGTTTCAATACAGCAAAAATg       >  2:353386/1‑149 (MQ=255)
                                                                                                                                        aCGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTaa                                         <  3:302946/114‑1 (MQ=255)
                                                                                                                                          gggAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAACCGGTTTCAATACAGCAAAAATgcgc    <  2:299309/149‑1 (MQ=255)
                                                                                                                                          gggAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAACCGGTTTCAATACAGCAAAAATgcgc    <  5:403350/149‑1 (MQ=255)
                                                                                                                                            gagaCGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGTTTGTTTTAATAGTGATGAAACCGGTTTCAATACAGCAAAAATGCGCgg  <  1:303417/149‑1 (MQ=255)
                                                                                                                                            gagaCGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAACCGGTTTCAATACAGCAAAAATGCGCgg  >  3:314529/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GCATGCCGATCAGTCGCGTCTGGCCGCGCTTTAGGCCGCGAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAACCGGTTTCAATACAGCAAAAATGCGCGG  >  NC_002947/3825820‑3826108

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑