Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,892,414 A→C 26.7% intergenic (+428/‑92) PP_RS08755 → / → napA tRNA‑Asp/Periplasmic nitrate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,892,4140AC26.7% 28.3 / 2.7 34intergenic (+428/‑92)PP_RS08755/napAtRNA‑Asp/Periplasmic nitrate reductase
Reads supporting (aligned to +/- strand):  ref base A (15/9);  new base C (7/3);  total (22/12)
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

GCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGCCCCTACTGCG  >  ACB122/1892293‑1892545
                                                                                                                         |                                                                                                                                    
gcgcCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCAGAAATTGCTGATCTACCTCATCCCGGGCTGGCAACGTCGCCAGCCCTCACCCCGTCACGAc                                                                                                                          >  1:96840/1‑134 (MQ=39)
    cgccCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCGGGGAACGGGGCCCGCCCTCCCCCCGTCACGACAAAgg                                                                                                                     >  7:54806/1‑135 (MQ=37)
      cccTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGCTCTACCTCATCCCGGGGTGTCCACGTCGCCAGCCCGCACCACTTCACGACAAAAgcc                                                                                                                   >  1:23367/1‑134 (MQ=39)
         tCAGGCCTTGCGATATCCCTGTGGGAGCGGGCCTGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCGGATCTACCTCCTACCGGGCGGTCCACGTCGACAGCCCGCCCCACTGCACGACAAAGGCGccc                                                                                                                >  8:170251/1‑134 (MQ=255)
                    gATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACATCATCCCGGGGTGTCAACCTCGAACGCCCTCACCACGTCAAGACAAAAGCGCCGTGTTTCCCCGc                                                                                                     >  8:24161/1‑135 (MQ=38)
                          ccTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTGAACGTCGCCAGCCCTCACCCCTTCACCACAAAGGCGCCGGGTTGCCCCCCTTGCa                                                                                                >  4:299332/1‑134 (MQ=39)
                          ccTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCGGGCAACGTCGACCGCCCTCACCACTTCACGACAAAAGCCCCGGGTTTGCCCGCTTgga                                                                                                >  7:178496/1‑132 (MQ=38)
                                 gAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCGACGTCGACAGCCCACCCCACATCACCACAAAGGCGCCCTGTTGCCCCGCTTGGATCCGCCTg                                                                                        >  7:169766/1‑135 (MQ=255)
                                 gAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGGGAACGTCGCCCGCCCGCACCCCTTCACGACAAAGGCGCCGTGTTGCCCCGCGTGGCTGAGCAtt                                                                                        >  6:99419/1‑134 (MQ=255)
                                      ggCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCGGTGAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGGATCAGCATGCCggg                                                                                   >  5:120424/1‑135 (MQ=255)
                                       gCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCgggg                                                                                  >  7:143551/1‑135 (MQ=255)
                                       gCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACATCATCCCGGGCTGTCAACGTCGACAGCCCGCCCCACTAAAAGACAAAGGCGCCGTGTTTGCCCGCTTGGTTAAGCATGCCcggg                                                                                  >  6:88600/1‑135 (MQ=255)
                                         aTGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGGTGTGAACGTCGCCCGCCCGCACCCCTTCACGACAAAGGCGCCGGGTTGCCCCCGGTGGATCACCATGCCGGGGcc                                                                                >  2:284461/1‑134 (MQ=255)
                                         aTGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACATCACGACAAAGGCGCCGTGTTTCCCCCCGTGCATCAGCAAGCCGGGGCa                                                                                >  2:33071/1‑135 (MQ=255)
                                                   aaGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTAAACCTCGACAGCCCGCCCCACATCACGACAAACGCGCCGTGTTGCCCCGCTTGCATCAGCAAGCCGGGGCAGCCGGGGccg                                                                      >  3:20599/1‑134 (MQ=255)
                                                       gCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGGCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTGGCATCAGCATGCCGGGGCAGCCGGCGCCttttt                                                                 <  8:220395/135‑1 (MQ=255)
                                                              taGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGCCAGCCCTCCCCCCTTCACGCCAAAGCCGCCGTGTTTCCCCGCTTGGCTCCTCATGCCGGGGCCGCCGGCGCCGTTTTTCGgtt                                                           >  6:8021/1‑134 (MQ=255)
                                                                    gCCCATAACTTGCCGAGCTCCCTCACCCCGGGCGGGCAACGTCGACCGCCCTCCCCACTTCACGACAAAGGCGCCGTGGTGCCCCGCTGGCATCAGCATGCGGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAg                                                    <  7:292585/135‑1 (MQ=255)
                                                                    gCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGCCAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTTCCCCGCGTGCATCAGCATGCCGGGGCAGCCGGCGCCGTTTTTCGGTTTGCACAg                                                    >  1:172512/1‑135 (MQ=255)
                                                                        aTAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGCCAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCCGCCCGCGCCGTTTTTGGGTTTTCACAGgcgc                                                >  6:86012/1‑135 (MQ=255)
                                                                        aTAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATAAGCATGCCGGGGCAGCCGGCGCCGTTTTGCGGTTTGCCCAGgcgc                                                >  8:89189/1‑135 (MQ=255)
                                                                        aTAAATTGCTGATCTACATCATCCCGGGCTGTAAACGTCGACAGCCCTGACCACTTCACCACAAAGGCGCAGGGTTGCCACCCATGGATCAACATGCCGGGGCACCCGGCCGCATTTTTCCTTTTACACCGGAgc                                                >  2:158343/1‑135 (MQ=255)
                                                                         tAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCCCCACTTCACGACAAACGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCCGCCGGCGCCCTTTTGCGTTTTACACACGCgcc                                               >  8:150911/1‑135 (MQ=255)
                                                                                 tGATCTACCTCACCCCGGGCTGTCAACGTCGACCGCCCTCACCACTTCACCACAAAGGCGCGGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCAt                                        <  7:150911/134‑1 (MQ=255)
                                                                                    tCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCCCCCCTTCACGAAAACGGCGCGGTGTGGCCCCGCTTGCATCTGCATGCGGGGGCAGCGGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGc                                     <  7:24161/134‑1 (MQ=255)
                                                                                               ccGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAACGGCGCCGTGTTGCCCCGCTTGCATCAGCAGGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGt                         <  7:20102/135‑1 (MQ=255)
                                                                                                      gacaccGTCGAAAGCCCTCACCACTTAACGAAAATGGCGGCTTGTTGCCCCGCTTGCATCAGCTTGGGGGGGCAGCGGGGGCTTTTTTGCGTTTTCCACTGGAGCCCGCCCATGGCAAACAGCTAAGTACGCAGc                  <  6:96447/130‑1 (MQ=255)
                                                                                                        cAACGGCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCAGCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGc                  <  8:159898/133‑1 (MQ=255)
                                                                                                         aaCGTCGACAGCCCTCACCACTTCACGACAACGGCGCCGTGTTGCCCCGCTGGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTc               <  1:3873/135‑1 (MQ=255)
                                                                                                          aCGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGc                                                                                                     >  7:271381/1‑49 (MQ=255)
                                                                                                          aCGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGc                                                                                                     <  8:271381/49‑1 (MQ=255)
                                                                                                            ggcGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTGGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGc            <  5:86012/133‑1 (MQ=255)
                                                                                                              ggACCGCCCTCCCCACTTCACCCCAAGGGCCCCGGGTTGCCCCGCTTGCATCAGCATGGGGGGGCAGCGGGCGCCTTTTTGCGAGTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGccc          <  3:117793/134‑1 (MQ=255)
                                                                                                              cGACCGCCCCCACCAATCCCCGACAACGGCGCCGTGTTGCCCCGCTGGCTGACGCAGGCGGGGTCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGccc          <  3:202474/135‑1 (MQ=255)
                                                                                                                      cgcaccccTTCCCGACAACGGCGCCGGGTTCCCCCGATTGCAGCAGCATGCCGGGGCAGCGGGGGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGCCCCTACTgcg  <  7:72615/133‑1 (MQ=255)
                                                                                                                         |                                                                                                                                    
GCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGCCCCTACTGCG  >  ACB122/1892293‑1892545

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