Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 749,989 A→G G727G (GGT→GGC ybgQ ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913749,9890AG100.0% 113.0 / NA 30G727G (GGT→GGCybgQputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/15);  total (15/15)

GCCAACCATACCCACGCTAATGCCGCTGTCATCCTGGCGGATATCCGCACCGAGGGGAGGGAACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAACCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTCATGTTCACCGCCACGGTGG  >  NC_000913/749860‑750080
                                                                                                                                 |                                                                                           
gCCAACCTTACCCACGCTAATGCCGCTGTCATCCTGGCGGATATCCGCACCGAGGGGAGGGAACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATcgcg                                                                                    <  2:542277/139‑1 (MQ=255)
    aCCATACCCACGCTAATGCCGCTGTCATCCTGGCGGATATCCGCACCGAGGGGAGGGAACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCtt                                                                                <  1:183484/139‑1 (MQ=255)
     ccATACCCACGCTAATGCCGCTGTCATCCTGGCGGATATCCGCACCGAGGGGAGGGAACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTc                                                                               >  1:14410/1‑139 (MQ=255)
               gCTAATGCCGCTGTCATCCTGGCGGATATCCGCACCGAGGGGAGGGAACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGtt                                                                     >  2:52719/1‑139 (MQ=255)
               gCTAATGCCGCTGTCATCCTGGCGGATATCCGCACCGAGGGGAGGGAACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGtt                                                                     >  2:509083/1‑139 (MQ=255)
               gCTAATGCCGCTGTCATCCTGGCGGATATCCGCACCGAGGGGAGGGAACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGtt                                                                     >  1:473191/1‑139 (MQ=255)
                cTAATGCCGCTGTCATCCTGGCGGATATCCGCACCGAGGGGAGGGAACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGttt                                                                    >  2:252461/1‑139 (MQ=255)
                   aTGCCGCTGTCATCCTGGCGGATATCCGCACCGAGGGGAGGGAACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCt                                                                 >  2:176927/1‑139 (MQ=255)
                                   ggCGGATATCCGCACCGAGGGGAGGGAACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGc                                                 <  2:473191/139‑1 (MQ=255)
                                      ggATATCCGCACCGAGGGGAGGGAACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAAc                                              >  2:118393/1‑139 (MQ=255)
                                          atCCGCACCGAGGGGAGGGAACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGtt                                          <  1:52719/139‑1 (MQ=255)
                                                aCCGAGGGGAGGGAACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCg                                    <  2:506953/139‑1 (MQ=255)
                                                    aggggagggAACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCgg                                <  1:509083/139‑1 (MQ=255)
                                                    aggggagggAACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCgg                                >  2:227345/1‑139 (MQ=255)
                                                     ggggagggAACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCggg                               <  2:14410/139‑1 (MQ=255)
                                                          gggaACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGt                          <  2:552801/139‑1 (MQ=255)
                                                            gaACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGtca                        >  1:87565/1‑139 (MQ=255)
                                                                  aCCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCaa                                                                             <  2:98384/80‑1 (MQ=255)
                                                                  aCCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCaa                                                                             >  1:98384/1‑80 (MQ=255)
                                                                  aCCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTCATg                  >  2:285/1‑139 (MQ=255)
                                                                    cGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTCATGtt                >  1:2637/1‑139 (MQ=255)
                                                                       tGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTCATGTTCAc             >  2:332804/1‑139 (MQ=255)
                                                                          cGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTCATGTTCACCGc          <  1:285/139‑1 (MQ=255)
                                                                           gTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTCATGTTCACCGcc         <  1:227345/139‑1 (MQ=255)
                                                                            ttGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTCATGTTCACCGCCa        <  2:2637/139‑1 (MQ=255)
                                                                               cGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGt                          >  2:297381/1‑118 (MQ=255)
                                                                               cGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGt                          <  1:297381/118‑1 (MQ=255)
                                                                                  atgatgACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTCATGTTCACCGCCACGGTgg  <  1:118393/139‑1 (MQ=255)
                                                                                              aCGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGtt                                          >  2:535140/1‑87 (MQ=255)
                                                                                              aCGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAGCCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGtt                                          <  1:535140/87‑1 (MQ=255)
                                                                                                                                 |                                                                                           
GCCAACCATACCCACGCTAATGCCGCTGTCATCCTGGCGGATATCCGCACCGAGGGGAGGGAACTGACCGCTGGCGTTGCGTATGATGACGTTAACGTCTTTACCGGAACGGGAAGCCAGTGATTTGTAACCTATCGCGCCTTCAATCCATGTTTCCTTGATGACGTTTTCTGCAACTGTTACGCCGTCGGGTAAGTCATTCATGTTCACCGCCACGGTGG  >  NC_000913/749860‑750080

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: