Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,663,482 A→C 8.2% K14T (AAG→ACG)  suhB → inositol‑phosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,663,4820AC 8.2% 104.0 / 2.6 51K14T (AAG→ACG) suhBinositol‑phosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base A (24/21);  new base C (3/2);  total (27/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 = 5.27e-01

TATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTACCCAC  >  NC_000913/2663345‑2663622
                                                                                                                                         |                                                                                                                                            
tattCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCgg                                                                                                                                         <  1:266124/143‑1 (MQ=255)
gatgcACGCGTCTTATCATAAAGCGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTGCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCgg                                                                                                                                         <  2:371049/139‑1 (MQ=255)
 attCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCCGTGAGAGAGACCGCTGCATCCGATGCTGAACATCGCCGTGCGCGCCGCGCGCAAGGCggg                                                                                                                                        >  1:626465/1‑143 (MQ=255)
 attCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCgat                                                                                                                                        <  1:689231/143‑3 (MQ=255)
   tCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCACCGTATTCCCGTTCTCTAACATCCAGTGAGAGAGACGGATGCATCCGTTGCTGAAGATCGGCGTGCGCGCAGCGCACACGGCGCGta                                                                                                                                      >  2:250944/1‑143 (MQ=255)
      cGCATCTTACCATAAACCGAAGACAGAGGCCGATCTCGCTGCTATACTCTGAGCCGCTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCgatgttca                                                                                                                                   <  1:502094/143‑8 (MQ=255)
          tCTTATCGTAACACGATGATAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaat                                                                                                                               <  2:418296/143‑1 (MQ=255)
                 aTACAACGACGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTCTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACACAGACGTGCGCGCAGCGCGCTAGGCGGGTAATTTAATTGCCaaa                                                                                                                        <  1:208286/143‑1 (MQ=255)
                  tAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGCGACCGATGCATGCGATGCTGAACATCGCAGTGCGCACAGCGCGCACGGCGGGTGATGCAATTGCCaaaa                                                                                                                       >  2:770164/1‑143 (MQ=255)
                  tAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCATTCACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAACTTAATTGCCaaaa                                                                                                                       >  1:797123/1‑143 (MQ=255)
                     aaCGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGACGGGAAATTTAATTGCCAAAAACt                                                                                                                    >  2:661394/1‑143 (MQ=255)
                              aGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCACAAACTATTCAAccc                                                                                                           >  1:781581/1‑143 (MQ=255)
                                aTGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAAGAGGGAATTTTAGTTGCCAAAAACTATGAAACCCCg                                                                                                         >  1:598137/1‑143 (MQ=255)
                                  gCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGa                                                                                                       >  1:724937/1‑143 (MQ=255)
                                       tctgGCTGCTATGCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGGGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTg                                                                                                  <  2:545116/143‑1 (MQ=255)
                                         tcGCTGCTATACTCTGCGCCGTTTTCCCGTTCATTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATAGCCAAAAACTATGAAACCCCGGACGCTGTa                                                                                                >  2:462489/1‑143 (MQ=255)
                                           gctgctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCAGCTGATGAAGAATATCGCCGTGCGCGCAGCGCGCAAGGCCGGTAATTTAATTGCCAAAACCTATGAAACCCCGGACGCTGgaga                                                                                              >  2:628813/1‑143 (MQ=255)
                                            ctgctATACTCCGCGCCGTTTTCCCGGTCTTTAACAGCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGaa                                                                                             <  2:826061/143‑1 (MQ=255)
                                                tataCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCACCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGa                                                                                         >  1:268870/1‑143 (MQ=255)
                                                            tggcgtCCCGTTCGTCAGCATCCCGTGAGAGAGTCCGAGGCATCCGATGCTGAACATCGCCGTGCGCGCTGCGCGCAAGGCGGTTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCa                                                                             <  2:299070/138‑1 (MQ=255)
                                                            ctttttCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGTATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCa                                                                             <  2:116993/141‑1 (MQ=255)
                                                             gTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAg                                                                            >  1:306090/1‑143 (MQ=255)
                                                                   ccGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGa                                                                      >  1:567212/1‑143 (MQ=255)
                                                                     gTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCAGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGAtt                                                                    >  1:746450/1‑143 (MQ=255)
                                                                       tCTTTAACATCCAGTGAGAGAGACGGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTc                                                                  <  2:128940/143‑1 (MQ=255)
                                                                         tttAACTTCCAGTGGGAGGGACCGAGGCATCCGATGCTGAACATCGCCGTGCGCGCACCGCGCAAGGCGGGTAATTTCATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGt                                                                <  2:195027/143‑1 (MQ=255)
                                                                         tttAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGAGCGCAGCGCTCAAGGCGAGTACCTTAGTGGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGt                                                                >  1:885262/1‑143 (MQ=255)
                                                                           tAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGGAGCGAGGCAGAAAGGCAGTAACGATTTccggg                                                              >  2:245126/1‑139 (MQ=255)
                                                                               aTCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCCGCGCGCAAGGCGGGTAATTTAATTGCCACAAACTATGAAACCCCGCACGCCGCAGAAGCGAGCCAGCAAGGCAGTAACGATTTCGTGACCaa                                                          >  2:795517/1‑143 (MQ=255)
                                                                                   aGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAACACTATGACACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGAATACGTGCCCAACGTa                                                      >  2:332002/1‑143 (MQ=255)
                                                                                         agagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAATCTAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATaaa                                                <  1:167692/143‑1 (MQ=255)
                                                                                           agaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGc                                              >  1:872694/1‑143 (MQ=255)
                                                                                           agaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGc                                              >  1:415974/1‑143 (MQ=255)
                                                                                               cGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGcc                                          <  1:159345/143‑1 (MQ=255)
                                                                                                gATGCATCCGATGCTGAACATCGCCGTGAGCGCAGGGCGCGAGGCGGGGAATTTAATTGCCAAAAACAACGAAACCCTTGACGCTGTAGAAGCGAGCCAGAAAGGCAGCAACGGTTTCGTGACCAACGTAGATGACGCTGCCg                                         >  2:265430/1‑143 (MQ=255)
                                                                                                      tCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCgg                                   >  1:100113/1‑143 (MQ=255)
                                                                                                      accGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCgg                                   <  1:504622/142‑1 (MQ=255)
                                                                                                       ccGATGCTGAGCATCGCCGTGCGCGCAGCGCGCCGGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGt                                  <  2:744537/143‑1 (MQ=255)
                                                                                                            gCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGCAACCCCGGACGCTGTAGAAGCGAGCCAGAAGGGCAGTAACGATTTCGTGACCAACGTAGATCAAGCTGGCGAAGCGGTGATTa                             >  1:384842/1‑143 (MQ=255)
                                                                                                               gAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCg                          >  1:361347/1‑143 (MQ=255)
                                                                                                               gAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCg                          >  1:780460/1‑143 (MQ=255)
                                                                                                                     cGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGa                    >  1:538614/1‑143 (MQ=255)
                                                                                                                       ccGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGTCACGGtt                  >  1:398829/1‑143 (MQ=255)
                                                                                                                              cgcagcgcggaAGGCGGGTAATTTAACTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATCTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAAt           <  2:81452/143‑1 (MQ=255)
                                                                                                                              cgcagcgcgcaATGCGGGTAATTTAGTGGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAAt           <  2:547861/143‑1 (MQ=255)
                                                                                                                              accagcccgcaAGGCGAGAAATTAAAGAGCCAAAAACTACAACGCGCCGGGCTTGATAAAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAAt           <  1:644457/141‑1 (MQ=255)
                                                                                                                                  gcgggcaAGGCGGGTGATTGGATTGCCAAAAACTATGAAGCGCCGGTCGTCGGAGAAGCGAGCCAGAAAGGCAGGAACGATTTCGTGAGCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTa       <  2:767581/143‑1 (MQ=255)
                                                                                                                                  gcgggcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATCATCGACACGATTCGTAAATCTTa       <  1:382003/143‑1 (MQ=255)
                                                                                                                                  gcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTAAGAAACCCCGGACGCTGAAGAAGACATCTATTATGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTa       <  1:201706/143‑1 (MQ=255)
                                                                                                                                       caATGCGCGTAATTTAATTGCCAAAAACTATGAAACCCGTGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACTTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTACCcac  <  1:389497/143‑1 (MQ=255)
                                                                                                                                       caAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTACCcac  <  1:813633/143‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                            
TATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTACCCAC  >  NC_000913/2663345‑2663622

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