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)

GTACCTGCTACTAACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAG  >  NC_000913/986927‑987376
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                
gTACCTGCTACTAACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAatat                                                                                                                                                                                                         >  1:80376/1‑251 (MQ=255)
      gCTACTAACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCGTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTCTCTCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAAc                                                                                                                                                                                                    <  1:85143/250‑1 (MQ=255)
           aaaaagagAAGGGACGTTCACTGCAAGAATATTTCGTTTCATATTTATTTATGACCCTCATGGTTTTTGTTATTACGCCTGCCACTGCCGGCAATAAGTTCTGTCAATAAAAATTTACGTAACTTTTGATGATAGGTTGTGGTGTTTATGTGTCTGCAGGTATCTTTCCATTCAAACTAAGTTTTCGTTCCTTGTAAAGGGTTACAAAGATAAATAGTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAg                                                                                                                                                                                              <  1:67899/247‑1 (MQ=255)
                      ggACGATCACTGCCAGAATATTGCGCTTCATAATTATTTATTACCCTCATTGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCAGTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAAc                                                                                                                                                                                                    <  1:34686/234‑1 (MQ=255)
                                                   atcatTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCAt                                                                                                                                                      >  1:4582/1‑251 (MQ=255)
                                                     catTATTTATTCCCTCTATTGTTTTTTTTAGGACTCCTGGCACTGACGTCAATAAGTTCTGTCAATTAAAATTTACGGAACTATTGATGAAAGTGTGGTGTCTGTATTTGTCTGGAGGCATCGTTCCATTCAAACTAACGTTTCTCTACCATGAAAGTGCTACAAAGATAAATATTGGGTTTCAAAAATAAAAAATATTTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTGCCATag                                                                                                                                                    <  1:26217/251‑1 (MQ=255)
                                                                                         ncTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATg                                                                                                                                                                                                        >  1:102074/2‑164 (MQ=255)
                                                                                                           nngTNNNNNCAANAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATg                                                                                                                                                                                                        <  2:102074/144‑1 (MQ=255)
                                                                                                                  gTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATCTTTGGTTTCAAAAAGAAAAAATATGTAACCCAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGAAATTTTAATCTATTTATATGATTTCCTTATATTTAAAtt                                                                                        >  1:129371/1‑251 (MQ=255)
                                                                                                                                                gagTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGa                                                           >  1:141471/1‑250 (MQ=255)
                                                                                                                                                           cTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTGCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTc                                                                  <  1:151713/232‑1 (MQ=255)
                                                                                                                                                                                  nccATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTGCAAAGATAAATATTTGGTTTCAAAACGCAAAAATATGTAACCACAAGTAAAATTTAAGGAACTTTGTGAAAACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAAt                                                                                         >  1:15538/2‑186 (MQ=255)
                                                                                                                                                                                              aaCGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTGAATCTATTTATATGAGTTCCTTATATTTAAATTAACTAAACGGAAATTTTGGTTCTGATGTAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATgn            <  1:19352/251‑2 (MQ=255)
                                                                                                                                                                                                ngTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGcc           >  1:39802/2‑250 (MQ=255)
                                                                                                                                                                                                         aacgttagtttgaatggaaagataaatatTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAAGTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATTATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGAGGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAg  <  1:156110/234‑1 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                
GTACCTGCTACTAACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAG  >  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: