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,2280CT100.0% 105.5 / NA 34*274* (TAG→TAAyfcODUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/18);  total (16/18)

GTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTGACCGAGAC  >  NC_000913/2449096‑2449359
                                                                                                                                    |                                                                                                                                   
gTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAAc                                                                                                                               >  2:986698/1‑139 (MQ=255)
    ggcgcgtttgttgTGGTGAATGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCAtt                                                                                                                          <  1:850756/132‑1 (MQ=255)
                        ttgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAAc                                                                                                                               <  2:838428/115‑1 (MQ=255)
                        ttgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAAc                                                                                                                               >  1:838428/1‑115 (MQ=255)
                        ttgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGt                                                                                                                      <  1:243825/124‑1 (MQ=255)
                        ttgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGt                                                                                                                      >  2:243825/1‑124 (MQ=255)
                                     gTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGa                                                                                          >  2:1197344/1‑139 (MQ=255)
                                           gccgccGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCt                                                                                    <  1:612380/139‑1 (MQ=255)
                                                       gCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGt                                                                                                                      <  1:1263032/93‑1 (MQ=255)
                                                       gCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGt                                                                                                                      >  2:1263032/1‑93 (MQ=255)
                                                        cccTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTa                                                                                                                     <  2:952528/93‑1 (MQ=255)
                                                        cccTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTa                                                                                                                     >  1:952528/1‑93 (MQ=255)
                                                                 gTTGAGCGCGTTGCGCCAGGTACGGGTAACGGCAGGGGCGAATGCTCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGGTTGGCTGCATACTTCTGGTATCAAgacg                                                              <  2:921995/139‑1 (MQ=255)
                                                                   tGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgca                                                            <  1:1197344/139‑1 (MQ=255)
                                                                       cgcgTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCg                                                                                                        >  1:1076746/1‑91 (MQ=255)
                                                                       cgcgTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCg                                                                                                        <  2:1076746/91‑1 (MQ=255)
                                                                           ttGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCata                                                    <  1:1416115/139‑1 (MQ=255)
                                                                                                    gggCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTgg                                                                                      >  1:134766/1‑80 (MQ=255)
                                                                                                    gggCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTgg                                                                                      <  2:134766/80‑1 (MQ=255)
                                                                                                    gggCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAg                           >  2:470080/1‑139 (MQ=255)
                                                                                                             ccccTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGAt                  <  1:470080/139‑1 (MQ=255)
                                                                                                              cccTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATa                                                                               >  2:39529/1‑77 (MQ=255)
                                                                                                              cccTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATa                                                                               <  1:39529/77‑1 (MQ=255)
                                                                                                                ctgcctgGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCa                                                                                 >  2:919694/1‑73 (MQ=255)
                                                                                                                ctgcctgGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCa                                                                                 <  1:919694/73‑1 (MQ=255)
                                                                                                                ctgcctgGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATa                                                                               >  1:858943/1‑75 (MQ=255)
                                                                                                                ctgcctgGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATa                                                                               <  2:858943/75‑1 (MQ=255)
                                                                                                                       gtgtATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgca                                                            >  2:765824/1‑87 (MQ=255)
                                                                                                                       gtgtATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgca                                                            <  1:765824/87‑1 (MQ=255)
                                                                                                                            tCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATGctg                                                                            >  1:849703/1‑66 (MQ=255)
                                                                                                                            tCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATActg                                                                            <  2:849703/66‑1 (MQ=255)
                                                                                                                            tCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGc             >  1:600526/1‑129 (MQ=255)
                                                                                                                            tCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGc             <  2:600526/129‑1 (MQ=255)
                                                                                                                             cTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTGACCGAGAc  >  2:1397136/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                   
GTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTGACCGAGAC  >  NC_000913/2449096‑2449359

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: