Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,719,768 G→A *97* (TAG→TAA yiaG → HTH_CROC1 family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,719,7680GA97.7% 135.4 / ‑6.7 44*97* (TAG→TAAyiaGHTH_CROC1 family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/22);  minor base T (0/1);  total (21/23)
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

GTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTA  >  NC_000913/3719640‑3719902
                                                                                                                                |                                                                                                                                      
gTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATcc                                                                                                                              >  1:136534/1‑139 (MQ=255)
 ttttGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCa                                                                                                                             <  1:676872/139‑1 (MQ=255)
        cGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTAtt                                                                                                                      <  2:136534/139‑1 (MQ=255)
         gTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTg                                                                                                                     >  2:6690/1‑139 (MQ=255)
          tATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGc                                                                                                                    >  2:243080/1‑139 (MQ=255)
          tATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGc                                                                                                                    >  1:706703/1‑139 (MQ=255)
           aTCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCt                                                                                                                   <  1:217152/139‑1 (MQ=255)
                tCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                              <  2:656498/139‑1 (MQ=255)
                  gCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAACCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACg                                                                                                            <  1:538759/139‑1 (MQ=255)
                                  gggAATCCAGACGCGTGAAGCCTTCAAGTGTCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTAtt                                                                                                                      <  2:773707/113‑1 (MQ=255)
                                  gggAATCCAGACGCGTGAAGCCTTCAAGTGTCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTAtt                                                                                                                      >  1:773707/1‑113 (MQ=255)
                                            aCGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                               >  1:204986/1‑94 (MQ=255)
                                            aCGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                               <  2:204986/94‑1 (MQ=255)
                                                 tGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGAt                                                                             <  1:6690/139‑1 (MQ=255)
                                                        ttCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTAtt                                                                                                                      >  2:206876/1‑91 (MQ=255)
                                                        ttCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTAtt                                                                                                                      <  1:206876/91‑1 (MQ=255)
                                                        ttCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCAt                                                                      >  2:152253/1‑139 (MQ=255)
                                                                ccGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGa                                                              >  1:778069/1‑139 (MQ=255)
                                                                   gaCTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTCTTTACGGTCCTGACGACAGGACCGTTTTCCAACCGATTAATCATAAATATGaaaa                                                           >  2:57918/2‑139 (MQ=255)
                                                                             gATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAAttgttg                                                 <  1:459592/139‑1 (MQ=255)
                                                                              aTGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGc                                                >  1:698694/1‑139 (MQ=255)
                                                                                            aGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGAc                                                                                            >  1:547879/1‑81 (MQ=255)
                                                                                            aGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGAc                                                                                            <  2:547879/81‑1 (MQ=255)
                                                                                              ccAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCg                                <  1:338887/139‑1 (MQ=255)
                                                                                                  cccGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAAc                                                                                 >  1:715810/1‑86 (MQ=255)
                                                                                                  cccGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAAc                                                                                 <  2:715810/86‑1 (MQ=255)
                                                                                                   ccGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCt                           >  2:645208/1‑139 (MQ=255)
                                                                                                        aTTAAGTAAGCAGTTGATGGAATATACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATg                      <  1:467803/139‑1 (MQ=255)
                                                                                                           aagtaagCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGcac                   <  1:803719/139‑1 (MQ=255)
                                                                                                                agCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                              >  2:267571/1‑43 (MQ=255)
                                                                                                                agCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                              <  1:267571/43‑1 (MQ=255)
                                                                                                                agCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGa                                                                              <  1:140393/75‑1 (MQ=255)
                                                                                                                agCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGa                                                                              >  2:140393/1‑75 (MQ=255)
                                                                                                                agCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCaa              <  1:369479/139‑1 (MQ=255)
                                                                                                                 gCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAAc             >  1:380681/1‑139 (MQ=255)
                                                                                                                     ttGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCt                           <  2:318340/121‑1 (MQ=255)
                                                                                                                     ttGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCt                           >  1:318340/1‑121 (MQ=255)
                                                                                                                      tGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGttt        >  1:496097/1‑139 (MQ=255)
                                                                                                                       gATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTg       >  1:576488/1‑139 (MQ=255)
                                                                                                                       gATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTg       <  1:436556/139‑1 (MQ=255)
                                                                                                                          ggAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGAt                                                                             <  1:712094/66‑1 (MQ=255)
                                                                                                                          ggAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGAt                                                                             >  2:712094/1‑66 (MQ=255)
                                                                                                                          ggAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACg    <  1:152253/139‑1 (MQ=255)
                                                                                                                            aaTAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTa  <  2:380681/139‑1 (MQ=255)
                                                                                                                                |                                                                                                                                      
GTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTA  >  NC_000913/3719640‑3719902

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