Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,957,560 +C intergenic (‑129/‑732) ppiC ← / → rep peptidyl‑prolyl cis‑trans isomerase C/rep helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,957,5581.C94.2% 175.0 / 6.5 69intergenic (‑127/‑734)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base . (4/0);  new base C (33/32);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCT‑‑GCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTC  >  NC_000913/3957494‑3957628
                                                                       |                                                                       
tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAgg                                                                 <  2:27679/76‑1 (MQ=255)
gaaTTGTCGGTGATGTGGTGGGTGTTACATATATGGCATAAGCGCGGAAAATCCGGGGTTCGCAGGTCCCTCAAgg                                                                 <  2:59471/75‑1 (MQ=255)
  aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAgggg                                                               <  1:4104952/76‑1 (MQ=255)
  aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAgggg                                                               <  2:5318377/76‑1 (MQ=255)
     gtcagtcaGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGgaga                                                           >  1:3954680/1‑76 (MQ=255)
       caggcaggGAGTGGCTGGTACATATAGGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGaaa                                                         <  2:5303664/76‑1 (MQ=255)
        agtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGaaaa                                                        <  1:4099186/76‑1 (MQ=255)
        agtcagGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGGCCCTCCAGGGGAGaaaa                                                        >  1:845718/1‑76 (MQ=255)
         gtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAAc                                                       <  2:2574174/76‑1 (MQ=255)
          tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCCCAGGTCCCCCCAGGGGAGAAAACt                                                      >  2:256689/1‑76 (MQ=255)
          tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACt                                                      >  1:2119854/1‑76 (MQ=255)
            agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                    <  2:1938691/76‑1 (MQ=255)
            agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                    <  1:3345907/76‑1 (MQ=255)
             gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                   <  1:1653847/76‑1 (MQ=255)
             gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                   <  1:4046882/76‑1 (MQ=255)
             gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGGCCCTCAAGGGGAGAAAAATTAg                                                   >  1:189100/1‑76 (MQ=255)
              ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAgg                                                  <  2:4395766/76‑1 (MQ=255)
               gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                    >  1:4297471/1‑73 (MQ=255)
               gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                    <  2:4297470/73‑1 (MQ=255)
               gAGTGGCTGGGACAGATATGGCAAAAGCCCGGAAAATCCCGGCTTCGCAGGCCCCTCCAGGGGAGAAAACTTAggg                                                 >  2:881598/1‑76 (MQ=255)
                aGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCCAGGGGGGGAAAACTAAggc                                                 >  1:850930/1‑75 (MQ=255)
                aGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGGTTCGCAGGGCCCTCCAGGGGAGAAAAATTAgggg                                                >  2:3322551/1‑75 (MQ=255)
                  tGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCct                                              >  2:5325025/1‑76 (MQ=255)
                     cTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTa                                           <  2:1672678/76‑1 (MQ=255)
                     cTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGAGCCTCTa                                           >  1:2017151/1‑76 (MQ=255)
                      tGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCACGGGCAGAAAACTTAGGGCCTCTag                                          >  2:5584087/1‑75 (MQ=255)
                      tGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                                          >  1:3977248/1‑76 (MQ=255)
                      tGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                                          >  1:5084326/1‑76 (MQ=255)
                      tGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                                          <  2:2501805/76‑1 (MQ=255)
                      tGGCACATAGAGGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCGAt                                          >  1:683832/1‑76 (MQ=255)
                          ccATATCTGGCAAACGCCCGGGAACTCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTATGGCCTCTATGGGc                                      >  1:3771546/2‑76 (MQ=255)
                          aCATATATGGCAAAAGCCCGGGAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGc                                      >  1:3739029/1‑76 (MQ=255)
                          aCATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGc                                      >  2:2977680/1‑76 (MQ=255)
                          aCATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGc                                      >  1:1973741/1‑76 (MQ=255)
                           cATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCt                                     <  1:1325010/76‑1 (MQ=255)
                            atatatGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACt                                                      >  1:5355472/1‑58 (MQ=255)
                            atatatGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACt                                                      <  2:5355471/58‑1 (MQ=255)
                            atatatGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTa                                    <  2:1624104/76‑1 (MQ=255)
                              atatGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCtcca                                  <  2:4230286/76‑4 (MQ=255)
                                atGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                                <  1:3927636/76‑1 (MQ=255)
                                 tGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAg                               >  2:4025058/1‑76 (MQ=255)
                                  ggCAAAATCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt                              <  2:4644780/76‑1 (MQ=255)
                                  ggCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACCAAGt                              >  2:5300317/1‑76 (MQ=255)
                                   gCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc                             <  1:4826092/76‑1 (MQ=255)
                                    cAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCa                            >  1:3013893/1‑76 (MQ=255)
                                     aaaaGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAt                           >  1:5347528/1‑76 (MQ=255)
                                      aaaGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAtt                          >  2:3810144/1‑76 (MQ=255)
                                       aaGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                         >  1:1262481/1‑76 (MQ=255)
                                        aGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                        >  2:4598516/1‑76 (MQ=255)
                                        aGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                        <  2:147653/76‑1 (MQ=255)
                                        aGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                        >  2:2443414/1‑76 (MQ=255)
                                         gCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                       <  1:4189929/76‑1 (MQ=255)
                                         gCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                       >  1:5625757/1‑76 (MQ=255)
                                          cccGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg                      <  2:4927619/76‑1 (MQ=255)
                                           ccGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgg                     >  2:934393/1‑76 (MQ=255)
                                            cGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                    <  2:5596078/76‑1 (MQ=255)
                                                           cttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttta       <  2:4041776/74‑1 (MQ=255)
                                                            ttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat      >  1:675557/2‑76 (MQ=255)
                                                            ttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat      >  2:3233589/2‑76 (MQ=255)
                                                            ttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat      >  1:2473354/2‑76 (MQ=255)
                                                            ttCGCAGGACCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCAACAAAGTCACTATCGCGGCTGGttatttcat      >  1:5314072/2‑76 (MQ=255)
                                                             tCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCa                            >  2:5066697/1‑53 (MQ=255)
                                                             tCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCa                            <  1:5066698/53‑1 (MQ=255)
                                                             tCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt     <  1:256689/76‑1 (MQ=255)
                                                             tCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt     <  2:890388/76‑1 (MQ=255)
                                                               gCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGt   >  1:5280272/1‑76 (MQ=255)
                                                               gCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGt   <  1:3163554/76‑1 (MQ=255)
                                                               gCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGt   <  2:396928/76‑1 (MQ=255)
                                                                cAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTc  >  2:4211582/1‑76 (MQ=255)
                                                                       |                                                                       
TAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCT‑‑GCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTC  >  NC_000913/3957494‑3957628

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