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.G95.7% 114.2 / ‑0.7 46intergenic (‑116/‑745)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (18/26);  total (19/27)
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
*NC_0009133,957,5480TG97.8% 116.6 / NA 46intergenic (‑117/‑744)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (18/27);  minor base C (1/0);  total (19/27)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTGCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGG  >  NC_000913/3957495‑3957610
                                                      |||                                                               
aaTTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                             >  1:2373864/1‑76 (MQ=255)
  ttGTCAGTCAGGGAGTGGCTGGTAGATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGga                                           <  2:3321418/76‑1 (MQ=255)
   tGTCAGTCAGGGAGTGGCTGGGACATATTAGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTACCTCCAGGGgag                                          >  1:4168813/1‑76 (MQ=255)
      cagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaa                                       <  1:4356466/76‑1 (MQ=255)
         tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACt                                    <  1:3855424/76‑1 (MQ=255)
           agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAAccta                                  <  1:4707900/76‑4 (MQ=255)
              gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAggg                               <  2:4554110/76‑1 (MQ=255)
              gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAggg                               >  1:3794618/1‑76 (MQ=255)
               aGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGc                              <  2:2121655/76‑1 (MQ=255)
               aGTGGCTGGGACATATATGGCAAAAACCCGGAAAAACCCGGCTTCGCAGGTCCCTCCAGGGGAGAAAAATTAGGGc                              >  1:3217191/1‑76 (MQ=255)
                gTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                             <  1:3110394/76‑1 (MQ=255)
                gTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                             >  1:2666640/1‑76 (MQ=255)
                 tGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                 >  1:3946159/1‑71 (MQ=255)
                 tGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                 <  2:3946159/71‑1 (MQ=255)
                 tGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCct                            >  2:3268205/1‑76 (MQ=255)
                 tGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCct                            <  2:2843061/76‑1 (MQ=255)
                 tGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCct                            <  2:1905488/76‑1 (MQ=255)
                  ggCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                           <  1:744419/76‑1 (MQ=255)
                   gCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctct                          <  2:3803222/76‑1 (MQ=255)
                     tGGTACATATATGGCAAGAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                        <  2:3761809/76‑1 (MQ=255)
                      gctACATATATGGCAAAAGCACGGAAAATCCGGGCGTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATg                       >  2:2490405/3‑76 (MQ=255)
                       gTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATgg                      <  2:704951/76‑1 (MQ=255)
                       gTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATgg                      >  2:1560184/1‑76 (MQ=255)
                       gTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATgg                      >  1:4823847/1‑76 (MQ=255)
                        tACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATggg                     >  2:4790288/1‑76 (MQ=255)
                        tACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATggg                     >  1:4058316/1‑76 (MQ=255)
                           atatctGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTGAGGGCCTCTATGGGCTa                  >  2:3075676/1‑76 (MQ=255)
                           atatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTa                  <  1:1969954/76‑1 (MQ=255)
                            tatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTACGGCCTCTATGGGCTAc                 >  2:2405603/1‑76 (MQ=255)
                             atatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                <  1:3861387/76‑1 (MQ=255)
                               atGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa              <  1:2297354/76‑1 (MQ=255)
                                tGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaacg             >  2:4538089/1‑74 (MQ=255)
                                 ggCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt            <  1:17239/76‑1 (MQ=255)
                                  gCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc           <  2:3031416/76‑1 (MQ=255)
                                  gCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc           <  2:61027/76‑1 (MQ=255)
                                    aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAt         <  2:270808/76‑1 (MQ=255)
                                    aaaaGCCCGGAAAATCCGGGCTTCGCAGGGCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTGAt         <  2:279557/76‑1 (MQ=255)
                                      aaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa       >  2:4742274/1‑76 (MQ=255)
                                      aaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa       <  2:1049561/76‑1 (MQ=255)
                                       aGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt      <  2:666846/76‑1 (MQ=255)
                                        gCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc     <  2:3719838/76‑1 (MQ=255)
                                        gCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc     <  2:2373863/76‑1 (MQ=255)
                                         cccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg    >  1:181306/1‑76 (MQ=255)
                                         cccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg    >  2:4770214/1‑76 (MQ=255)
                                         cccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg    >  1:3479906/1‑76 (MQ=255)
                                           cGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg  <  1:3236507/76‑1 (MQ=255)
                                                      |||                                                               
AATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTGCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGG  >  NC_000913/3957495‑3957610

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