Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 987,152 G→T intergenic (‑170/+433) ompF ← / ← asnS outer membrane porin F/asparagine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,1520GT100.0% 36.0 / NA 14intergenic (‑170/+433)ompF/asnSouter membrane porin F/asparagine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (7/7);  total (7/7)

GTACCTGCTACTAACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTC‑ATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTG‑AAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAG  >  NC_000913/986927‑987376
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                
gTACCTGCTACTAACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTC‑ATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTG‑AAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAatat                                                                                                                                                                                                         >  1:80376/1‑251 (MQ=255)
      gCTACTAACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTC‑ATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCGTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTCTCTCTACCGTG‑AAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAAc                                                                                                                                                                                                    <  1:85143/250‑1 (MQ=255)
           aaaaagagAAGGGACGTTCACTGCAAGAATATTTCGTTTCATATTTATTTATGACCCTC‑ATGGTTTTTGTTATTACGCCTGCCACTGCCGGCAATAAGTTCTGTCAATAAAAATTTACGTAACTTTTGATGATAGGTTGTGGTGTTTATGTGTCTGCAGGTATCTTTCCATTCAAACTAAGTTTTCG‑TTCCTTGTAAAGGGTTACAAAGATAAATAGTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAg                                                                                                                                                                                              <  1:67899/247‑1 (MQ=255)
                      ggACGATCACTGCCAGAATATTGCGCTTCATAATTATTTATTACCCTC‑ATTGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCAGTCAAACTAACGTTTCGCTACCGTG‑AAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAAc                                                                                                                                                                                                    <  1:34686/234‑1 (MQ=255)
                                                   atcatTATTTATTACCCTC‑ATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTG‑AAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCAt                                                                                                                                                      >  1:4582/1‑251 (MQ=255)
                                                     catTATTTATT‑CCCTCTATTGTTTTTTTTAGGACTCCTGGCACTGACGTCAATAAGTTCTGTCAATTAAAATTTACGGAACTATTGATGAAAGTGTGGTGTCTGTATTTGTCTGGAGGCATCGTTCCATTCAAACTAACGTTTCTCTACCATG‑AAAGTGCTACAAAGATAAATATTGGGTTTCAAAAATAAAAAATATTTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTGCCATag                                                                                                                                                    <  1:26217/251‑1 (MQ=255)
                                                                                         ncTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTG‑AAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATg                                                                                                                                                                                                        >  1:102074/2‑164 (MQ=255)
                                                                                                           nngTNNNNNCAANAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTG‑AAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATg                                                                                                                                                                                                        <  2:102074/144‑1 (MQ=255)
                                                                                                                  gTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTG‑AAAGTGCTACAAAGATAAATCTTTGGTTTCAAAAAGAAAAAATATGTAACCCAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGAAATTTTAATCTATTTATATGATTTCCTTATATTTAAAtt                                                                                        >  1:129371/1‑251 (MQ=255)
                                                                                                                                                gagTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTG‑AAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGa                                                           >  1:141471/1‑250 (MQ=255)
                                                                                                                                                           cTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTG‑AAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTGCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTc                                                                  <  1:151713/232‑1 (MQ=255)
                                                                                                                                                                                  nccATTCAAACTAACGTTTCGCTACCGTG‑AAAGTGCTGCAAAGATAAATATTTGGTTTCAAAACGCAAAAATATGTAACCACAAGTAAAATTTAAGGAACTTTGTGAAAACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAAt                                                                                         >  1:15538/2‑186 (MQ=255)
                                                                                                                                                                                              aaCGTTTCGCTACCGTG‑AAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTGAATCTATTTATATGAGTTCCTTATATTTAAATTAACTAAACGGAAATTTTGGTTCTGATGTAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATgn            <  1:19352/251‑2 (MQ=255)
                                                                                                                                                                                                ngTTTCGCTACCGTG‑AAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGcc           >  1:39802/2‑250 (MQ=255)
                                                                                                                                                                                                         aacgttagtttgaatggaaagataaatatTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAAGTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATTATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGAGGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAg  <  1:156110/234‑1 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                
GTACCTGCTACTAACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTC‑ATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTG‑AAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAG  >  NC_000913/986927‑987376

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