Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 2,946,944 G→A 6.3% D241D (GAC→GAT amiC ← N‑acetylmuramoyl‑L‑alanine amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090482,946,9440GA 6.3% 515.4 / 15.5 160D241D (GAC→GATamiCN‑acetylmuramoyl‑L‑alanine amidase
Reads supporting (aligned to +/- strand):  ref base G (78/71);  new base A (5/5);  total (83/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 167 total aligned reads displayed.

GATAGAGACAAACAGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCGAGC  >  AP009048/2946875‑2947012
                                                                     |                                                                    
gATAGAGACAAACAGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGt                                                                     >  1:771815/1‑71 (MQ=255)
 aTAGAGACAAACAGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATg                                                                      >  1:4687213/1‑69 (MQ=255)
   agagACAAACAGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGt                                                                     >  1:5030855/1‑68 (MQ=255)
      gaCAAACAGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCAt                                                               <  1:745242/71‑1 (MQ=255)
           aCAGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTgcgc                                                          >  1:1339633/1‑71 (MQ=255)
           aCAGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTgcgc                                                          <  1:5463728/71‑1 (MQ=255)
             aGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGt                                                        >  1:6329864/1‑71 (MQ=255)
               gTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCa                                                      <  1:4215510/71‑1 (MQ=255)
               gTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCa                                                      <  1:885085/71‑1 (MQ=255)
               gTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATCTCTTCAtt                                                              <  1:3678484/63‑1 (MQ=255)
                tCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCAtt                                                              <  1:4244236/62‑1 (MQ=255)
                tCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGcgcg                                                         <  1:2120383/67‑1 (MQ=255)
                  aGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGTAATGAAGATGTCTTCATTGCGCGTCATg                                                    <  1:5126200/70‑1 (MQ=255)
                  aGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATg                                                    <  1:5874364/70‑1 (MQ=255)
                  aGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATg                                                    >  1:1655163/1‑70 (MQ=255)
                        cTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCt                                             <  1:3239745/71‑1 (MQ=255)
                        cTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCt                                             <  1:7806875/71‑1 (MQ=255)
                        cTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCt                                             <  1:5154877/71‑1 (MQ=255)
                         tGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCtt                                            >  1:6251323/1‑71 (MQ=255)
                          gTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCtt                                            >  1:7705240/1‑70 (MQ=255)
                           tttCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTc                                           <  1:1788476/70‑1 (MQ=255)
                            ttCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCtt                                                                  >  1:7066208/1‑46 (MQ=255)
                            ttCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCa                                                      >  1:3085061/1‑58 (MQ=255)
                             tCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATg                                                                      <  1:3556997/41‑1 (MQ=255)
                             tCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATg                                                                      <  1:296289/41‑1 (MQ=255)
                              cTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCata                                        <  1:4810583/70‑1 (MQ=255)
                               tGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATg                                                    <  1:915798/57‑1 (MQ=255)
                               tGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTa                                                  <  1:6814461/59‑1 (MQ=255)
                                gggCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTg                                     >  1:2818235/1‑71 (MQ=255)
                                 ggCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTg                                     >  1:4780500/1‑70 (MQ=255)
                                 ggCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTg                                     <  1:393954/70‑1 (MQ=255)
                                 ggCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTg                                     <  1:1390474/70‑1 (MQ=255)
                                    ttttGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCat                                         <  1:2111745/63‑1 (MQ=255)
                                    ttttGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCat                                         <  1:1946143/63‑1 (MQ=255)
                                     tttGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCctc                                >  1:7585155/1‑71 (MQ=255)
                                     tttGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCctc                                >  1:925210/1‑71 (MQ=255)
                                      ttGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCctct                               >  1:1785090/1‑71 (MQ=255)
                                       tGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTgcgc                                                          <  1:1388639/43‑1 (MQ=255)
                                       tGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCa                                                      >  1:5174684/1‑47 (MQ=255)
                                       tGCTACGAGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCa                                                      >  1:7525904/1‑47 (MQ=255)
                                         cTACGCGCACTTGGAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCttt                             <  1:242724/70‑1 (MQ=255)
                                         cTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCttt                             <  1:2922658/70‑1 (MQ=255)
                                         cTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATAGCCCTCttt                             <  1:6168494/70‑1 (MQ=255)
                                           aCGCGCACTTGCAATGGAATGAAGATGTCTTCATTgcgc                                                          >  1:5600683/1‑39 (MQ=255)
                                           aCGCGCACTTGCAATGGAATGAAGATGTCTTCATTgcgc                                                          >  1:1876456/1‑39 (MQ=255)
                                           aCGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGc                                    <  1:5094792/61‑1 (MQ=255)
                                             gcgcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGa                        <  1:4703887/71‑1 (MQ=255)
                                             gcgcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGa                        <  1:4031885/71‑1 (MQ=255)
                                              cgcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGAt                       <  1:3234223/71‑1 (MQ=255)
                                               gcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTtctc                          <  1:1877889/67‑1 (MQ=255)
                                               gcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATc                      <  1:6532932/71‑1 (MQ=255)
                                                cACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCata                                        >  1:1394959/1‑52 (MQ=255)
                                                cACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCat                                         >  1:7358818/1‑51 (MQ=255)
                                                cACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCa                     <  1:6507623/71‑1 (MQ=255)
                                                 aCTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGccc                                  >  1:3457132/1‑57 (MQ=255)
                                                 aCTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGccc                                  >  1:2698513/1‑57 (MQ=255)
                                                 aCTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGAt                       >  1:3405241/1‑68 (MQ=255)
                                                 aCTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGAt                       >  1:5517602/1‑68 (MQ=255)
                                                  cTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGt                                                   >  1:1650113/1‑39 (MQ=255)
                                                  cTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCa                     <  1:6222746/69‑1 (MQ=255)
                                                   ttGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCa                     >  1:7565462/1‑68 (MQ=255)
                                                   ttGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCa                     >  1:3192337/1‑68 (MQ=255)
                                                    tGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCat                                         <  1:3136256/47‑1 (MQ=255)
                                                    tGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCat                                         <  1:2443690/47‑1 (MQ=255)
                                                    tGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCaga                   >  1:4441014/1‑69 (MQ=255)
                                                     gcaatgcaatGAAGATGTCTTCATTTCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagaga                 <  1:7372228/70‑1 (MQ=255)
                                                     gCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagagag                >  1:3876816/1‑71 (MQ=255)
                                                     gCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagag                  >  1:1441822/1‑69 (MQ=255)
                                                      cAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAgc               <  1:2913797/71‑1 (MQ=255)
                                                         tGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCatat                                       <  1:4643698/44‑1 (MQ=255)
                                                         tGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAgcgc             >  1:2323170/1‑70 (MQ=255)
                                                          ggAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCatat                                       <  1:6907926/43‑1 (MQ=255)
                                                          ggAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCa            >  1:4065604/1‑70 (MQ=255)
                                                          gcaaTGAAGACATCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTtct                           <  1:2768879/53‑1 (MQ=38)
                                                           gAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGAt                       >  1:397980/1‑58 (MQ=255)
                                                           gAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCaaa          >  1:5861192/1‑71 (MQ=255)
                                                           gAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCaaa          <  1:5259272/71‑1 (MQ=255)
                                                           gAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCaaa          <  1:4959796/71‑1 (MQ=255)
                                                           gAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCaa           <  1:7628913/70‑1 (MQ=255)
                                                            aaTGAAGACATCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCg                         >  1:1211694/1‑55 (MQ=39)
                                                             aTGAAGATGTCTTCATTGCGCGTCATGTACACCTTCatat                                       >  1:2645161/1‑40 (MQ=255)
                                                             aTGAAGACATCTTCATTGCGCGTCATGTACACCTTCATAtt                                      <  1:7499178/41‑1 (MQ=37)
                                                             aTGAAGACATCTTCATTGCGCGTCATGTACACCTTCATAtt                                      <  1:5405705/41‑1 (MQ=37)
                                                             aTGAAGACATCTTCATTGCGCGTCATGTACACCTTCATAtt                                      <  1:180057/41‑1 (MQ=37)
                                                              tGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagaga                 >  1:4092638/1‑61 (MQ=255)
                                                              tGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagaga                 >  1:1438045/1‑61 (MQ=255)
                                                              tGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagaga                 >  1:5958207/1‑61 (MQ=255)
                                                               gAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagaga                 >  1:7464936/1‑60 (MQ=255)
                                                               gAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagaga                 >  1:3661578/1‑60 (MQ=255)
                                                               gAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagaga                 >  1:2538192/1‑60 (MQ=255)
                                                               gAAGACATCTTCATTGCGCGTCATGTACACCTTCATATTg                                     >  1:1743867/1‑40 (MQ=37)
                                                               gAAGACATCTTCATTGCGCGTCATGTACACCTTCATATTGc                                    >  1:940505/1‑41 (MQ=38)
                                                               gAAGACATCTTCATTGCGCGTCATGTACACCTTCATATTGc                                    >  1:2844653/1‑41 (MQ=38)
                                                                aaGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGc      <  1:7077391/70‑1 (MQ=255)
                                                                aaGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGc      <  1:6482090/70‑1 (MQ=255)
                                                                aaGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGc      <  1:3848191/70‑1 (MQ=255)
                                                                   aTGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagaga                 <  1:3375845/56‑1 (MQ=255)
                                                                    tGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCGAGc  >  1:1559268/1‑70 (MQ=255)
                                                                     gTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCaa           >  1:2663677/1‑60 (MQ=255)
                                                                     gTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCaa           >  1:1233553/1‑60 (MQ=255)
                                                                     |                                                                    
GATAGAGACAAACAGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCGAGC  >  AP009048/2946875‑2947012

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