Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,147,548 T→G 10.7% intergenic (+8/‑73) rpmF → / → plsX 50S ribosomal subunit protein L32/putative phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,147,5480TG10.7% 94.9 / 2.9 45intergenic (+8/‑73)rpmF/plsX50S ribosomal subunit protein L32/putative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (18/21);  new base G (3/3);  total (21/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.40e-01

GCATGCGTCGTTCCCATGACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCA  >  NC_000913/1147407‑1147689
                                                                                                                                             |                                                                                                                                             
gCATGCGTCGTTCCCATGACGCGCCGACCTCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCCAGGTCGTCGCCACGGAATCACGCATc                                                                                                                                              >  2:722185/1‑143 (MQ=255)
       tcgtTCCCATGACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTTCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGtga                                                                                                                                       <  1:339213/143‑1 (MQ=255)
        cgtTCCCATGCCGCGCTGACCGGAGTCACCCGCCTGTCTGTAGACAAACCTTCTGGTGCAGAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGAGCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGtgat                                                                                                                                      <  2:497713/143‑1 (MQ=255)
               aTGACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCtt                                                                                                                               <  1:283402/143‑1 (MQ=255)
                tGACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTa                                                                                                                              >  1:152252/1‑143 (MQ=255)
                 gACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCGCGCATCTGCGTGATGAAGCTTAg                                                                                                                             >  1:127696/1‑143 (MQ=255)
                       cTGACCGGAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAACCACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGa                                                                                                                       <  1:24786/143‑1 (MQ=255)
                              gagtgcccAGCCTGTCGTGAGCCAAAACTTCTGGTGAAAAACACCTGCGCCACCACATCACTGCCGACGGTTACGGCCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGGGATGAAGCTTAGTGAGGATTTTccc                                                                                                                <  2:249827/137‑1 (MQ=255)
                                gTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTCGTGAGGATTTTCCCCa                                                                                                              >  1:636533/1‑143 (MQ=255)
                                 tCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACCTCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGACGCTTAGTGAGGATTTTcccccc                                                                                                             >  1:50722/1‑141 (MQ=255)
                                          tgtctgtAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTggg                                                                                                    >  1:551329/1‑143 (MQ=255)
                                               gtAGACAAAGCTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCGAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAg                                                                                               <  2:158986/143‑1 (MQ=255)
                                               gtAGACAAAACTTCTGGTGAAAAACCCCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAg                                                                                               <  1:65982/143‑1 (MQ=255)
                                                tAGACAAAACTTCTGGTGAAAAACACCTGCGGCACCACAGCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCAGCTGCGTGGTGAAGCTTAGTTAGGATTTTCCCCAGGCAACTGGGGAAAGa                                                                                              <  2:453893/143‑1 (MQ=255)
                                                 ggACAAAACTTCTGGGGAAAAACACCCGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCAGCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGAc                                                                                             <  2:350993/142‑1 (MQ=255)
                                                 aGACAAAACTTCTGGTGAAACACCCCTGCGTCACCACATCACTGCCGACGGTTACTCCCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGAc                                                                                             <  1:656910/143‑1 (MQ=255)
                                                  gACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGAcc                                                                                            <  2:254847/143‑1 (MQ=255)
                                                   aCAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCa                                                                                           >  1:287312/1‑143 (MQ=255)
                                                        aCTTCTGGTGAAAACCGCCTGCGTCACCACATCCCTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCg                                                                                      <  1:291677/143‑1 (MQ=255)
                                                           tCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGggc                                                                                   <  2:773613/143‑1 (MQ=255)
                                                           gcgggTGAAAAACACTTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCAGCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGggc                                                                                   <  2:748758/140‑1 (MQ=255)
                                                                tGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGGAAGCCCAAACCGGGCGGcga                                                                              >  1:824753/1‑143 (MQ=255)
                                                                tGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGcga                                                                              >  1:459641/1‑143 (MQ=255)
                                                                 gAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTCGTGAGGATTTTCCCCAGGCAACTGGCGAAAGACCAAACCAGgcggcgcc                                                                             >  2:438972/1‑141 (MQ=255)
                                                                    aaaCACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGcgacggt                                                                          >  2:449675/1‑141 (MQ=255)
                                                                         ccTGCGTCCCCGCACCCTCGCCGCCTGTTACAACCGCGGCCCCAAGGTCATCTCCAAGTAATCACGCATCTGCGTGATCAACCCTACTGAGGATTTTCCCCCAGAAAATTAGGAATGACTAAACCGCGCGGCGACCATACCtt                                                                     <  2:292127/143‑1 (MQ=255)
                                                                               tCACCACATCACCGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACg                                                               <  1:649515/143‑1 (MQ=255)
                                                                                         aCTTCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCGGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTg                                                     <  2:23432/143‑1 (MQ=255)
                                                                                             ccGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGt                                                 >  1:65091/1‑143 (MQ=255)
                                                                                               ggcGGTTACGACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTa                                               <  1:185551/141‑1 (MQ=255)
                                                                                               gACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAAGCGGGCGGCGACGATACCGTGACCCGGCTAACCCTGGCGtga                                               >  2:29910/1‑141 (MQ=255)
                                                                                               gACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTa                                               <  1:125472/143‑1 (MQ=255)
                                                                                                aCGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAg                                              <  2:287312/143‑1 (MQ=255)
                                                                                                  ggTTCCTGCCGCGGCCGCAAGGTCAGCGCTAGGTAAGCACGCGTCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGAt                                            <  2:28456/143‑1 (MQ=255)
                                                                                                    ttagttccGTGGCTGCGAGGTCCTCGCTAAGTAATCACGCATCTGCGTGTGGAAGCTGGGTGATGATTTTCCCTAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGt                                          <  2:622901/137‑1 (MQ=255)
                                                                                                       ctacCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGACGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGAAGATACCTTGACACGTCTAACCCTGGCGTTAGATGTcct                                       >  2:578311/1‑141 (MQ=255)
                                                                                                                ccGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGa                              >  1:447886/1‑143 (MQ=255)
                                                                                                                ccGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGa                              >  1:189172/1‑143 (MQ=255)
                                                                                                                  gCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGAtt                            >  1:740237/1‑143 (MQ=255)
                                                                                                                    aaGGTCATCGCTAAGTAATCACGCATCTGCGGGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGAtttt                          <  2:160000/143‑1 (MQ=255)
                                                                                                                               tAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGt               <  2:152252/143‑1 (MQ=255)
                                                                                                                                 cttcaTCACGCAGATGCGTGATTAAGCTTAGTGAGGATTTTCCCCAGGCAACTGTGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGa             >  1:760094/5‑143 (MQ=255)
                                                                                                                                 cttcaTCACGCAGATGCGTGATTAAGCTTAGTGAGGATTTTCCCCAGGCAACTGTGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCGTCCGTGa             >  1:439327/5‑143 (MQ=255)
                                                                                                                                   tcaTCACGCAGATGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACa           >  1:442424/3‑143 (MQ=255)
                                                                                                                                            aTCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGAAGATACCTGGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGtgcatgca  >  1:517285/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                             
GCATGCGTCGTTCCCATGACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCA  >  NC_000913/1147407‑1147689

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