Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,957,548 1 bp→GG intergenic (‑117/‑744) 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,5471.G89.3% 135.8 / 12.3 56intergenic (‑116/‑745)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base . (5/1);  new base G (25/25);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,957,5480TG98.2% 141.7 / NA 56intergenic (‑117/‑744)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (29/26);  minor base C (1/0);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCT‑‑GCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGG  >  NC_000913/3957494‑3957610
                                                       ||||                                                                  
tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAgg                                               <  2:27679/76‑1 (MQ=255)
gaaTTGTCGGTGATGTGGTGGGTGTTACATATATGGCATAAGCGCGGAAAATCCGGGGTTCGCAGGTCCCTCAAgg                                               <  2:59471/75‑1 (MQ=255)
  aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAgggg                                             <  2:5318377/76‑1 (MQ=255)
  aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAgggg                                             <  1:4104952/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:4046882/76‑1 (MQ=255)
             gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                 <  1:1653847/76‑1 (MQ=255)
             gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGGCCCTCAAGGGGAGAAAAATTAg                                 >  1:189100/1‑76 (MQ=255)
              ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAgg                                <  2:4395766/76‑1 (MQ=255)
               gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                  <  2:4297470/73‑1 (MQ=255)
               gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                  >  1:4297471/1‑73 (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:5084326/1‑76 (MQ=255)
                      tGGTACATATATGGCAAAAGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                        >  1:3977248/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:147653/76‑1 (MQ=255)
                                        aGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt      >  2:4598516/1‑76 (MQ=255)
                                        aGCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt      >  2:2443414/1‑76 (MQ=255)
                                         gCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc     >  1:5625757/1‑76 (MQ=255)
                                         gCCCGGAAAATCCGG‑‑GCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc     <  1:4189929/76‑1 (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)
                                                       ||||                                                                  
TAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCT‑‑GCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGG  >  NC_000913/3957494‑3957610

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