Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,979,525 T→C 8.2% intergenic (‑310/‑228) insA ← / → uspC IS1 repressor TnpA/universal stress protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,979,5250TC 8.2% 152.9 / 4.0 49intergenic (‑310/‑228)insA/uspCIS1 repressor TnpA/universal stress protein
Reads supporting (aligned to +/- strand):  ref base T (20/25);  new base C (2/2);  total (22/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 = 5.93e-01

TGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCAGACTAATA  >  NC_000913/1979379‑1979660
                                                                                                                                                  |                                                                                                                                       
tGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACttaatgttaa                                                                                                                                       <  4:24404/149‑1 (MQ=255)
     acaGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGaaa                                                                                                                                  >  4:465557/1‑149 (MQ=255)
      caGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAAt                                                                                                                                 <  1:52034/149‑1 (MQ=255)
          cAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGttt                                                                                                                             <  2:26841/149‑1 (MQ=255)
           aaaaTAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGtttt                                                                                                                            >  3:6291/1‑149 (MQ=255)
           aaaaTAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGtttt                                                                                                                            >  1:75577/1‑149 (MQ=255)
            aaaTAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTa                                                                                                                           <  1:220609/149‑1 (MQ=255)
            aaaTAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTa                                                                                                                           <  3:353578/149‑1 (MQ=255)
             aaTAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTAc                                                                                                                          >  2:117686/1‑149 (MQ=255)
               tAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACat                                                                                                                        <  4:372158/149‑1 (MQ=255)
                       cGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAAc                                                                                                                <  3:465557/149‑1 (MQ=255)
                            ataatCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACtt                                                                                                           >  2:371811/1‑149 (MQ=255)
                                tcatcTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAGTACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGttt                                                                                                       >  1:138244/1‑149 (MQ=255)
                                   tcTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTa                                                                                                                           >  3:300429/1‑126 (MQ=255)
                                   tcTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACCAATTAAATACAAAGAAAGGCAACAATCAACTTCATGTTAATGAAATGTTTTa                                                                                                                           <  4:300429/126‑1 (MQ=255)
                                    cTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATg                                                                                                   <  3:54349/149‑1 (MQ=255)
                                               aatCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAg                                                                                        <  3:260490/149‑1 (MQ=255)
                                                 tCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTa                                                                                      <  2:75577/149‑1 (MQ=255)
                                                     tttCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACt                                                                                  >  1:85240/1‑149 (MQ=255)
                                                           gCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTc                                                                            <  2:151478/149‑1 (MQ=255)
                                                            cAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAATAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCt                                                                           <  2:85240/149‑1 (MQ=255)
                                                                 ccAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATaa                                                                      <  1:309919/149‑1 (MQ=255)
                                                                                   aaCAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAAt                                                         >  4:245085/1‑144 (MQ=255)
                                                                                   aaCAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAAt                                                         <  3:245085/144‑1 (MQ=255)
                                                                                   aaCAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTatat                                                    <  1:117686/149‑1 (MQ=255)
                                                                                    aCAGACACTTAGCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAAAGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGCATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATg                                                   >  1:388060/1‑149 (MQ=255)
                                                                                      aGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATtat                                                      <  4:84255/144‑1 (MQ=255)
                                                                                      aGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATtat                                                      >  3:84255/1‑144 (MQ=255)
                                                                                          acTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTCAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCAcc                                             >  4:283451/1‑149 (MQ=255)
                                                                                          acTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTCAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCAcc                                             >  2:105299/1‑149 (MQ=255)
                                                                                                        aaaCTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTc                                                                            <  4:120624/104‑1 (MQ=255)
                                                                                                        aaaCTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTc                                                                            >  3:120624/1‑104 (MQ=255)
                                                                                                        aaaCTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTgattgatt                               >  1:321211/1‑149 (MQ=255)
                                                                                                         aaCTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGatta                              >  4:419958/1‑149 (MQ=255)
                                                                                                          aCTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTCAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGattaa                             <  3:283451/149‑1 (MQ=255)
                                                                                                          aCTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTCAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGattaa                             <  1:105299/149‑1 (MQ=255)
                                                                                                              aTTAAATACAAAGAAAGGCAACAATCAACTTAATGTTGATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTAt                         >  1:286374/1‑149 (MQ=255)
                                                                                                               ttAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTc                                                                            <  2:188954/97‑1 (MQ=255)
                                                                                                               ttAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTc                                                                            >  1:188954/1‑97 (MQ=255)
                                                                                                                    tACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACt                                                                                  >  1:458359/1‑86 (MQ=255)
                                                                                                                    tACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACt                                                                                  <  2:458359/86‑1 (MQ=255)
                                                                                                                    tACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTAc                   >  4:8130/1‑149 (MQ=255)
                                                                                                                     aCAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACt                  <  3:419958/149‑1 (MQ=255)
                                                                                                                          gaaagGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACtaataa             >  4:424825/1‑149 (MQ=255)
                                                                                                                            aagGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATc           <  2:354087/149‑1 (MQ=255)
                                                                                                                             agGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCa          <  2:388060/149‑1 (MQ=255)
                                                                                                                               gCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCAGa        <  4:6291/149‑1 (MQ=255)
                                                                                                                                ctacaaTCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCAGAc       <  2:138244/147‑1 (MQ=255)
                                                                                                                                    aaTCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCAGACTAat   <  2:286374/149‑1 (MQ=255)
                                                                                                                                     aTCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCAGACTAata  >  4:240168/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                       
TGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCAGACTAATA  >  NC_000913/1979379‑1979660

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