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.C98.1% 131.5 / NA 52intergenic (‑127/‑734)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base C (22/29);  minor base A (1/0);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTGCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTC  >  NC_000913/3957495‑3957628
                                                                   |                                                                      
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:61027/76‑1 (MQ=255)
                                  gCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc                             <  2:3031416/76‑1 (MQ=255)
                                    aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAt                           <  2:270808/76‑1 (MQ=255)
                                    aaaaGCCCGGAAAATCCGGGCTTCGCAGGGCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTGAt                           <  2:279557/76‑1 (MQ=255)
                                      aaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                         <  2:1049561/76‑1 (MQ=255)
                                      aaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTa                         >  2:4742274/1‑76 (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:3479906/1‑76 (MQ=255)
                                         cccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg                      >  1:181306/1‑76 (MQ=255)
                                         cccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCg                      >  2:4770214/1‑76 (MQ=255)
                                           cGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                    <  1:3236507/76‑1 (MQ=255)
                                                        gcttCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttt        <  1:3509446/73‑1 (MQ=255)
                                                           tCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACTAAGTCATTATCCGGGGTGGtacttttatt     >  2:2411217/1‑76 (MQ=255)
                                                           tCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt     >  1:4857294/1‑76 (MQ=255)
                                                           tCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttatt     >  1:4872242/1‑76 (MQ=255)
                                                              cAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTc  >  2:3753731/1‑76 (MQ=255)
                                                              cAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTc  <  1:457233/76‑1 (MQ=255)
                                                                   |                                                                      
AATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCTGCTCGCAGGTCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTC  >  NC_000913/3957495‑3957628

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