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.C93.9% 161.9 / 2.8 49coding (592/720 nt)yggNDUF2884 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (21/25);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.63e-01

GATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGGCG  >  NC_000913/3100902‑3101165
                                                                                                                                   |                                                                                                                                     
gATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTCCACTCGGTTTGGATTGAGGATTGCAGcccccccca                                                                                                                              >  2:514902/1‑137 (MQ=255)
     ttttAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCt                                                                                                                           >  1:520210/1‑139 (MQ=255)
     ttttAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCt                                                                                                                           >  2:409196/1‑139 (MQ=255)
       ttAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTcc                                                                                                                        <  1:514902/139‑1 (MQ=255)
            ttCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGc                                                                                                                    <  2:1341686/139‑1 (MQ=255)
            ttCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGc                                                                                                                    <  1:442249/139‑1 (MQ=255)
             tCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGca                                                                                                                   <  1:409196/139‑1 (MQ=255)
               ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGcaca                                                                                                                 >  1:1061762/1‑139 (MQ=255)
               ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGcaca                                                                                                                 >  1:416468/1‑139 (MQ=255)
                 gACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACAtt                                                                                                               <  2:1061762/139‑1 (MQ=255)
                       ggCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTc                                                                                                                         >  2:1519745/1‑123 (MQ=255)
                       ggCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTc                                                                                                                         <  1:1519745/123‑1 (MQ=255)
                                 cTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGATATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGc                                                                                                                                      >  2:1630751/1‑99 (MQ=255)
                                 cTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGcaca                                                                                                                 >  1:737567/1‑121 (MQ=255)
                                 cTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGcaca                                                                                                                 <  2:737567/121‑1 (MQ=255)
                                 cTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGc                                                                                                                                      <  1:1630751/99‑1 (MQ=255)
                                    tCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                                            <  2:520210/139‑1 (MQ=255)
                                              tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCa                                                                                  >  2:1136779/1‑139 (MQ=255)
                                              tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCa                                                                                  <  1:1339894/139‑1 (MQ=255)
                                                     gcgGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCttt                                                                           <  2:578339/139‑1 (MQ=255)
                                                      cgGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTc                                                                          >  1:1248996/1‑139 (MQ=255)
                                                      cgGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTc                                                                          >  1:1194191/1‑139 (MQ=255)
                                                       ggCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGcccccc                                                                                                                                 >  1:1326587/1‑82 (MQ=255)
                                                       ggCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGcccccc                                                                                                                                 <  2:1326587/82‑1 (MQ=255)
                                                         cTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcgcg                                                                       <  2:1248996/139‑1 (MQ=255)
                                                          tACAAATATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGccccccccagg                                                                                                                            >  1:952518/1‑80 (MQ=255)
                                                          tACAAATATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGccccccccagg                                                                                                                            <  2:952518/84‑5 (MQ=255)
                                                              aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  1:1359305/139‑1 (MQ=255)
                                                              aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  2:416468/139‑1 (MQ=255)
                                                                 aTCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTTGGTTTGGATTGAGGATTGCAGccccccccagg                                                                                                                            <  2:1005529/77‑5 (MQ=255)
                                                                 aTCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTACAGccccccccagg                                                                                                                            >  1:1005529/1‑73 (MQ=255)
                                                                   cTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACAtt                                                                                                               >  2:1558463/1‑89 (MQ=255)
                                                                   cTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACAtt                                                                                                               <  1:1558463/89‑1 (MQ=255)
                                                                   cTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCAt                                                             >  1:979176/1‑139 (MQ=255)
                                                                      tGCCAAACTGCTGGAAATCTTTTTCCTGGTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTaa                                                          >  1:46735/1‑139 (MQ=255)
                                                                        ccAAACTGCTGGAAATCTTTTTCCTGGTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATg                                                        <  2:46735/139‑1 (MQ=255)
                                                                         cAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACAtt                                                                                                               <  2:1072075/83‑1 (MQ=255)
                                                                         cAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACAtt                                                                                                               >  1:1072075/1‑83 (MQ=255)
                                                                                ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTg                                                >  2:1373136/1‑139 (MQ=255)
                                                                                ctgGAAATCCTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTGATTAATGCTGTCCTg                                                >  1:1134129/1‑139 (MQ=255)
                                                                                          ttttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  >  2:1144410/1‑111 (MQ=255)
                                                                                          ttttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  1:1144410/111‑1 (MQ=255)
                                                                                                tGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCAtt                                <  2:1194191/139‑1 (MQ=255)
                                                                                                  cTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATCGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAGTTCCACCCATTGc                              <  2:1134129/139‑1 (MQ=255)
                                                                                                       tCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        <  1:20427/124‑1 (MQ=255)
                                                                                                       tCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        >  2:20427/1‑124 (MQ=255)
                                                                                                            tCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGTGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACt                    <  2:979176/139‑1 (MQ=255)
                                                                                                            tCGGTTTGGATTGAGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACt                    >  1:399920/1‑139 (MQ=255)
                                                                                                                        gAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCgcgc                                                                      <  1:1087643/76‑1 (MQ=255)
                                                                                                                        gAGCATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcgcga                                                                      >  2:1087643/1‑75 (MQ=255)
                                                                                                                         aGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCa                                      >  2:714048/1‑108 (MQ=255)
                                                                                                                         aGGATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCa                                      <  1:714048/108‑1 (MQ=255)
                                                                                                                         aGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTc      >  2:152956/1‑139 (MQ=255)
                                                                                                                          ggATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  2:338514/79‑1 (MQ=255)
                                                                                                                          ggATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  >  1:338514/1‑79 (MQ=255)
                                                                                                                          ggATTGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTc      >  1:435100/1‑139 (MQ=255)
                                                                                                                              tGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGgcg  <  1:1373136/139‑1 (MQ=255)
                                                                                                                              tGCAGCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGgcg  <  1:1136779/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                     
GATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGGCG  >  NC_000913/3100902‑3101165

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