Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,565,966 C→T *427* (TAG→TAA yjiN ← zinc‑type alcohol dehydrogenase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,565,9660CT97.1% 106.7 / ‑6.1 35*427* (TAG→TAAyjiNzinc‑type alcohol dehydrogenase‑like protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/16);  minor base G (1/0);  total (19/16)
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

TCGGTGACAAGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAAC  >  NC_000913/4565834‑4566089
                                                                                                                                    |                                                                                                                           
tCGGTGACAAGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTaaaaa                                                                                                                       <  1:256695/139‑1 (MQ=255)
tCGGTGACAAGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTaaaaa                                                                                                                       <  1:256696/139‑1 (MQ=255)
      aCAATTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCg                                                                                                                 <  1:508086/139‑1 (MQ=255)
         aGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGGTTTTTAACAATTGCTGAgg                                                                                                              >  2:408369/1‑139 (MQ=255)
          gTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAggg                                                                                                             >  1:364632/1‑139 (MQ=255)
                           ccccGTACCACGGGGGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGc                                                                                            >  1:437007/1‑139 (MQ=255)
                                  ccACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACa                                                                                     <  1:614214/139‑1 (MQ=255)
                                     cGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCa                                                                                  <  2:61730/139‑1 (MQ=255)
                                               tGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAAt                                                                        <  2:364632/139‑1 (MQ=255)
                                                 aaaaGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTa                                                                      >  2:312208/1‑139 (MQ=255)
                                                  aaaGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAg                                                                     >  2:653630/1‑139 (MQ=255)
                                                     gagGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGcc                                                                                                                  >  1:421419/1‑91 (MQ=255)
                                                     gagGCCACATTATTACCATTCTTTTAGTATCGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGcc                                                                                                                  <  2:421419/91‑1 (MQ=255)
                                                          cacaTTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATAc                                                             <  1:312208/139‑1 (MQ=255)
                                                                       ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAg                                                >  2:663046/1‑139 (MQ=255)
                                                                       ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAg                                                >  2:504115/1‑139 (MQ=255)
                                                                       ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAg                                                >  1:163412/1‑139 (MQ=255)
                                                                       ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAg                                                >  1:699334/1‑139 (MQ=255)
                                                                       ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAg                                                >  2:272192/1‑139 (MQ=255)
                                                                         cTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCg                                              >  1:199664/1‑139 (MQ=255)
                                                                         cTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCg                                              >  1:750930/1‑139 (MQ=255)
                                                                                tATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGtt                                       <  2:437007/139‑1 (MQ=255)
                                                                                  tGGTTTATTTGATCTCTCTCGCGGTTCGCCACTTAGTCTTTCATGGATTTTTAAAAATTGCTGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGa                                     <  2:490954/139‑1 (MQ=255)
                                                                                  tGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGa                                     <  1:663046/139‑1 (MQ=255)
                                                                                    gTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGaca                                   >  2:423371/1‑139 (MQ=255)
                                                                                     tttATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTGTCATGGATTTGTAAAAATTGCCGAGGCGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGacac                                  >  1:490954/1‑139 (MQ=255)
                                                                                          ttGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                             <  2:699334/139‑1 (MQ=255)
                                                                                              tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACt                         >  2:390285/1‑139 (MQ=255)
                                                                                                   ctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGa                    >  1:358528/1‑139 (MQ=255)
                                                                                                            aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg           <  1:653630/139‑1 (MQ=255)
                                                                                                              ccACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGat         >  1:442819/1‑139 (MQ=255)
                                                                                                               cACTTAGTTTTTCATGGGTTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGata        <  1:408369/139‑1 (MQ=255)
                                                                                                                     gTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTAGACACGAATAAACTGCAGATCTTTGCCGATATTTAAc  <  1:390285/139‑1 (MQ=255)
                                                                                                                           cATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCaa                                                    <  2:471627/83‑1 (MQ=255)
                                                                                                                           cATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCaa                                                    >  1:471627/1‑83 (MQ=255)
                                                                                                                                    |                                                                                                                           
TCGGTGACAAGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAAC  >  NC_000913/4565834‑4566089

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