Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,101,039 (C)7→8 coding (585/720 nt) yggN ← DUF2884 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,101,0321.C97.1% 110.9 / ‑3.2 34coding (592/720 nt)yggNDUF2884 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (13/20);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

AATAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCC  >  NC_000913/3100895‑3101138
                                                                                                                                           |                                                                                                           
aaTAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑c                                                                                                           <  2:185235/139‑1 (MQ=255)
    gagGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑ccccc                                                                                                       >  2:1340394/1‑139 (MQ=255)
                aaaTTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTccc                                                                                            <  2:1281307/139‑1 (MQ=255)
                 aaTTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCa                                                                                           <  2:664079/139‑1 (MQ=255)
                    tCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGca                                                                                        <  2:1321888/139‑1 (MQ=255)
                      ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGcaca                                                                                      >  1:265852/1‑139 (MQ=255)
                              ggCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTaa                                                                              >  1:216976/1‑139 (MQ=255)
                              ggCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTaa                                                                              >  2:1013508/1‑139 (MQ=255)
                                   tGCGGCTATCTTCCATAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑ccccccccagg                                                                                                 <  1:861945/114‑5 (MQ=255)
                                   tGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑ccccccccagg                                                                                                 >  2:861945/1‑110 (MQ=255)
                                   tGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cccccccca                                                                                                   <  1:474143/112‑3 (MQ=255)
                                   tGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cccccccca                                                                                                   >  2:474143/1‑110 (MQ=255)
                                           tCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                 <  1:640546/139‑1 (MQ=255)
                                           tCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTGTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                 <  1:1340394/139‑1 (MQ=255)
                                              tCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCtt                                                              <  2:265852/139‑1 (MQ=255)
                                                     tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTTCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCCAGGCTTCCCAGCACAtt                                                                                    <  2:1209740/111‑1 (MQ=255)
                                                     tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCCAGGCTTCCCAGCACAtt                                                                                    >  1:1209740/1‑111 (MQ=255)
                                                     tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACAt                                                                                     <  1:638604/109‑1 (MQ=255)
                                                     tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACAt                                                                                     >  2:638604/1‑109 (MQ=255)
                                                     tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCa                                                       >  2:436060/1‑139 (MQ=255)
                                                     tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCa                                                       >  1:1297173/1‑139 (MQ=255)
                                                       acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAcc                                                     >  2:568937/1‑139 (MQ=255)
                                                               gCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcgc                                             >  1:7170/1‑139 (MQ=255)
                                                                cTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcgcg                                            <  1:1013508/139‑1 (MQ=255)
                                                                cTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcgcg                                            <  2:216976/139‑1 (MQ=255)
                                                                     aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                       <  2:7170/139‑1 (MQ=255)
                                                                          cTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCAt                                  <  1:523318/139‑1 (MQ=255)
                                                                             tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                 >  1:597575/1‑105 (MQ=255)
                                                                             tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                 <  2:597575/105‑1 (MQ=255)
                                                                                       ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTg                     >  2:504916/1‑139 (MQ=255)
                                                                                         ggAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTa                   <  2:1033928/139‑1 (MQ=255)
                                                                                         ggAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTa                   <  1:568937/139‑1 (MQ=255)
                                                                                              tCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTccc                                                                                            >  2:1035225/1‑61 (MQ=255)
                                                                                              tCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTccc                                                                                            <  1:1035225/61‑1 (MQ=255)
                                                                                               cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc             <  1:436060/139‑1 (MQ=255)
                                                                                                     ccTGCTTTTTCCCCTCGGTTTTGATTGAGGATTGCACC‑‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCa       <  1:746410/139‑1 (MQ=255)
                                                                                                         cTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGc   <  1:504916/139‑1 (MQ=255)
                                                                                                          tttttCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGcc  <  2:1057772/139‑1 (MQ=255)
                                                                                                              tCCACTCGGTTTGGATTGAGGATTGCAG‑‑cccccccca                                                                                                   >  2:287842/1‑35 (MQ=255)
                                                                                                              tCCACTCGGTTTGGATTGAGGATTGCAG‑‑cccccccca                                                                                                   <  1:287842/37‑3 (MQ=255)
                                                                                                                                 ggATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTa                   >  2:353652/1‑99 (MQ=255)
                                                                                                                                 ggATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTa                   <  1:353652/99‑1 (MQ=255)
                                                                                                                                           |                                                                                                           
AATAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCC  >  NC_000913/3100895‑3101138

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