Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 2,946,943 T→C 6.2% I242V (ATC→GTC)  amiC ← N‑acetylmuramoyl‑L‑alanine amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090482,946,9430TC 6.2% 462.7 / 16.6 162I242V (ATC→GTC) amiCN‑acetylmuramoyl‑L‑alanine amidase
Reads supporting (aligned to +/- strand):  ref base T (78/74);  new base C (5/5);  total (83/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

Only 100 of 163 total aligned reads displayed.

GATAGAGACAAACAGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCGAGC  >  AP009048/2946875‑2947012
                                                                    |                                                                     
gATAGAGACAAACAGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGt                                                                     >  1:7070854/1‑71 (MQ=255)
 aTAGAGACAAACAGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATg                                                                      >  1:4687213/1‑69 (MQ=255)
  tAGAGACAAACAGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCt                                                                   >  1:363227/1‑71 (MQ=255)
   agagACAAACAGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGt                                                                     >  1:5030855/1‑68 (MQ=255)
    gagaCAAACAGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCtt                                                                  >  1:4113602/1‑70 (MQ=255)
             aGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGt                                                        >  1:396469/1‑71 (MQ=255)
             aGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGt                                                        >  1:6329864/1‑71 (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)
                tCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGcgcg                                                         <  1:2530530/67‑1 (MQ=255)
                 cAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTgc                                                            >  1:4731231/1‑63 (MQ=255)
                 cAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATg                                                    >  1:4019887/1‑71 (MQ=255)
                  aGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATg                                                    <  1:5874364/70‑1 (MQ=255)
                  aGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATg                                                    >  1:1655163/1‑70 (MQ=255)
                   gCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCa                                                                >  1:5504507/1‑57 (MQ=255)
                        cTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCt                                             <  1:7806875/71‑1 (MQ=255)
                        cTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCt                                             <  1:3239745/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:5868137/70‑1 (MQ=255)
                           tttCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTc                                           <  1:1788476/70‑1 (MQ=255)
                           tttCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCa                                          >  1:4688788/1‑71 (MQ=255)
                            ttCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCtt                                                                  >  1:924617/1‑46 (MQ=255)
                            ttCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCtt                                                                  >  1:7066208/1‑46 (MQ=255)
                            ttCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCtt                                                                  >  1:1372469/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)
                                 ggCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTg                                     <  1:393954/70‑1 (MQ=255)
                                 ggCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTg                                     >  1:4780500/1‑70 (MQ=255)
                                 ggCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGc                                    >  1:1393913/1‑71 (MQ=255)
                                    ttttGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCa                                                                >  1:2170854/1‑40 (MQ=255)
                                    ttttGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCat                                         <  1:1946143/63‑1 (MQ=255)
                                     tttGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCAtt                                                              <  1:1802385/41‑1 (MQ=255)
                                     tttGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCAtt                                                              <  1:4087380/41‑1 (MQ=255)
                                     tttGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCAtt                                                              <  1:7540263/41‑1 (MQ=255)
                                     tttGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCAtt                                                              <  1:2077256/41‑1 (MQ=255)
                                     tttGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATAtt                                      >  1:1419920/1‑65 (MQ=255)
                                     tttGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCctc                                >  1:925210/1‑71 (MQ=255)
                                       tGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTgcgc                                                          <  1:1388639/43‑1 (MQ=255)
                                         cTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCttt                             <  1:2922658/70‑1 (MQ=255)
                                         cTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATAGCCCTCttt                             <  1:6168494/70‑1 (MQ=255)
                                           aCGCGCACTTGCAATGGAATGAAGATGTCTTCATTgcgc                                                          >  1:5600683/1‑39 (MQ=255)
                                           aCGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGc                                    <  1:6834015/61‑1 (MQ=255)
                                              cgcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGAt                       <  1:3234223/71‑1 (MQ=255)
                                               gcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTtctc                          <  1:1877889/67‑1 (MQ=255)
                                               gcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATc                      <  1:6532932/71‑1 (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:777953/1‑68 (MQ=255)
                                                 aCTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGAt                       >  1:5517602/1‑68 (MQ=255)
                                                 aCTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGAt                       >  1:1512695/1‑68 (MQ=255)
                                                  cTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGt                                                   >  1:1650113/1‑39 (MQ=255)
                                                  cTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCa                     <  1:6222746/69‑1 (MQ=255)
                                                   ttGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCa                     >  1:7180538/1‑68 (MQ=255)
                                                   ttGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCa                     >  1:4487655/1‑68 (MQ=255)
                                                   ttGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCa                     >  1:3192337/1‑68 (MQ=255)
                                                   ttGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCa                     >  1:7565462/1‑68 (MQ=255)
                                                   ttGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCa                     >  1:5103479/1‑68 (MQ=255)
                                                    tGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCat                                         <  1:2443690/47‑1 (MQ=255)
                                                    tGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCat                                         <  1:3136256/47‑1 (MQ=255)
                                                     gCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagag                  <  1:1762078/69‑1 (MQ=255)
                                                      cAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAgc               <  1:2913797/71‑1 (MQ=255)
                                                         tGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCa            <  1:2570590/71‑1 (MQ=255)
                                                          ggAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCa            >  1:4065604/1‑70 (MQ=255)
                                                          gcaaTGAAGACATCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTtct                           <  1:2768879/53‑1 (MQ=38)
                                                           gAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCaaa          <  1:5259272/71‑1 (MQ=255)
                                                           gAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCaaa          <  1:5186724/71‑1 (MQ=255)
                                                           gAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCaa           <  1:7628913/70‑1 (MQ=255)
                                                            aaTGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAAc         <  1:7623789/71‑1 (MQ=255)
                                                            aaTGAAGACATCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCg                         >  1:1211694/1‑55 (MQ=39)
                                                             aTGAAGATGTCTTCATTGCGCGTCATGTACACCTTCatat                                       >  1:7438120/1‑40 (MQ=255)
                                                             aTGAAGACATCTTCATTGCGCGTCATGTACACCTTCATAtt                                      <  1:5405705/41‑1 (MQ=37)
                                                             aTGAAGACATCTTCATTGCGCGTCATGTACACCTTCATAtt                                      <  1:180057/41‑1 (MQ=37)
                                                              tGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagaga                 >  1:5958207/1‑61 (MQ=255)
                                                              tGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagaga                 >  1:1438045/1‑61 (MQ=255)
                                                              tGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACgg       >  1:7510257/1‑71 (MQ=255)
                                                              tGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACgg       >  1:4273651/1‑71 (MQ=255)
                                                              tGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACgg       <  1:2697772/71‑1 (MQ=255)
                                                               gAAGATGTCTTCATTGCGCGTCATGTACAGCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGc      >  1:2844595/1‑71 (MQ=255)
                                                               gAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagaga                 >  1:1258598/1‑60 (MQ=255)
                                                               gAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagaga                 >  1:7464936/1‑60 (MQ=255)
                                                               gAAGACATCTTCATTGCGCGTCATGTACACCTTCATATTg                                     >  1:1743867/1‑40 (MQ=37)
                                                               gAAGACATCTTCATTGCGCGTCATGTACACCTTCATATTGc                                    >  1:2844653/1‑41 (MQ=38)
                                                               gAAGACATCTTCATTGCGCGTCATGTACACCTTCATATTGc                                    >  1:940505/1‑41 (MQ=38)
                                                                aaGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAgcg              >  1:7455457/1‑62 (MQ=255)
                                                                aaGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCaaa          >  1:4913284/1‑66 (MQ=255)
                                                                aaGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGc      <  1:6482090/70‑1 (MQ=255)
                                                                aaGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGc      <  1:3848191/70‑1 (MQ=255)
                                                                aaGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGc      <  1:2410060/70‑1 (MQ=255)
                                                                aaGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCg     <  1:3855187/71‑1 (MQ=255)
                                                                  gATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCaaa          >  1:5004811/1‑64 (MQ=255)
                                                                   aTGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagaga                 <  1:5147243/56‑1 (MQ=255)
                                                                   aTGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagaga                 <  1:3392247/56‑1 (MQ=255)
                                                                   aTGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagaga                 <  1:3375845/56‑1 (MQ=255)
                                                                    |                                                                     
GATAGAGACAAACAGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCGAGC  >  AP009048/2946875‑2947012

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