Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,655,347 C→T *109* (TAG→TAA ynfA ← UPF0060 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,655,3470CT100.0% 121.2 / NA 38*109* (TAG→TAAynfAUPF0060 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/20);  total (18/20)

CGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACG  >  NC_000913/1655217‑1655479
                                                                                                                                  |                                                                                                                                    
cGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTgcg                                                                                                                              >  2:161764/1‑139 (MQ=255)
         ccGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGccc                                                                                                                    <  1:669391/139‑1 (MQ=255)
             aGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaa                                                                                                                <  2:309705/139‑1 (MQ=255)
               tGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCATCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaaca                                                                                                              <  2:572607/139‑1 (MQ=255)
               tGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaaca                                                                                                              <  2:318055/139‑1 (MQ=255)
                 aTCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaaca                                                                                                              <  2:361660/137‑1 (MQ=255)
                 aTCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaaca                                                                                                              >  1:361660/1‑137 (MQ=255)
                  tCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaa                                                                                                                <  1:245901/134‑1 (MQ=255)
                  tCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaa                                                                                                                >  2:245901/1‑134 (MQ=255)
                              gCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCg                                                                                               <  1:161764/139‑1 (MQ=255)
                                                          aaaaGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGt                                                                                                                                 >  1:53486/1‑78 (MQ=255)
                                                          aaaaGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGt                                                                                                                                 <  2:53486/78‑1 (MQ=255)
                                                               aTGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATaa                                                              >  2:204918/1‑139 (MQ=255)
                                                               aTGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATaa                                                              >  1:206085/1‑139 (MQ=255)
                                                                tGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATaaa                                                             <  2:443924/139‑1 (MQ=255)
                                                                        aaTACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGt                                                     >  1:304687/1‑139 (MQ=255)
                                                                               tttGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGc                                              >  1:266037/1‑139 (MQ=255)
                                                                                ttGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCa                                                                                      <  1:336731/98‑1 (MQ=255)
                                                                                ttGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCa                                                                                      >  2:336731/1‑98 (MQ=255)
                                                                                          acacTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAAcg                                   >  1:517310/1‑139 (MQ=255)
                                                                                                ccTCTTATCGTATTGATAAAATGGATTTATGTTCTTGCGCGCGCCACAGCCCGCAACAATGATCAACATGCCGTAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCAGCCACAACGCGCAGc                             <  2:526093/139‑1 (MQ=255)
                                                                                                 ctctTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGcc                            >  2:406738/1‑139 (MQ=255)
                                                                                                    ttatcttatTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACg                                               <  2:457086/117‑1 (MQ=255)
                                                                                                    ttatcttatTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACg                                               >  1:457086/1‑117 (MQ=255)
                                                                                                      atcttatTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaaca                                                                                                              <  1:329301/52‑1 (MQ=255)
                                                                                                      atcttatTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGcaaca                                                                                                              >  2:329301/1‑52 (MQ=255)
                                                                                                          tatTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACg                   <  1:293838/139‑1 (MQ=255)
                                                                                                           atTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCAt                                                                >  1:266171/1‑93 (MQ=255)
                                                                                                           atTGATAAAATGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCAt                                                                <  2:266171/93‑1 (MQ=255)
                                                                                                                 aaaaTGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCa            >  1:418496/1‑139 (MQ=255)
                                                                                                                   aaTGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATTAACGCACCCGTCCAGTCATAAAGAGTCAGTTTc                                                  <  1:249219/99‑1 (MQ=255)
                                                                                                                   aaTGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATTAACGCACCCGTCCAGTCATAAAGAGTCAGTTTc                                                  >  2:249219/1‑99 (MQ=255)
                                                                                                                   aaTGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTg              <  2:69741/135‑1 (MQ=255)
                                                                                                                   aaTGGATTTATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTg              >  1:69741/1‑135 (MQ=255)
                                                                                                                           tATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACg  <  2:451487/139‑1 (MQ=255)
                                                                                                                           tATGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACg  <  1:204918/139‑1 (MQ=255)
                                                                                                                             tGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTc                                                  <  2:207497/89‑1 (MQ=255)
                                                                                                                             tGTTCTTACGTGCGCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTc                                                  >  1:207497/1‑89 (MQ=255)
                                                                                                                                  |                                                                                                                                    
CGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACG  >  NC_000913/1655217‑1655479

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: