Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,386,978 C→T W798* (TGG→TAG)  mscM ← miniconductance mechanosensitive channel MscM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,386,9780CT98.2% 196.8 / ‑6.5 57W798* (TGG→TAG) mscMminiconductance mechanosensitive channel MscM
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (29/27);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGC  >  NC_000913/4386842‑4387115
                                                                                                                                        |                                                                                                                                         
gatgatAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACtaca                                                                                                                                        <  1:827715/140‑1 (MQ=255)
    atAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGcac                                                                                                                                    <  1:182859/140‑1 (MQ=255)
      aaaCACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACaa                                                                                                                                  <  1:731736/140‑1 (MQ=255)
                  ggCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCa                                                                                                                      <  2:888738/140‑1 (MQ=255)
                    cAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGca                                                                                                                                     <  1:870719/123‑1 (MQ=255)
                    cAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGca                                                                                                                                     >  2:870719/1‑123 (MQ=255)
                      aTCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGc                                                                                                                  >  1:149120/1‑140 (MQ=255)
                                 accgAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCa                                                                                                       <  1:952448/140‑1 (MQ=255)
                                     aGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAAt                                                                                                   >  2:172463/1‑140 (MQ=255)
                                     aGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAAt                                                                                                   >  2:483035/1‑140 (MQ=255)
                                        gTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGa                                                                                                >  1:672503/1‑140 (MQ=255)
                                            ttGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAgcgc                                                                                            >  2:805308/1‑140 (MQ=255)
                                              ggCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCAc                                                                                          <  1:483035/140‑1 (MQ=255)
                                               gCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCAcc                                                                                         <  2:54876/140‑1 (MQ=255)
                                                cTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCa                                                                                        <  1:172463/140‑1 (MQ=255)
                                                cTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCa                                                                                        >  1:71155/1‑140 (MQ=255)
                                                 tCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAg                                                                                       >  1:456117/1‑140 (MQ=255)
                                                  ccAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGc                                                                                      >  2:600759/1‑140 (MQ=255)
                                                   cAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGcc                                                                                     >  2:770273/1‑140 (MQ=255)
                                                    aGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCg                                                                                    >  1:476275/1‑140 (MQ=255)
                                                       cTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCaa                                                                                 <  1:805308/140‑1 (MQ=255)
                                                             aCGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTg                                                                           <  2:456117/140‑1 (MQ=255)
                                                               gCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTgcg                                                                         <  1:600759/140‑1 (MQ=255)
                                                                    gTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGa                                                                              >  1:805163/1‑130 (MQ=255)
                                                                    gTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGa                                                                              <  2:805163/130‑1 (MQ=255)
                                                                    gTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTg                                                                    <  1:436532/140‑1 (MQ=255)
                                                                     tACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGc                                                                                   >  1:446374/1‑124 (MQ=255)
                                                                     tACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGc                                                                                   <  2:446374/124‑1 (MQ=255)
                                                                     tACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGa                                                                   >  2:370641/1‑140 (MQ=255)
                                                                         gTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGAATTGCGCACTGATGGc                                                               >  1:14967/1‑140 (MQ=255)
                                                                          tGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCa                                                              >  1:242923/1‑140 (MQ=255)
                                                                          tGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCa                                                              >  2:182680/1‑140 (MQ=255)
                                                                           ggAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCAt                                                             >  2:960544/1‑140 (MQ=255)
                                                                           ggAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCAt                                                             >  2:213450/1‑140 (MQ=255)
                                                                           ggAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCAt                                                             <  2:149120/140‑1 (MQ=255)
                                                                            gAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATc                                                            >  2:106336/1‑140 (MQ=255)
                                                                             aGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCg                                                           >  2:345661/1‑140 (MQ=255)
                                                                               gTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCgag                                                         >  2:613767/1‑140 (MQ=255)
                                                                                  aCATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATc                                                      >  1:726254/1‑140 (MQ=255)
                                                                                         acaGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCg                                               >  2:659806/1‑140 (MQ=255)
                                                                                               aaaTATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCa                                         <  1:370641/140‑1 (MQ=255)
                                                                                                aaTATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCAt                                        >  2:408496/1‑140 (MQ=255)
                                                                                                aaTATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCAt                                        >  1:399427/1‑140 (MQ=255)
                                                                                                    ttttCGAGGAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGAc                                    <  2:476275/140‑1 (MQ=255)
                                                                                                            gAAGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTg                            <  2:672503/140‑1 (MQ=255)
                                                                                                             aaGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGc                           <  1:213450/140‑1 (MQ=255)
                                                                                                              aGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCt                          <  2:242923/140‑1 (MQ=255)
                                                                                                              aGCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCt                          <  1:454154/140‑1 (MQ=255)
                                                                                                               gCCGAAAGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTc                         >  2:761460/1‑140 (MQ=255)
                                                                                                                   aaaGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTc                     <  1:106336/140‑1 (MQ=255)
                                                                                                                     aGCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGAATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCg                   <  2:14967/140‑1 (MQ=255)
                                                                                                                      gCGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCgg                  <  1:770273/140‑1 (MQ=255)
                                                                                                                       cGGAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTACggg                 <  1:790380/140‑1 (MQ=255)
                                                                                                                        ggAATGGATTTCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGc                <  2:399427/140‑1 (MQ=255)
                                                                                                                                tttCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGtga        >  2:561739/1‑140 (MQ=255)
                                                                                                                                 ttCTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGtgat       >  1:363159/1‑140 (MQ=255)
                                                                                                                                   cTGACTACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGtgatga     <  1:182680/140‑1 (MQ=255)
                                                                                                                                      aCCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATgc  <  2:892800/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                         
GATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGC  >  NC_000913/4386842‑4387115

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