Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,601,825 G→C 7.7% P88A (CCC→GCC) ‡ ftsX ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,601,8250GC 7.7% 181.7 / 5.5 52P88A (CCC→GCC) ‡ftsXputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (25/23);  new base C (2/2);  total (27/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.87e-01

CGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGCGAGGAT  >  NC_000913/3601708‑3601970
                                                                                                                     |                                                                                                                                                 
cGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCagag                                                                                                                                        >  3:657877/1‑129 (MQ=255)
cGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCagag                                                                                                                                        <  4:657877/129‑1 (MQ=255)
                     ccAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGaa                                                                                                                                      >  2:681488/1‑110 (MQ=255)
                     ccAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGaa                                                                                                                                      <  1:681488/110‑1 (MQ=255)
                         cGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCag                                                                                                                                          >  4:541591/1‑102 (MQ=255)
                         cGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCag                                                                                                                                          <  3:541591/102‑1 (MQ=255)
                                                tttGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGcc                                                                     >  4:826369/1‑148 (MQ=255)
                                                        gACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCag                                                                                                                                          <  2:784354/71‑1 (MQ=255)
                                                        gACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCag                                                                                                                                          >  1:784354/1‑71 (MQ=255)
                                                             ataataCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGa                                                                >  1:44067/1‑140 (MQ=255)
                                                             ataataCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGa                                                                <  2:44067/140‑1 (MQ=255)
                                                               aataCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAgcg                                                      <  3:826369/148‑1 (MQ=255)
                                                                 taCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATgcgc                                                                      <  2:279241/130‑1 (MQ=255)
                                                                 taCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATgcgc                                                                      >  1:279241/1‑130 (MQ=255)
                                                                     gCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTgg                                                                                                       <  2:506697/93‑1 (MQ=255)
                                                                     gCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTgg                                                                                                       >  1:506697/1‑93 (MQ=255)
                                                                                   ttAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACgg                                                                                                >  3:1132675/1‑86 (MQ=255)
                                                                                   ttAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACgg                                                                                                <  4:1132675/86‑1 (MQ=255)
                                                                                    tAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGaa                                                               >  3:192917/1‑118 (MQ=255)
                                                                                    tAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGaa                                                               <  4:192917/118‑1 (MQ=255)
                                                                                      aCGTTTTTGTACACCATATATCAGACGCTGGGCAGCGTCAGAGAAATGGCGATa                                                                                                                             <  1:868275/54‑1 (MQ=255)
                                                                                      aCGTTTTTGTACACCATATATCAGACGCTGGGCAGCGTCAGAGAAATGGCGATa                                                                                                                             >  2:868275/1‑54 (MQ=255)
                                                                                      aCGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCtttt                                                                                       <  2:422869/92‑1 (MQ=255)
                                                                                      aCGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCtttt                                                                                       >  1:422869/1‑92 (MQ=255)
                                                                                         tttttGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATa                                                         >  3:513167/1‑119 (MQ=255)
                                                                                         tttttGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATa                                                         <  4:513167/119‑1 (MQ=255)
                                                                                                    cATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGtt                                                                                              <  2:942813/71‑1 (MQ=255)
                                                                                                    cATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGtt                                                                                              >  1:942813/1‑71 (MQ=255)
                                                                                                      tataACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTgg                                                                                                       >  3:110904/1‑60 (MQ=255)
                                                                                                      tataACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTgg                                                                                                       <  4:110904/60‑1 (MQ=255)
                                                                                                      tataACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTc                                                                                      >  3:123892/1‑77 (MQ=255)
                                                                                                      tataACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTc                                                                                      <  4:123892/77‑1 (MQ=255)
                                                                                                       ataACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCGCCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGt                                  >  2:285616/1‑128 (MQ=255)
                                                                                                       ataACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCAt                                                          <  1:413629/104‑1 (MQ=255)
                                                                                                       ataACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCAt                                                          >  2:413629/1‑104 (MQ=255)
                                                                                                       ataACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGt                                  <  1:285616/128‑1 (MQ=255)
                                                                                                        taACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTg                          >  2:923959/1‑135 (MQ=255)
                                                                                                        taACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTg                          <  1:923959/135‑1 (MQ=255)
                                                                                                         aaCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCAt                                                          <  4:673193/102‑1 (MQ=255)
                                                                                                         aaCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCAt                                                          >  3:673193/1‑102 (MQ=255)
                                                                                                         aaCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGc                         <  1:22032/135‑1 (MQ=255)
                                                                                                         aaCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGc                         >  2:22032/1‑135 (MQ=255)
                                                                                                           ctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATa                                                         >  3:1079239/3‑101 (MQ=255)
                                                                                                           ctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATa                                                         <  4:1079239/99‑1 (MQ=255)
                                                                                                           cAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTg                                                                                            <  4:189299/66‑1 (MQ=255)
                                                                                                           cAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTg                                                                                            >  3:189299/1‑66 (MQ=255)
                                                                                                           cAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTg          >  3:460856/1‑148 (MQ=255)
                                                                                                            aGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTgg         >  2:427088/1‑148 (MQ=255)
                                                                                                             gACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGc        >  3:585455/1‑148 (MQ=255)
                                                                                                             gACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCtt                                                                                         >  2:1013227/1‑67 (MQ=255)
                                                                                                             gACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCtt                                                                                         <  1:1013227/67‑1 (MQ=255)
                                                                                                              aCGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGCg       <  1:427088/148‑1 (MQ=255)
                                                                                                                   gggCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGCGAGGAt  <  4:460856/148‑1 (MQ=255)
                                                                                                                     gCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCaa                                                                           <  3:586388/73‑1 (MQ=255)
                                                                                                                     gCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCaa                                                                           >  4:586388/1‑73 (MQ=255)
                                                                                                                     |                                                                                                                                                 
CGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGCGAGGAT  >  NC_000913/3601708‑3601970

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: ‑