Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,954,256 T→A 7.5% I272N (ATT→AAT)  cmoB → tRNA (cmo5U34)‑carboxymethyltransferase, carboxy‑SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,954,2560TA 7.5% 176.6 / 5.9 53I272N (ATT→AAT) cmoBtRNA (cmo5U34)‑carboxymethyltransferase, carboxy‑SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (22/27);  new base A (2/2);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-02

CTGGAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCG  >  NC_000913/1954111‑1954401
                                                                                                                                                 |                                                                                                                                                 
cTGGAAACGCTGGTTATTGATGGTGATGAGAACTCGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCTCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGc                                                                                                                                                >  1:136518/1‑149 (MQ=255)
cTGGAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGc                                                                                                                                                <  4:395101/149‑1 (MQ=255)
      aCGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGAtgt                                                                                                                                          >  3:63284/1‑149 (MQ=255)
       cGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCAATGCGGAtatc                                                                                                                                         <  3:47219/149‑4 (MQ=255)
        gCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGa                                                                                                                                        <  4:380233/149‑1 (MQ=255)
        gCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGa                                                                                                                                        <  2:147237/149‑1 (MQ=255)
          tGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGc                                                                                                                                      <  4:57147/149‑1 (MQ=255)
            gTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGt                                                                                                                                    <  2:136518/149‑1 (MQ=255)
             ttATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGtt                                                                                                                                   <  3:250748/149‑1 (MQ=255)
                 tGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTacca                                                                                                                               <  1:170149/149‑1 (MQ=255)
                 tGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTacca                                                                                                                               <  3:74722/149‑1 (MQ=255)
                        gatgaGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGa                                                                                                                                        <  1:277023/133‑1 (MQ=255)
                        gatgaGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGa                                                                                                                                        >  2:277023/1‑133 (MQ=255)
                              aaCACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGc                                                                                                                                                >  4:317004/1‑119 (MQ=255)
                              aaCACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagca                                                                                                                  <  4:40923/149‑1 (MQ=255)
                              aaCACGCTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGc                                                                                                                                                <  3:317004/119‑1 (MQ=255)
                                cacGGTGCTGGTACCAGGCGATCGTTACTCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagcagc                                                                                                                >  3:31774/1‑149 (MQ=255)
                                cacGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagcagc                                                                                                                <  1:80484/149‑1 (MQ=255)
                                          gTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGa                                                                                                      >  4:79418/1‑149 (MQ=255)
                                            aCCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGa                                                                                                                                        <  1:362540/113‑1 (MQ=255)
                                            aCCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGa                                                                                                                                        >  2:362540/1‑113 (MQ=255)
                                             ccAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatg                                                                                                   <  4:63284/149‑1 (MQ=255)
                                              cAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatgg                                                                                                  >  1:12051/1‑149 (MQ=255)
                                               aGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatgga                                                                                                 <  3:153654/149‑1 (MQ=255)
                                                 gCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatggatg                                                                                               <  3:79418/149‑1 (MQ=255)
                                                   gATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGt                                                                                             >  3:151500/1‑149 (MQ=255)
                                                    aTCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGtt                                                                                            <  4:151500/149‑1 (MQ=255)
                                                    aTCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGtt                                                                                            >  4:429235/1‑149 (MQ=255)
                                                      cGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTAc                                                                                          <  4:57701/149‑1 (MQ=255)
                                                                aaaTGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTg                                                                                <  4:20324/149‑1 (MQ=255)
                                                                   tGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGcc                                                                             <  2:12051/149‑1 (MQ=255)
                                                                            tCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGAt                                                                          <  1:371651/143‑1 (MQ=255)
                                                                            tCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGAt                                                                          >  2:371651/1‑143 (MQ=255)
                                                                            tCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTc                                                                    >  1:262398/1‑149 (MQ=255)
                                                                                       ccTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGa                                                         >  3:327542/1‑149 (MQ=255)
                                                                                            cgcgCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCgg                                                    >  2:340344/1‑149 (MQ=255)
                                                                                                  ggcgctgAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTaaaa                                              >  2:429215/1‑149 (MQ=255)
                                                                                                      ctgAAAAGCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTCCCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGt                                          <  1:457215/149‑1 (MQ=255)
                                                                                                           aaaCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAg                                     <  3:429235/149‑1 (MQ=255)
                                                                                                           aaaCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAg                                     <  3:84185/149‑1 (MQ=255)
                                                                                                            aaCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAgg                                    <  1:340344/149‑1 (MQ=255)
                                                                                                              cTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGtt                                  >  1:276374/1‑149 (MQ=255)
                                                                                                                           tgtgGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTaa                     >  3:156003/1‑149 (MQ=255)
                                                                                                                                ttttGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAAcgcg                >  2:85055/1‑149 (MQ=255)
                                                                                                                                      tcacatCCGCAATGCGGATATGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGc          >  3:98700/5‑149 (MQ=255)
                                                                                                                                      aGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGc          >  4:8467/1‑149 (MQ=255)
                                                                                                                                       cacatCCGCAATGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGaa                                      >  2:394829/4‑120 (MQ=255)
                                                                                                                                       cacatCCGCAATGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGaa                                      <  1:394829/117‑1 (MQ=255)
                                                                                                                                        atatCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGt                                          <  2:152900/115‑1 (MQ=255)
                                                                                                                                        atatCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGt                                          >  1:152900/1‑115 (MQ=255)
                                                                                                                                         tatCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGa       <  4:327542/149‑1 (MQ=255)
                                                                                                                                            ccGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTg    >  4:31087/1‑149 (MQ=255)
                                                                                                                                              gCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTgcg  >  1:282109/1‑149 (MQ=255)
                                                                                                                                                 ttGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTaaaa                                              <  4:322678/102‑1 (MQ=255)
                                                                                                                                                 ttGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTaaaa                                              >  3:322678/1‑102 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CTGGAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCG  >  NC_000913/1954111‑1954401

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