Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,600,215 +C intergenic (‑26/+23) yjjA ← / ← dnaC DUF2501 domain‑containing protein YjjA/DNA replication protein DnaC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,600,2141.C97.1% 121.9 / ‑4.0 35intergenic (‑25/+24)yjjA/dnaCDUF2501 domain‑containing protein YjjA/DNA replication protein DnaC
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (16/18);  total (16/19)
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

GTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCG  >  NC_000913/4599971‑4600439
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                  
gtgcgtgtggtcaccggtatcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGtat                                                                                                                                                                                                                             >  1:186818‑M2/238‑251 (MQ=255)
           caccggtatcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGaa                                                                                                                                                                                                                   <  1:21237‑M2/24‑1 (MQ=255)
                  atcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCCGAAACGGGGAt                                                                                                                                                                                                           >  2:253790‑M2/220‑251 (MQ=255)
                          cggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCCGAAACGGGGATTAATActc                                                                                                                                                                                                   >  2:163825‑M2/212‑251 (MQ=255)
                             ttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCttt                                                                                                                                                                                                >  2:381776‑M2/209‑251 (MQ=255)
                               ttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTa                                                                                                                                                                                               <  1:219341‑M2/44‑1 (MQ=255)
                                 gtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCt                                                                                                                                                                                            >  1:330636‑M2/205‑251 (MQ=255)
                                            actggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGc                                                                                                                                                                                  <  1:113540‑M2/57‑1 (MQ=255)
                                                          tcggcccgctgacgtgggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATcc                                                                                                                                                                    <  2:26387‑M2/71‑1 (MQ=255)
                                                                          ggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCca                                                                                                                                                   <  1:37801‑M2/88‑1 (MQ=255)
                                                                                 ggatttccgccgatcagcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTg                                                                                                                                            <  2:336376‑M2/95‑1 (MQ=255)
                                                                                         gccgatcaacggcgccctgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctggTGTGTCGCCAGTATACCGCCAGAAACGGGGAGTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCTCCCACAAACTGTTACCCAg                                                                                                                                    <  2:309536‑M2/103‑1 (MQ=255)
                                                                                          ccgatcagcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcgtcgtgcatgccggaTGTGTCGCCAGTATACCGCCAGAAGCGGGGATTAATACGCTTGACCTGTTACCCGGCTACGGGAGCGATCCCAGTTGAAGATCACCCACAAACGGTTACCCAgg                                                                                                                                   <  2:410976‑M2/104‑1 (MQ=255)
                                                                                               catcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCAt                                                                                                                              <  2:409074‑M2/109‑1 (MQ=255)
                                                                                                   ggcgcactgttgccaagcgcaaacttcttcttaatccgcatgtggctgttaaagagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagagtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAACCGGGGATTCATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACgg                                                                                                                          >  1:391013‑M2/139‑251 (MQ=255)
                                                                                                       cactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaggctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACCGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTcc                                                                                                                       >  1:24944‑M2/135‑250 (MQ=255)
                                                                                                                               tcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACt                                                                                                 <  1:383287‑M2/138‑1 (MQ=255)
                                                                                                                                   aatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTTCCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACATACTGGTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCa                                                                                           <  1:32659‑M2/144‑1 (MQ=255)
                                                                                                                                            gtggctgttaatgtgaaaagcgtaaccgaaggaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcggtcattcctgaTGTGTCGCCAGTATGCCGCCAGAAACGGGGATTAATACTCTTTTCCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACTGTCCATCACGCGTTCGCCCAGCAACGTGGTCATCTCTTCCat                                                                                 <  2:299563‑M2/154‑1 (MQ=255)
                                                                                                                                                         agaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTc                                                                    <  2:186818‑M2/167‑1 (MQ=255)
                                                                                                                                                                  taaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattccttaTGTGTCGCAAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATccc                                                            <  2:385918‑M2/175‑1 (MQ=255)
                                                                                                                                                                   agccgaatgaacgaataaagggagggaaagctccctttattattgttagcaaagtgtgcttcgttccttcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCCACGGCAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGTTCCATCACGCGTTCGCCCAGCACCTTGGTCAACTCTTCCATATTACTGTTGGTCAGCATCCCgg                                                          >  2:154060‑M2/75‑251 (MQ=255)
                                                                                                                                                                        aatgaacgaataaagggagcgaaagctccctttattattgttagaaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGAAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGATCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTg                                                       >  1:284366‑M2/70‑249 (MQ=255)
                                                                                                                                                                            aacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGtt                                                 >  1:1297‑M2/66‑251 (MQ=255)
                                                                                                                                                                                       ggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGaagaag                                        >  1:45754‑M2/55‑249 (MQ=255)
                                                                                                                                                                                         agcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAgcg                                    >  2:79466‑M2/53‑251 (MQ=255)
                                                                                                                                                                                            gaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAgcggcg                                 >  1:398476‑M2/50‑251 (MQ=255)
                                                                                                                                                                                                gctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATcg                             <  2:261605‑M2/206‑1 (MQ=255)
                                                                                                                                                                                                        tattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCt                      <  1:79466‑M2/213‑1 (MQ=255)
                                                                                                                                                                                                              tgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTATTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTtgat               >  1:70359‑M2/32‑251 (MQ=255)
                                                                                                                                                                                                               gttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTtgat               <  1:253790‑M2/220‑1 (MQ=255)
                                                                                                                                                                                                                     aaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACt        >  2:91411‑M2/25‑251 (MQ=255)
                                                                                                                                                                                                                          gtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACttttt    <  2:269081‑M2/231‑1 (MQ=255)
                                                                                                                                                                                                                          gtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACttttt    >  1:269081‑M2/20‑250 (MQ=255)
                                                                                                                                                                                                                            gcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCg  <  1:381776‑M2/233‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                  
GTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCG  >  NC_000913/4599971‑4600439

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