Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,449,228 C→T *274* (TAG→TAA yfcO ← DUF2544 family putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,449,2280CT96.5% 86.7 / ‑6.2 29*274* (TAG→TAAyfcODUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/13);  minor base G (0/1);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTGACCGAGA  >  NC_000913/2449094‑2449358
                                                                                                                                      |                                                                                                                                  
cGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACa                                                                                                                                >  2:449232/1‑139 (MQ=255)
    gATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTc                                                                                                                            <  1:887757/139‑1 (MQ=255)
        cggcgTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTc                                                                                                                        <  2:1094791/139‑1 (MQ=255)
                                       gTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGa                                                                                         >  2:391382/1‑139 (MQ=255)
                                        tGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGAt                                                                                        <  1:465815/139‑1 (MQ=255)
                                          gCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTg                                                                                      <  1:391382/139‑1 (MQ=255)
                                                   ggAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAAc                                                                                                                              <  2:632924/90‑1 (MQ=255)
                                                   ggAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAAc                                                                                                                              >  1:632924/1‑90 (MQ=255)
                                                                 gagTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAaca                                                               >  2:451332/1‑139 (MQ=255)
                                                                   gTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacg                                                             <  1:449232/139‑1 (MQ=255)
                                                                         cgcgTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGAc                                                       >  2:832505/1‑139 (MQ=255)
                                                                           cgTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCa                                                     <  1:832505/139‑1 (MQ=255)
                                                                             ttGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCata                                                   <  1:451332/139‑1 (MQ=255)
                                                                                     ggTACGGGTAATGGCAGGAGCGAATGCCCCTGCCTGGTGTATCTAAGGGGTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATActgc                                                                          <  1:63090/108‑1 (MQ=255)
                                                                                     ggTACGGGTAATGGCAGAGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATActgc                                                                          >  2:63090/1‑108 (MQ=255)
                                                                                               aTGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCAt                                                                                                                          >  1:596266/1‑50 (MQ=255)
                                                                                               aTGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCAt                                                                                                                          <  2:596266/50‑1 (MQ=255)
                                                                                                 ggCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTa                                            <  1:292196/126‑1 (MQ=255)
                                                                                                 ggCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTa                                            >  2:292196/1‑126 (MQ=255)
                                                                                                     ggggCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgca                                                           >  1:352361/1‑107 (MQ=255)
                                                                                                     ggggCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgca                                                           <  2:352361/107‑1 (MQ=255)
                                                                                                      gggCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCg                                                                                                       <  2:1010820/62‑1 (MQ=255)
                                                                                                      gggCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCg                                                                                                       >  1:1010820/1‑62 (MQ=255)
                                                                                                      gggCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTtataa                          >  2:230870/1‑138 (MQ=255)
                                                                                                      gggCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAg                          >  1:354323/1‑139 (MQ=255)
                                                                                                      gggCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAg                          >  1:163760/1‑139 (MQ=255)
                                                                                                                 cctgcctgGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGAttt               >  1:334191/1‑139 (MQ=255)
                                                                                                                              tCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTGACCgaga  >  2:782361/1‑139 (MQ=255)
                                                                                                                                   gggTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATa                                                                              <  2:676434/58‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                  
CGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTGACCGAGA  >  NC_000913/2449094‑2449358

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