Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,909,328 G→A 9.2% Q385* (CAG→TAG) ‡ aprDA ← alkaline protease secretion ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,909,3280GA 9.2% 149.3 / 4.5 54Q385* (CAG→TAG) ‡aprDAalkaline protease secretion ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (24/25);  new base A (2/3);  total (26/28)
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

GCTGTCCACCTCGGCAAAGCGGGCGATGTTTTCGGCGATGGTGCCTTCGAACAGCTCCACGTCCTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTGTCACCCGGGGCCAGGCTGAAGCTCACCCCGCGTACGATGCTGGTGCTGCTGCCTGGGGCG  >  NC_002947/2909179‑2909464
                                                                                                                                                     |                                                                                                                                         
gCTGTCCACCTCGGCAAAGCGGGCGATGTTTTCGGCGATGGTGCCTTCGAACAGCTCCACGTCCTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTg                                                                                                                                           >  2:1154515/1‑150 (MQ=255)
 cTGTCCACCTCGGCAAAGCGGGCGATGTTTTCGGCGATGGTGCCTTCGAACAGCTCCACGTCCTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTgg                                                                                                                                          <  2:453485/150‑1 (MQ=255)
       aCCTCGGCAAAGCGGGCGATGTTTTCGGCGATGGTGCCTTCGAACAGCTCCACGTCCTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGggccggc                                                                                                                                    <  1:1272916/150‑1 (MQ=255)
       aCCTCGGCAAAGCGGGCGATGTTTTCGGCGATGGTGCCTTCGAACAGCTCCACGTCCTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGggccggc                                                                                                                                    >  1:688960/1‑150 (MQ=255)
            ggCAAAGCGGGCGATGTTTTCGGCGATGGTGCCTTCGAACAGCTCCACGTCCTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCCAggccggcaaggt                                                                                                                               <  2:592946/150‑6 (MQ=255)
            ggCAAAGCGGGCGATGTTTTCGGCGATGGTGCCTTCGAACAGCTCCACGTCCTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCCAggccggcaag                                                                                                                                 >  2:483484/1‑145 (MQ=255)
            ggCAAAGCGGGCGATGTTTTCGGCGATGGTGCCTTCGAACAGCTCCACGTCCTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCCAggccggcaag                                                                                                                                 <  1:483484/148‑4 (MQ=255)
              cAAAGCGGGCGAGGTTTTCGGCGAGGGGGCCTTCGAACAGCGCCACGGCCTGCGGCAAGTAGCCCAGCCACGGGCCCGGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGc                                                                                                                             <  1:160842/150‑1 (MQ=255)
                     ggCGATGTTTTCGGCGATGGTGCCTTCGAACAGCTCCACGTCCTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCNCCGTCCNGGCGCACCTTGCNGGCCTGGGCCGGCCACACGCCTACCaa                                                                                                                      <  1:1154515/150‑1 (MQ=255)
                            ttttCGGCGATGGTGCCTTCGAACAGCTCCACGTCCTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAggcgg                                                                                                               <  1:1066655/150‑1 (MQ=255)
                            ttttCGGCGATGGTGCCTTCGAACAGCTCCACGTCCTGCGGCAGGTAGCCCAGCCACGAGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAggcgg                                                                                                               >  1:880731/1‑150 (MQ=255)
                                 ggCGATGGTGCCTTCGAACAGCTCCACGTCCTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAg                                                                                                          <  1:758758/150‑1 (MQ=255)
                                            cTTCGAACAGCTCCACGTCCTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCCAggccggcaag                                                                                                                                 <  2:1036825/116‑4 (MQ=255)
                                            cTTCGAACAGCTCCACGTCCTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCCAggccggcaag                                                                                                                                 >  1:1036825/1‑113 (MQ=255)
                                                gAACAGCTCCACGTCCTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGNGCCGGCNACACGCCTACCANCAGGCGGGCGAGGGTAGACTTGCccga                                                                                           <  1:860754/150‑1 (MQ=255)
                                                      cTCCACGTCCTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCccgaggccga                                                                                    <  2:1306906/150‑1 (MQ=255)
                                                      cTCCACGTCCTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCccgaggccga                                                                                    >  1:711657/1‑150 (MQ=255)
                                                           cGTCCTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGc                                                                               <  1:1297216/150‑1 (MQ=255)
                                                              ccTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGa                                                                            <  2:571474/150‑1 (MQ=255)
                                                                         tAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGAACGAGGCCGAGGGGCCGATGACACCCAgg                                                                 <  1:336389/150‑1 (MQ=255)
                                                                          aGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACa                                                                                                                    >  1:759311/1‑99 (MQ=255)
                                                                          aGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACa                                                                                                                    <  2:759311/99‑1 (MQ=255)
                                                                           gCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGGCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCt                                                               <  1:1245279/150‑1 (MQ=255)
                                                                            cccGGCCTCTGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCtt                                                              >  2:55787/1‑149 (MQ=255)
                                                                            cccAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGAAGGCACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTg                                                              >  1:601950/1‑150 (MQ=255)
                                                                             ccAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAgg                                                                                                         >  1:486440/1‑107 (MQ=255)
                                                                             ccAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAgg                                                                                                         <  2:486440/107‑1 (MQ=255)
                                                                             ccAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGCCAGCACCCGATAATACCAACATGCGGGGGGGGGTAGACTTCTCGCTGCCGGAGCGGGCGAGGACACCCATGCTGt                                                             >  1:26370/1‑150 (MQ=255)
                                                                              cAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTGTc                                                            >  1:270647/1‑150 (MQ=255)
                                                                               aGCCACGAGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTGTCa                                                           <  2:880731/150‑1 (MQ=255)
                                                                                       gCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTGTCACCCGGGGc                                                   >  2:87187/1‑150 (MQ=255)
                                                                                       gCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTGTCACCCGGGGc                                                   >  1:1405479/1‑150 (MQ=255)
                                                                                        cccAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTa                                                                                                                          >  1:572721/1‑79 (MQ=255)
                                                                                        cccAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTa                                                                                                                          <  2:572721/79‑1 (MQ=255)
                                                                                         ccAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATg                                                                          >  2:1354512/1‑125 (MQ=255)
                                                                                         ccAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATg                                                                          <  1:1354512/125‑1 (MQ=255)
                                                                                         ccAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTGTCACCCGGGGCCa                                                 <  2:802397/150‑1 (MQ=255)
                                                                                         ccAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTGTCACCCGGGGCCa                                                 <  2:711657/150‑1 (MQ=255)
                                                                                           aGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAg                                                                  >  1:233577/1‑131 (MQ=255)
                                                                                           aGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAg                                                                  <  2:233577/131‑1 (MQ=255)
                                                                                                ggCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGAt                                                                           <  1:39283/117‑1 (MQ=255)
                                                                                                ggCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGAt                                                                           >  2:39283/1‑117 (MQ=255)
                                                                                                  ccTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTATCACCCGGGGCCAGGCTGa                                           <  2:160996/147‑1 (MQ=255)
                                                                                                  ccTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTATCACCCGGGGCCAGGCTGa                                           >  1:160996/1‑147 (MQ=255)
                                                                                                    ttGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTGTCACCCGGGGCCAGGCTGAAGCTc                                      >  2:108056/1‑150 (MQ=255)
                                                                                                        tCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACCCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGCGCGGCCGATGACACCCAGCCTCTCACCCCGGGCCCGGCTCAAGCCCAccc                                  >  2:1392484/1‑150 (MQ=255)
                                                                                                              tGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTGTCACCCGGGGCCAGGCTGAAGCTCACCCCGCGTa                            >  1:372247/1‑150 (MQ=255)
                                                                                                                           gcCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCa                                                                   <  2:273687/98‑1 (MQ=255)
                                                                                                                           gcCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCa                                                                   >  1:273687/1‑98 (MQ=255)
                                                                                                                              gTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTGTCACCCGGGGCCAGGCTGAAGCTCACCCCGCGTACGATGCTGGtgctgct            >  1:28381/1‑150 (MQ=255)
                                                                                                                                  aGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTGTCACCCGGGGCCAGGCTGAAGCTCACCCCGCGTACGATGCTGGTGCTGCTGCCt        <  1:108056/150‑1 (MQ=255)
                                                                                                                                    gcgcACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTGTCACCCGGGGCCAGGCTGAAGCTCACCCCGCGTACGATGCTGGTGCTGCTGCCTgg      >  2:214703/1‑150 (MQ=255)
                                                                                                                                    gcgcACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTGTCACCCGGGGCCAGGCTGAAGCTCACCCCGCGTACGATGCTGGTGCTGCTGCCTgg      >  2:1354084/1‑150 (MQ=255)
                                                                                                                                     cgcACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTGTCACCCGGGGCCAGGCTGAAGCTCACCCCGCGTACGATGCTGGTGCTGCTGCCTggg     <  1:87187/150‑1 (MQ=255)
                                                                                                                                     cgcACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTGTCACCCGGGGCCAGGCTGAAGCTCACCCCGCGTACGATGCTGGTGCTGCTGCCTggg     <  2:1286426/150‑1 (MQ=255)
                                                                                                                                        aCCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTGTCACCCGGGGCCAGGCTGAAGCTCACCCCGCGTACGATGCTGGTGCTGCTGCCTGGGgcg  >  1:631629/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                         
GCTGTCCACCTCGGCAAAGCGGGCGATGTTTTCGGCGATGGTGCCTTCGAACAGCTCCACGTCCTGCGGCAGGTAGCCCAGCCACGGGCCCAGTTCGGCCTTGTTCCAGGTGAAGATATCGGCGCCGTCCAGGCGCACCTTGCCGGCCTGGGCCGGCCACACGCCTACCAACAGGCGGGCGAGGGTAGACTTGCCCGAGGCCGAGGGGCCGATGACACCCAGGCTGTCACCCGGGGCCAGGCTGAAGCTCACCCCGCGTACGATGCTGGTGCTGCTGCCTGGGGCG  >  NC_002947/2909179‑2909464

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: