Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,134,229 T→C 8.2% Y33Y (TAT→TAC)  yghT → putative ATP‑binding protein YghT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,134,2290TC 8.2% 104.1 / 2.9 50Y33Y (TAT→TAC) yghTputative ATP‑binding protein YghT
Reads supporting (aligned to +/- strand):  ref base T (23/21);  new base C (2/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01

CAAATAAATATTTTTAATGGATATTTAAATTAAAGGATATATTCATGCAGTCAATAACACCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACATAAAACAATTGAACGAAATCAGGTGAAAACAGCAAAAAAATTGCCTGGACC  >  NC_000913/3134087‑3134369
                                                                                                                                              |                                                                                                                                            
cAAATAAATATTTTTAATGTATATTGAAATTACCGACTACAACCATGCAGTCAATAACACCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAatat                                                                                                                                              <  2:332915/143‑1 (MQ=255)
 aaataaatATTTTTAATGGATATTTAAATTAAAGGATATATTCATGCAGTCAATAACACCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGCGCATCTTACTACCGGTGTCGATAAATATg                                                                                                                                             >  1:581251/1‑143 (MQ=255)
 aaataaatATGTTTAATGGATATTTAAATTAAAGGATATATTCATGCAGTCAATAACACCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATg                                                                                                                                             <  1:452171/143‑1 (MQ=255)
    taaataTTTTTAATGGATATTTAAATTAAAGGATATATTCATGCAGTTAATAACACCTCCATTCATTGCCGTTCTTGGTAGCGATGATTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGtg                                                                                                                                          <  1:270072/143‑1 (MQ=255)
      aataTTTTTAATGGATATTTAAATTAAAGGATATATTCATGCAGTCAATAACACCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGtgct                                                                                                                                        <  1:97878/143‑1 (MQ=255)
         atTTTTAATGGATATTTAAATTAAAGGATATATTCATGCAGTCAATAACACCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCACGCAAGTCAACGGTGTGTGAACATCTTATTACCGGTGTCGAAAAAAATGGTGCTGcc                                                                                                                                     >  1:212117/1‑143 (MQ=255)
             ttAATGGATATTTAAATTAAAGGATATATTCATGCAGTCAATAACACCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGaaa                                                                                                                                 <  1:532393/143‑1 (MQ=255)
                      atTTAAATTAAAGGATATATTCATGCAGTCAATAACACCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGaattaatt                                                                                                                        >  1:440915/1‑143 (MQ=255)
                      atTTAAATTAAAGGATATATTCATGCAGTCAATAACACCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTAGAAAACTATGGTGCTGCCGAAAGAGATCAtt                                                                                                                        >  1:55872/1‑143 (MQ=255)
                         taaattaaaGGATATATTCATGCAGTCAATAACACCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCCCTGTCGAACAACATGGCCCTGCCGAAACAGTCCATTTac                                                                                                                     >  1:326330/1‑142 (MQ=255)
                          aaagtgaaGGTTATATTCATGCAGTCAATAACACCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAgg                                                                                                                    <  2:171000/137‑1 (MQ=255)
                              taaaGGATATATTCATGCAGTCAATAACACCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGaaaa                                                                                                                <  1:186130/143‑1 (MQ=255)
                                        atTCATGCAGTCAATAACACCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATACGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGaa                                                                                                      >  1:748870/1‑143 (MQ=255)
                                        atTCATGCAGTCAATAACACCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATACGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGaa                                                                                                      >  1:338378/1‑143 (MQ=255)
                                               cAGTCAATAACACCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTGGGAGAACCGGCGGGATATgtcgg                                                                                               >  1:64896/1‑143 (MQ=255)
                                                  tCAATAACACCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGATAAGTATGGTGCTGCCGTAGGAGTTTATTTAGGTAAACAGACCGAAAATgtaggtcg                                                                                            >  1:249198/1‑143 (MQ=255)
                                                  tCAATAACACCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATgtcggtcg                                                                                            >  1:587141/1‑143 (MQ=255)
                                                    aaTAACACCTGCATTAATTGCAGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTg                                                                                          <  1:191068/143‑1 (MQ=255)
                                                      tagcacCTCCATTAATTGCCGTTATTGTTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCa                                                                                        <  2:719559/140‑1 (MQ=255)
                                                      gaaCCCCTCGATTAATTGCGGTTATGGGTAGCGATGGTGCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCa                                                                                        <  2:719497/142‑1 (MQ=255)
                                                       aaCATCTCTCTTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAg                                                                                       <  1:605405/143‑1 (MQ=255)
                                                        acacCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGCGATGTAGTATCCGATTTAAAGAGCTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGt                                                                                      >  1:349040/1‑143 (MQ=255)
                                                             tCCATTAATTGCCGTTATTGGTAGCGATGGTTAAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACaa                                                                                 >  2:33297/1‑143 (MQ=255)
                                                                           ttATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGACAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATACCCGTTCATg                                                                   >  2:622043/1‑143 (MQ=255)
                                                                                     cGATGGTTCAGCCAAGTCAACGGTGTGTGAACATATTGATACCGTTGTCGAAAAATCTGGTGCTGCCGAGATTCTTCATTTAGGAAAACAGGCCGAAAATGTCAGCAGTGCAGTGAAAAAATTAACGTTGATGGGAAAATGCt                                                         >  2:53813/1‑143 (MQ=255)
                                                                                      gATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAACATTACCGTTGATGGGAAAATCCtt                                                        >  1:593521/1‑143 (MQ=255)
                                                                                              aGGCCAGTCAACGGTGTGTGATCGTCTTGTTACCGTTGTCGAAAAATTGGTTTCTGCCTAAAGAGTTCATTTCGGAAGACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGTTGGGAAAATCCTTACATaaaa                                                <  2:48153/143‑1 (MQ=255)
                                                                                                       aaCGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAACTTACCGTTGATGGGAAAATCCTTACATAAAACAAATGAAc                                       >  2:576235/1‑143 (MQ=255)
                                                                                                        aCGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGTGTTCACCGCGGATTGTTGGCCGAAAACCCCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACATAAAACAATTGAACg                                      >  1:174812/1‑143 (MQ=255)
                                                                                                         cGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGCAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGCAAAATCCTTCCCTACAACAATTGAAcca                                     >  2:720368/1‑141 (MQ=255)
                                                                                                              tttgAACTTTTTATTACCGTTGTCGAAAAATATGGTGCGGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACATAAAACAATTGAACGAAATCa                                <  1:45169/141‑1 (MQ=255)
                                                                                                              tgttAAAATTTTATTATAGTTGTCGAAAATTATGGTGCTGCTGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGGCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACATAAAACAATTGAACGAAATCa                                <  1:772849/143‑1 (MQ=255)
                                                                                                                 gAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAACATCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACCTAAAACAATTGAACAAACTCAAGt                             >  2:137016/1‑143 (MQ=255)
                                                                                                                  aaCATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACATAACACAATTGAACGAAATCAGGTg                            >  2:560317/1‑143 (MQ=255)
                                                                                                                  aaCATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACATAAAACAATTGAACGAAATCAGGTg                            <  2:587141/143‑1 (MQ=255)
                                                                                                                    cATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGGAAATCCTTACATAAAACAATTGAAGGAAATCAGGTGaa                          >  1:357149/1‑143 (MQ=255)
                                                                                                                     aacTTATTACCGTTTTCGAAAAATATGGTGCTGCAGAAAGAGTTCATTTAGGAAAACAAGCCGGAAATAACGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCATTACATAAAACAATTGAACGAAAACAGGTGaaa                         >  2:272707/3‑143 (MQ=255)
                                                                                                                     aTCTTATTACCGTTGTCGAAGAATATGGTGCTGCCGAAAGAGTTCATTTAGCAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACATAAAACAATTGAACGAAATCAGGTGaaa                         >  2:396737/1‑143 (MQ=255)
                                                                                                                     aTCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACATAAAACAATTGAACGAAATCAGGTGaaa                         >  2:624169/1‑143 (MQ=255)
                                                                                                                      tCTTATTACCGTAGTGGAAAAACATGGGGCCGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACATAAAACAATTGAACGAAATCAGGTGaaaa                        <  1:123705/143‑1 (MQ=255)
                                                                                                                        ttattaCCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACATAAAACAATTGAACGAAATCAGGTGAAAACa                      <  1:465798/143‑1 (MQ=255)
                                                                                                                         tattaCCGTTGTCGAAAAATACGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACATAAAACAATTGAACGAAATCAGGTGAAAACAg                     <  2:338378/143‑1 (MQ=255)
                                                                                                                         tattaCCGTTGTCGAAAAATACGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACATAAAACAATTGAACGAAATCAGGTGAAAACAg                     <  2:748870/143‑1 (MQ=255)
                                                                                                                         tattaCCGTTGTCGAAAAAGATGGTGCTGCCGAAAGAGTTCATTTAGGAAAAAAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCAGACATAAAACAATTGAACGAAATCAGGTGAAAACAg                     >  2:39319/1‑143 (MQ=255)
                                                                                                                               cGTTGTCGAAACATGTCGTGCTGTCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACATAAAACAATTGAACGAAATCAGGTGAAAACAGCaaaaa               <  2:702833/143‑1 (MQ=255)
                                                                                                                                 ttGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACATAAAACAATTGAACGAAATCAGGTGAAAACAGCaaaaaaa             >  1:673406/1‑143 (MQ=255)
                                                                                                                                   gacGAAAAATATGGGGCAGCAGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACATAAAACAATTGAACGAAATCAGGTGAAAACAGCAAAAAAAtt           <  1:116857/141‑1 (MQ=255)
                                                                                                                                    tCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACATAAAACAATTGAACGAAATCAGGTGAAAACAGCAAAAAAATTg          <  2:144318/143‑1 (MQ=255)
                                                                                                                                    tCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACATAAAACAATTGAACGAAATCAGGTGAAAACAGCAAAAAAATTg          <  1:653855/143‑1 (MQ=255)
                                                                                                                                     ctaagaaTATGGTGCTGCCGAAAGAGTTCATTTAGTAAATGAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCTTTTATGGGAAAATCCTTACATAAAACAATTGAACGAAATCAGTTGAAAACAGCAAAAAAATTGc         <  2:234000/138‑1 (MQ=255)
                                                                                                                                      gAAAAATATGGTGCTGTCGAAAAAGTTCATTTAGGAAAACAGTCCGGAAATGTCGGTCGTGCAGTGACAAAAGTACCGTTGATGGGAAAATCCTTACATAAAACAATTGAACGAAATCAGGTGAAAACAGCAAAAAAATTGcc        <  1:456365/143‑1 (MQ=255)
                                                                                                                                            tatGGTGCTGCCGAAAGAGTTCATTTAAGAAAACAGTGCTATGATGTCCGTCGTGCAGTGACAAAATTACGGACGACGGGATAATCCTCACGTAGGACAAGTGAACGAAATCAGGGGAAAACAGCAAGAGCATTGCCTGTAcc  >  2:506117/1‑143 (MQ=255)
                                                                                                                                            tatGGTGCTGCCGAAAGAGTACATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAAGCCTTACATAAAACAATTGAACGAAATCAGGAGAAAACAGCAAAAAAAATGCCTGCAcc  >  2:397139/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                            
CAAATAAATATTTTTAATGGATATTTAAATTAAAGGATATATTCATGCAGTCAATAACACCTCCATTAATTGCCGTTATTGGTAGCGATGGTTCAGGCAAGTCAACGGTGTGTGAACATCTTATTACCGTTGTCGAAAAATATGGTGCTGCCGAAAGAGTTCATTTAGGAAAACAGGCCGGAAATGTCGGTCGTGCAGTGACAAAATTACCGTTGATGGGAAAATCCTTACATAAAACAATTGAACGAAATCAGGTGAAAACAGCAAAAAAATTGCCTGGACC  >  NC_000913/3134087‑3134369

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