Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,388,502 G→C 14.3% A290G (GCC→GGC)  mscM ← mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,388,5020GC14.3% 193.7 / 23.5 63A290G (GCC→GGC) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base C (4/5);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

CAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAAT  >  NC_000913/4388362‑4388636
                                                                                                                                            |                                                                                                                                      
cAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAGCGTCAGGCTGCCAGccag                                                                                                                                 >  1:578788/1‑146 (MQ=255)
 aaGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCAGGCTGCCAGCCAGa                                                                                                                                <  2:502440/148‑1 (MQ=255)
      gttgtGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAGCGTCAGGCTGCCAGCCAGACgtta                                                                                                                           <  2:578788/148‑4 (MQ=255)
           ggTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGt                                                                                                                          <  2:365055/144‑1 (MQ=255)
           ggTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGt                                                                                                                          >  1:365055/1‑144 (MQ=255)
                           gCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCAGGCTGCCAGccag                                                                                                                                 >  2:368413/1‑119 (MQ=255)
                           gCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCAGGCTGCCAGccag                                                                                                                                 <  1:368413/121‑3 (MQ=255)
                                  tGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTgc                                                                                               >  2:328914/1‑148 (MQ=255)
                                     cACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGcc                                                                                            <  3:144739/148‑1 (MQ=255)
                                         tGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatccat                                                                                                      <  3:8744/134‑1 (MQ=255)
                                         tGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatccat                                                                                                      >  4:8744/1‑134 (MQ=255)
                                               cGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGa                                                                                                                                <  3:213776/102‑1 (MQ=255)
                                               cGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGa                                                                                                                                >  4:213776/1‑102 (MQ=255)
                                                    cgcTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCAGGCTGCCAGCCAGa                                                                                                                                <  4:164073/97‑1 (MQ=255)
                                                    cgcTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCAGGCTGCCAGCCAGa                                                                                                                                >  3:164073/1‑97 (MQ=255)
                                                            cGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAg                                                                     >  2:729300/1‑148 (MQ=255)
                                                              aGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTg                                                                                                                                 >  2:1013556/1‑86 (MQ=255)
                                                              aGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTg                                                                                                                                 <  1:1013556/86‑1 (MQ=255)
                                                              aGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatccat                                                                                                      <  2:320368/113‑1 (MQ=255)
                                                              aGCAGATTGGACGATCCCAGCCATTGCGACTGGTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatccat                                                                                                      >  1:320368/1‑113 (MQ=255)
                                                                    ttGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGt                                                                                                                          <  1:8815/87‑1 (MQ=255)
                                                                    ttGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGt                                                                                                                          >  2:8815/1‑87 (MQ=255)
                                                                    ttGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCg                                                                                                                      >  2:200048/1‑91 (MQ=255)
                                                                    ttGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCg                                                                                                                      <  1:200048/91‑1 (MQ=255)
                                                                    ttGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGt                                                             >  4:773806/1‑148 (MQ=255)
                                                                      ggACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACgag                                                                                                            <  4:480727/99‑1 (MQ=255)
                                                                      ggACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACgag                                                                                                            >  3:480727/1‑99 (MQ=255)
                                                                           aTCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCaa                                                                                                                >  2:1075000/1‑90 (MQ=255)
                                                                           aTCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCaa                                                                                                                <  1:1075000/90‑1 (MQ=255)
                                                                           aTCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAg                                                                               >  1:822362/1‑123 (MQ=255)
                                                                           aTCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAg                                                                               <  2:822362/123‑1 (MQ=255)
                                                                             cccAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCaa                                                                                 >  2:446196/1‑119 (MQ=255)
                                                                             cccAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCaa                                                                                 <  1:446196/119‑1 (MQ=255)
                                                                               cAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCAGGCTGCCAGccaga                                                                                                                                >  3:192004/1‑67 (MQ=255)
                                                                               cAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCAGGCTGCCAGCCAGa                                                                                                                                <  4:192004/70‑1 (MQ=255)
                                                                                aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACgaga                                                                                                           >  2:61984/1‑90 (MQ=255)
                                                                                aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACgaga                                                                                                           <  4:795593/90‑1 (MQ=255)
                                                                                aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACgaga                                                                                                           >  4:489190/1‑90 (MQ=255)
                                                                                aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACgaga                                                                                                           <  1:53892/90‑1 (MQ=255)
                                                                                aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACgaga                                                                                                           <  1:61984/90‑1 (MQ=255)
                                                                                aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACgaga                                                                                                           <  3:489190/90‑1 (MQ=255)
                                                                                aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACgaga                                                                                                           >  2:53892/1‑90 (MQ=255)
                                                                                aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACgaga                                                                                                           >  3:795593/1‑90 (MQ=255)
                                                                                aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAAtt                                                 >  2:254114/1‑148 (MQ=255)
                                                                                 gCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTg                                                <  1:729300/148‑1 (MQ=255)
                                                                                  ccATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGc                                                                                                  >  4:89816/1‑97 (MQ=255)
                                                                                  ccATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGc                                                                                                  <  3:89816/97‑1 (MQ=255)
                                                                                      tGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGcga                                           <  3:773806/148‑1 (MQ=255)
                                                                                               tCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGccgccg                                                                         <  2:75184/109‑1 (MQ=255)
                                                                                               tCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGccgccg                                                                         >  1:75184/1‑109 (MQ=255)
                                                                                                           ttCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCgg                      <  1:328914/148‑1 (MQ=255)
                                                                                                               aCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTg                                                                                       <  3:618612/79‑1 (MQ=255)
                                                                                                               aCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTg                                                                                       >  4:618612/1‑79 (MQ=255)
                                                                                                               aCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGt                                                                    >  1:1054092/1‑98 (MQ=255)
                                                                                                               aCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGt                                                                    <  2:1054092/98‑1 (MQ=255)
                                                                                                                cGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTgc                                               >  1:241549/1‑118 (MQ=255)
                                                                                                                cGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTgc                                               <  2:241549/118‑1 (MQ=255)
                                                                                                                             tGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGAtt                                                                                    >  2:35825/1‑68 (MQ=255)
                                                                                                                             tGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGAtt                                                                                    <  1:35825/68‑1 (MQ=255)
                                                                                                                                aaCGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCttt                                 >  3:342589/1‑116 (MQ=255)
                                                                                                                                aaCGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCttt                                 <  4:342589/116‑1 (MQ=255)
                                                                                                                                  cGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCgg                      >  4:692989/1‑125 (MQ=255)
                                                                                                                                  cGTCTGGCTGGCAGCCTGACGCTGTCGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCgg                      <  3:692989/125‑1 (MQ=255)
                                                                                                                                    tcaggctgccAGCCAGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccg                                                                            <  1:66800/66‑1 (MQ=38)
                                                                                                                                    tcaggctgccAGCCAGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccg                                                                            >  2:66800/4‑69 (MQ=38)
                                                                                                                                     ctggctggCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAAt  <  4:546375/142‑1 (MQ=255)
                                                                                                                                     ctggctggCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAAt  >  3:546375/1‑142 (MQ=255)
                                                                                                                                       ggctggcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAAttt                                                       >  1:872044/1‑87 (MQ=255)
                                                                                                                                       ggctggcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAAttt                                                       <  2:872044/87‑1 (MQ=255)
                                                                                                                                        gctggcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTg                                                >  2:20289/1‑93 (MQ=255)
                                                                                                                                        gctggcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATACGCGGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTg                                                <  1:20289/93‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                      
CAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAAT  >  NC_000913/4388362‑4388636

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