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.C90.0% 92.7 / 3.7 31coding (592/720 nt)yggNDUF2884 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (12/15);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

AATAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGGCG  >  NC_000913/3100895‑3101165
                                                                                                                                          |                                                                                                                                     
aaTAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGc                                                                                                                                      <  1:371135/139‑1 (MQ=255)
      ggATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGccccccc                                                                                                                                >  1:472275/1‑139 (MQ=255)
         ttattTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAgg                                                                                                                             <  2:1052531/139‑1 (MQ=255)
          tattTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGccccccccagg                                                                                                                            >  2:454997/1‑135 (MQ=255)
          tattTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGccccccccagg                                                                                                                            <  1:840602/139‑5 (MQ=255)
                    tCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGca                                                                                                                   <  2:295713/139‑1 (MQ=255)
                                   tGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGt                                                                                                    >  2:99548/1‑139 (MQ=255)
                                   tGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGcccccccca                                                                                                                              >  1:132382/1‑110 (MQ=255)
                                   tGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGcccccccca                                                                                                                              <  2:132382/112‑3 (MQ=255)
                                          aTCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGcac                                                                                                                  >  2:395489/1‑118 (MQ=255)
                                          aTCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGcac                                                                                                                  <  1:395489/118‑1 (MQ=255)
                                                     tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCa                                                                                  >  2:749673/1‑139 (MQ=255)
                                                                cTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcgcg                                                                       <  2:5712/139‑1 (MQ=255)
                                                                  aCAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGTATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGcc                                                                     <  2:472275/139‑1 (MQ=255)
                                                                     aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  1:454997/139‑1 (MQ=255)
                                                                     aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  1:512286/139‑1 (MQ=255)
                                                                     aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  2:39430/139‑1 (MQ=255)
                                                                     aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  1:760684/139‑1 (MQ=255)
                                                                            ttGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCttt                                                                           <  1:208974/123‑1 (MQ=255)
                                                                            ttGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCttt                                                                           >  2:208974/1‑123 (MQ=255)
                                                                              gCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGcaca                                                                                                                 <  1:70716/83‑1 (MQ=255)
                                                                              gCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGcaca                                                                                                                 >  2:70716/1‑83 (MQ=255)
                                                                                                         cTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                                            <  2:657735/77‑1 (MQ=255)
                                                                                                         cTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                                            >  1:657735/1‑77 (MQ=255)
                                                                                                          tttttCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGcc                             <  1:643336/139‑1 (MQ=255)
                                                                                                              tCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTcc                                       >  1:140963/1‑125 (MQ=255)
                                                                                                              tCCACTCGGTTTGGATTGAGGATTGCAGACCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTcc                                       <  2:140963/125‑1 (MQ=255)
                                                                                                                   tCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGc                                                                                                                    >  1:581438/1‑43 (MQ=255)
                                                                                                                   tCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGc                                                                                                                    <  2:581438/43‑1 (MQ=255)
                                                                                                                   tCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACt                    >  1:34902/1‑139 (MQ=255)
                                                                                                                          ggATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCt             <  1:907995/139‑1 (MQ=255)
                                                                                                                                     tGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTc                                                   <  2:21288/90‑1 (MQ=255)
                                                                                                                                     tGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTc                                                   >  1:21288/1‑90 (MQ=255)
                                                                                                                                     tGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGgcg  >  2:211423/1‑139 (MQ=255)
                                                                                                                                     tGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGgcg  >  1:436336/1‑139 (MQ=255)
                                                                                                                                     tGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGcc                                                                     >  2:295727/1‑71 (MQ=255)
                                                                                                                                     tGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGcc                                                                     <  1:295727/71‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                     
AATAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGGCG  >  NC_000913/3100895‑3101165

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: