Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 428,188 G→T R181L (CGC→CTC)  secD → SecYEG protein translocase auxillary subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913428,1880GT95.2% 206.8 / 2.7 62R181L (CGC→CTC) secDSecYEG protein translocase auxillary subunit
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (33/26);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

CAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCTGCGTGCGG  >  NC_000913/428091‑428283
                                                                                                 |                                                                                               
cAAAATATCGATAGCCTACGCAGTCACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCctcg                                                                                               <  2:488502/100‑1 (MQ=255)
   aaTATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATg                                                                                            >  2:1200711/1‑100 (MQ=255)
   aaTATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATg                                                                                            <  1:290995/100‑1 (MQ=255)
    atatCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGc                                                                                           >  1:981539/1‑100 (MQ=255)
    agatCGATAGCCTACGGAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATg                                                                                            <  2:427265/97‑1 (MQ=255)
     tatCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCt                                                                                          <  2:994292/100‑1 (MQ=255)
      atCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTa                                                                                         >  2:717184/1‑100 (MQ=255)
         gATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACATTCCTCGATGCTAAAg                                                                                      >  2:920571/1‑100 (MQ=255)
             gCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCg                                                                                  >  1:376358/1‑100 (MQ=255)
               cTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtg                                                                                <  1:1040479/100‑1 (MQ=255)
               cTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtg                                                                                <  2:703668/100‑1 (MQ=255)
                tACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtga                                                                               <  1:570269/100‑1 (MQ=255)
                  cGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtgatg                                                                             >  1:92956/1‑100 (MQ=255)
                    cAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGaa                                                                           >  1:1171254/1‑100 (MQ=255)
                    cAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGaa                                                                           >  2:391244/1‑100 (MQ=255)
                      gTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGc                                                                         >  1:1262553/1‑100 (MQ=255)
                          ccTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCAtt                                                                     >  1:209390/1‑100 (MQ=255)
                          ccTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCAtt                                                                     <  1:904548/100‑1 (MQ=255)
                            tGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGc                                                                   >  2:848386/1‑100 (MQ=255)
                            tGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGc                                                                   <  2:509339/100‑1 (MQ=255)
                                 gAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATc                                                              <  2:1262544/100‑1 (MQ=255)
                                   aaaGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTg                                                            <  1:538462/100‑1 (MQ=255)
                                      gggCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGc                                                         >  2:849614/1‑100 (MQ=255)
                                         cATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCa                                                        >  1:1233409/1‑98 (MQ=255)
                                          aTCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAgc                                                     <  1:355103/100‑1 (MQ=255)
                                            cccGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAg                                                      <  1:1407056/97‑1 (MQ=255)
                                                tataCCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATc                                               >  2:577204/1‑100 (MQ=255)
                                                    ccACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGa                                           >  2:1143018/1‑100 (MQ=255)
                                                    ccACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGAGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGa                                           >  2:854898/1‑100 (MQ=255)
                                                     cACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGAc                                          <  1:502762/100‑1 (MQ=255)
                                                      aCTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGAcc                                         <  2:1174130/100‑1 (MQ=255)
                                                         gTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTgg                                      >  2:632463/1‑100 (MQ=255)
                                                         gTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTg                                       <  1:522588/99‑1 (MQ=255)
                                                          ttCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGACATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGt                                     >  1:1102867/1‑100 (MQ=255)
                                                           tCGTAAAGAATCCAACTACGGCCTGAGGATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTg                                    <  2:937316/100‑1 (MQ=255)
                                                           tCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTg                                    <  2:962732/100‑1 (MQ=255)
                                                            cGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGa                                   >  1:519570/1‑100 (MQ=255)
                                                              tAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGAtt                                 >  2:584386/1‑100 (MQ=255)
                                                               aaaGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTa                                >  2:1356516/1‑100 (MQ=255)
                                                                aaGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTag                               >  1:505202/1‑100 (MQ=255)
                                                                  gAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagca                             >  2:768487/1‑100 (MQ=255)
                                                                  gAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagca                             >  1:795265/1‑100 (MQ=255)
                                                                   aaaaCAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagcag                            >  2:1011284/1‑100 (MQ=255)
                                                                    aaaCAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagcagc                           >  2:86904/1‑100 (MQ=255)
                                                                    aaaCAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagcagc                           <  1:768491/100‑1 (MQ=255)
                                                                     aacaacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGcc                          >  1:1253436/1‑100 (MQ=255)
                                                                     aacaacTACGACCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGTGCGAGCGCCATCCGGACCTGGTGATTAGCAGcc                          <  2:1034051/100‑1 (MQ=255)
                                                                      acaacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGCCCTGGTGATTcgcagcaa                         >  1:289477/1‑98 (MQ=255)
                                                                          cTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgc                     <  1:391245/100‑1 (MQ=255)
                                                                            aCGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcag                   <  2:1396583/100‑1 (MQ=255)
                                                                             cGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcagc                  <  2:1330857/100‑1 (MQ=255)
                                                                              ggCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcagca                 >  1:144213/1‑100 (MQ=255)
                                                                              ggCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcagca                 >  1:144238/1‑100 (MQ=255)
                                                                               gCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCaa                <  2:609086/100‑1 (MQ=255)
                                                                                ccTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAAc               <  2:727641/100‑1 (MQ=255)
                                                                                 cTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAAcc              <  2:1105508/100‑1 (MQ=255)
                                                                                       aTCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgc        <  1:920579/100‑1 (MQ=255)
                                                                                        tCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgcg       >  1:1292530/1‑100 (MQ=255)
                                                                                        tCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgcg       >  1:741451/1‑100 (MQ=255)
                                                                                          aCTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgcgtg     >  1:706199/1‑100 (MQ=255)
                                                                                           cTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgcgtgc    <  2:748409/100‑1 (MQ=255)
                                                                                            tttCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgcgtgcg   <  1:118618/100‑1 (MQ=255)
                                                                                             ttCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCTGCGTGCgg  >  1:1429617/1‑100 (MQ=255)
                                                                                                 |                                                                                               
CAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCTGCGTGCGG  >  NC_000913/428091‑428283

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