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

AATAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCC  >  NC_000913/3100895‑3101138
                                                                                                                                           |                                                                                                           
aaTAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑c                                                                                                           <  2:185235/139‑1 (MQ=255)
    gagGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑ccccc                                                                                                       >  2:1340394/1‑139 (MQ=255)
                aaaTTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTccc                                                                                            <  2:1281307/139‑1 (MQ=255)
                 aaTTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCa                                                                                           <  2:664079/139‑1 (MQ=255)
                    tCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGca                                                                                        <  2:1321888/139‑1 (MQ=255)
                      ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGcaca                                                                                      >  1:265852/1‑139 (MQ=255)
                              ggCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTaa                                                                              >  1:216976/1‑139 (MQ=255)
                              ggCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTaa                                                                              >  2:1013508/1‑139 (MQ=255)
                                   tGCGGCTATCTTCCATAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑ccccccccagg                                                                                                 <  1:861945/114‑5 (MQ=255)
                                   tGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑ccccccccagg                                                                                                 >  2:861945/1‑110 (MQ=255)
                                   tGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑cccccccca                                                                                                   <  1:474143/112‑3 (MQ=255)
                                   tGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑cccccccca                                                                                                   >  2:474143/1‑110 (MQ=255)
                                           tCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                 <  1:640546/139‑1 (MQ=255)
                                           tCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTGTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                 <  1:1340394/139‑1 (MQ=255)
                                              tCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCtt                                                              <  2:265852/139‑1 (MQ=255)
                                                     tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTTCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCCAGGCTTCCCAGCACAtt                                                                                    <  2:1209740/111‑1 (MQ=255)
                                                     tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCCAGGCTTCCCAGCACAtt                                                                                    >  1:1209740/1‑111 (MQ=255)
                                                     tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACAt                                                                                     <  1:638604/109‑1 (MQ=255)
                                                     tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACAt                                                                                     >  2:638604/1‑109 (MQ=255)
                                                     tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCa                                                       >  2:436060/1‑139 (MQ=255)
                                                     tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCa                                                       >  1:1297173/1‑139 (MQ=255)
                                                       acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAcc                                                     >  2:568937/1‑139 (MQ=255)
                                                               gCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcgc                                             >  1:7170/1‑139 (MQ=255)
                                                                cTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcgcg                                            <  1:1013508/139‑1 (MQ=255)
                                                                cTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcgcg                                            <  2:216976/139‑1 (MQ=255)
                                                                     aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                       <  2:7170/139‑1 (MQ=255)
                                                                          cTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCAt                                  <  1:523318/139‑1 (MQ=255)
                                                                             tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                 >  1:597575/1‑105 (MQ=255)
                                                                             tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                 <  2:597575/105‑1 (MQ=255)
                                                                                       ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTg                     >  2:504916/1‑139 (MQ=255)
                                                                                         ggAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTa                   <  2:1033928/139‑1 (MQ=255)
                                                                                         ggAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTa                   <  1:568937/139‑1 (MQ=255)
                                                                                              tCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTccc                                                                                            >  2:1035225/1‑61 (MQ=255)
                                                                                              tCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTccc                                                                                            <  1:1035225/61‑1 (MQ=255)
                                                                                               cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc             <  1:436060/139‑1 (MQ=255)
                                                                                                     ccTGCTTTTTCCCCTCGGTTTTGATTGAGGATTGCACC‑‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCa       <  1:746410/139‑1 (MQ=255)
                                                                                                         cTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGc   <  1:504916/139‑1 (MQ=255)
                                                                                                          tttttCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGcc  <  2:1057772/139‑1 (MQ=255)
                                                                                                              tCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑cccccccca                                                                                                   >  2:287842/1‑35 (MQ=255)
                                                                                                              tCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑cccccccca                                                                                                   <  1:287842/37‑3 (MQ=255)
                                                                                                                                 ggATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTa                   >  2:353652/1‑99 (MQ=255)
                                                                                                                                 ggATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTa                   <  1:353652/99‑1 (MQ=255)
                                                                                                                                           |                                                                                                           
AATAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑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: ‑