Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,471,361 T→G 13.3% intergenic (‑32/+39) tufA ← / ← fusA translation elongation factor EF‑Tu 1/protein chain elongation factor EF‑G, GTP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,471,3610TG13.3% 131.4 / 12.7 45intergenic (‑32/+39)tufA/fusAtranslation elongation factor EF‑Tu 1/protein chain elongation factor EF‑G, GTP‑binding
Reads supporting (aligned to +/- strand):  ref base T (17/22);  new base G (3/3);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAACGCAGC  >  NC_000913/3471213‑3471508
                                                                                                                                                    |                                                                                                                                                   
ggTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCtt                                                                                                                                                     >  3:259074/1‑149 (MQ=255)
     tAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGga                                                                                                                                                <  3:277529/149‑1 (MQ=255)
      aGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgag                                                                                                                                               >  4:248694/1‑149 (MQ=255)
         cAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgagaga                                                                                                                                            <  3:334619/149‑1 (MQ=255)
          aGTACGGTGGTGATTGCAGCTGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGggagag                                                                                                                                           <  1:421435/149‑1 (MQ=21)
          aGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGagagag                                                                                                                                           <  4:259074/149‑1 (MQ=255)
          aGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGagagag                                                                                                                                           <  4:168452/149‑1 (MQ=255)
             aCGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAc                                                                                                                                        <  4:32857/149‑1 (MQ=255)
                     gATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACttt                                                                                                                                >  3:152774/1‑149 (MQ=255)
                      aTTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTg                                                                                                                               <  1:388864/149‑1 (MQ=255)
                           aGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTAtt                                                                                                                          >  3:427641/1‑149 (MQ=255)
                                tCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGGTCTCCCCTGCa                                                                                                                                                   >  4:388396/1‑119 (MQ=255)
                                tCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCa                                                                                                                                                   <  3:388396/119‑1 (MQ=255)
                                           ttACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTAtt                                                                                                          <  4:72141/149‑1 (MQ=255)
                                                   gTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACgg                                                                                                  <  1:386366/149‑1 (MQ=255)
                                                   gTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACgg                                                                                                  <  1:366176/149‑1 (MQ=255)
                                                      aaCGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCt                                                                                               >  1:147136/1‑149 (MQ=255)
                                                        cGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTc                                                                                             <  4:427641/149‑1 (MQ=255)
                                                                aTAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACgg                                                                                     >  1:292365/1‑149 (MQ=255)
                                                                     aCCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGa                                                                                <  4:215364/149‑1 (MQ=255)
                                                                                               gTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATa                                                      <  4:44110/149‑1 (MQ=255)
                                                                                                 tCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACt                                                    >  3:387053/1‑149 (MQ=255)
                                                                                                    aaTTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGAGTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCa                                                 <  4:261159/149‑1 (MQ=255)
                                                                                                      tttttCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAgg                                               >  3:423226/1‑149 (MQ=255)
                                                                                                        tttCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGaa                                             <  2:292365/149‑1 (MQ=255)
                                                                                                                 acacGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTAc                                                                                                       >  1:463991/1‑82 (MQ=255)
                                                                                                                 acacGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTAc                                                                                                       <  2:463991/82‑1 (MQ=255)
                                                                                                                    cGGCTATATTCCTTTCTATAGTGCTCTCCCTTTCAGGAGAGAGCATGGGACTTTGGTATTAACCCTTAGGCTTATTTACTACTGGCTTCAATTACGGCCTGAGCAACATTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTa                                 <  2:147136/149‑1 (MQ=255)
                                                                                                                         atatTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTatgatg                            <  4:387053/149‑1 (MQ=255)
                                                                                                                            ttCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTa                                 <  4:86278/141‑1 (MQ=255)
                                                                                                                            ttCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTa                                 >  3:86278/1‑141 (MQ=255)
                                                                                                                            ttCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCAc                         >  3:445684/1‑149 (MQ=255)
                                                                                                                                 tACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCt                    >  2:16716/1‑149 (MQ=255)
                                                                                                                                  aCTATAGTGCTCTCCCCTTCAGGAGAGAACACGGGACTTTGGTATTAACCCTTAgg                                                                                                                <  1:70895/56‑1 (MQ=255)
                                                                                                                                  aCTATAGTGCTCTCCCCTTCAGGAGAGAACACGGGACTTTGGTATTAACCCTTAgg                                                                                                                >  2:70895/1‑56 (MQ=255)
                                                                                                                                    tataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTg                 >  3:170474/1‑149 (MQ=255)
                                                                                                                                     ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTgg                <  3:248694/149‑1 (MQ=255)
                                                                                                                                       cgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTc              >  4:96049/2‑149 (MQ=255)
                                                                                                                                       cgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTc              <  3:96049/148‑1 (MQ=255)
                                                                                                                                       aGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTg                 >  3:5071/1‑146 (MQ=255)
                                                                                                                                       aGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTg                 <  4:5071/146‑1 (MQ=255)
                                                                                                                                         tGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCag            >  2:256461/1‑149 (MQ=255)
                                                                                                                                           ctctctCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCagag          >  3:189481/1‑149 (MQ=255)
                                                                                                                                              tcCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAAc       >  4:110479/1‑149 (MQ=255)
                                                                                                                                                cccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAAtt                                           <  4:24790/111‑1 (MQ=255)
                                                                                                                                                cccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAAtt                                           >  3:24790/1‑111 (MQ=255)
                                                                                                                                                   ttCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAACGCAGc  >  2:226716/1‑149 (MQ=255)
                                                                                                                                                   ttCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAACGCAGc  >  3:110258/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAACGCAGC  >  NC_000913/3471213‑3471508

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