Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,327,860 A→G E242E (GAA→GAG crfC → clamp‑binding sister replication fork colocalization protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,327,8600AG97.8% 148.8 / ‑4.5 45E242E (GAA→GAGcrfCclamp‑binding sister replication fork colocalization protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (21/23);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

TGGACGTCGATTTTCCTTTTTCTGCTTACGCCGCCATTGAGCATATTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTT  >  NC_000913/4327724‑4327984
                                                                                                                                        |                                                                                                                            
tGGACGTCGATTTTCCTTTTTCTGCTTACGCCGCCATTGAGCATATTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGcc                                                                                                                           >  2:226945/1‑140 (MQ=255)
      tCGATTTTCCTTTTTCTGCTTACGCCGCCATTGAGCATATTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATCCCCCCGGGCCAAATGAGGc                                                                                                                            <  2:625471/133‑1 (MQ=255)
      tCGATTTTCCTTTTTCTGCTTACGCCGCCATTGAGCATATTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGc                                                                                                                            >  1:625471/1‑133 (MQ=255)
      tCGATTTTCCTTTTTCTGCTTACGCCGCCATTGAGCATATTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCaa                                                                                                                     >  2:219660/1‑140 (MQ=255)
      tCGATTTTCCTTTTTCTGCTTACGCCGCCATTGAGCATATTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCaa                                                                                                                     >  1:98388/1‑140 (MQ=255)
             tCCTTTTTCTGCTTACGCCGCCATTGAGCATATTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAAc                                                                                                                    <  1:1016740/134‑1 (MQ=255)
             tCCTTTTTCTGCTTACGCCGCCATTGAGCATATTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAAc                                                                                                                    >  2:1016740/1‑134 (MQ=255)
               cTTTTTCTGCTTACGCCGCCATTGAGCATATTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTg                                                                                                            <  2:98388/140‑1 (MQ=255)
                    tCTGCTTACGCCGCCATTGAGCATATTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCaaaa                                                                                                       <  1:1031426/140‑1 (MQ=255)
                    tCTGCTTACGCCGCCATTGAGCATATTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCaaaa                                                                                                       >  1:342805/1‑140 (MQ=255)
                               cgccATTGAGCATATTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAAcc                                                                                            <  1:219660/140‑1 (MQ=255)
                                  cATTGAGCATATTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACcagc                                                                                         >  1:386980/1‑140 (MQ=255)
                                   aTTGAGCATATTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACcagca                                                                                        >  2:317387/1‑140 (MQ=255)
                                       aGCATATTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTg                                                                                    >  2:283969/1‑140 (MQ=255)
                                       aGCATATTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTg                                                                                    >  1:1112187/1‑140 (MQ=255)
                                          atatTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCa                                                                                 <  2:912968/140‑1 (MQ=255)
                                          atatTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATACCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCa                                                                                 <  2:1058320/140‑1 (MQ=255)
                                               cccGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACgcgc                                                                            <  2:1182709/140‑1 (MQ=255)
                                                   tGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTcg                                                                        >  1:698173/1‑140 (MQ=255)
                                                   tGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTcg                                                                        >  1:117889/1‑140 (MQ=255)
                                                   tGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTcg                                                                        >  1:656469/1‑140 (MQ=255)
                                                    gATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTcgg                                                                       <  1:226945/140‑1 (MQ=255)
                                                    gATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTcgg                                                                       <  2:689820/140‑1 (MQ=255)
                                                     aTTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCagcagc                                                                                      <  2:375073/124‑1 (MQ=255)
                                                     aTTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCagcagc                                                                                      >  1:375073/1‑124 (MQ=255)
                                                       tGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTcggcgg                                                                    >  2:882547/1‑140 (MQ=255)
                                                           gTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACt                                                                <  2:117889/140‑1 (MQ=255)
                                                             ggAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTgg                                                              <  1:317387/140‑1 (MQ=255)
                                                              gAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGc                                                             >  1:354932/1‑140 (MQ=255)
                                                                 tttGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGt                                                          >  1:1177903/1‑140 (MQ=255)
                                                                  ttGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTg                                                         <  2:698173/140‑1 (MQ=255)
                                                                    gTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCt                                                       >  1:27898/1‑140 (MQ=255)
                                                                    gTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCt                                                       <  2:656469/140‑1 (MQ=255)
                                                                         tCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGAGGTGGAtt                                                  >  2:321621/1‑140 (MQ=255)
                                                                            aGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATtata                                               <  2:1177903/140‑1 (MQ=255)
                                                                                 ggCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCaa                                          <  2:354932/140‑1 (MQ=255)
                                                                                                 cGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCg                          <  1:126365/140‑1 (MQ=255)
                                                                                                   gTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGGTTATAAGCAAATGAAAAAGAACTcccct                        >  1:862930/1‑137 (MQ=255)
                                                                                                   gTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGGGCTGGGGTTTAAGCCACTGGAAAAGATCTCCgct                        >  2:341508/1‑138 (MQ=255)
                                                                                                   gTCAGTCGACGTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGAt                        <  1:829524/140‑1 (MQ=255)
                                                                                                      aGTTGACGTTACTGGATCCCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTAGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGaa                     <  2:108797/140‑1 (MQ=255)
                                                                                                             gTTACTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTcc              <  1:38311/140‑1 (MQ=255)
                                                                                                                 cTGGATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGa          <  2:27898/140‑1 (MQ=255)
                                                                                                                   ggATACCCCCGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAgg        <  1:42051/140‑1 (MQ=255)
                                                                                                                         ccccGGGCCAAATGAGGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGGGCTGGGTTATAAGCAACTGAAAAAGATCTCCCATGAAAAGGGCCGGGAggggggtt  >  1:949714/1‑134 (MQ=255)
                                                                                                                                        |                                                                                                                            
TGGACGTCGATTTTCCTTTTTCTGCTTACGCCGCCATTGAGCATATTCCCGTGATTGAAGTGGAGTTTGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTT  >  NC_000913/4327724‑4327984

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: