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% 51.2 / NA 16intergenic (‑170/+433)ompF/asnSouter membrane porin F/asparagine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (4/12);  total (4/12)

TAGATTTCTGCAGCGTTTGCAGTACCTGCTACTAACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAAT  >  NC_000913/986906‑987385
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
taGATTTCTGCAGCGTTTGCAGTACCTGCTACTAACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTCTTGATGCGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACAGTGAAAGTGCTACAAAGATAAATAttt                                                                                                                                                                                                                                       >  1:334227/1‑251 (MQ=255)
            gCGTTTGCAGTACCTGCTACTAACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAg                                                                                                                                                                                                                           <  1:15012/251‑1 (MQ=255)
                           gCTACTAACAGAGCAGGGACGATCCCTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAAcc                                                                                                                                                                                                            <  1:353640/251‑1 (MQ=255)
                                             acgatcgtttccAGAATATTGCGCTTCATCATTATTTATTACCTTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTGTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGa                                                                                                                                                                                          <  2:36426/241‑1 (MQ=255)
                                                                       cataattattTATTACCCTCATGGTTTTTTTTAGGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCa                                                                                                                                                                <  1:177925/247‑1 (MQ=255)
                                                                           attattAATTACCATAGTGGTTTATTTTATGACCCCTGCCACTGCGGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACTTATCGCTTCCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATGTAAGGAACTTTGTGAACACCTTCATATTTCCATaga                                                                                                                                                            <  2:348600/251‑1 (MQ=255)
                                                                              atttatttCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGAc                                                                                                                                                         <  1:38072/251‑1 (MQ=255)
                                                                                       ccTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTgatgat                                                                                                                                                 >  2:36792/1‑250 (MQ=255)
                                                                                          cATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTGCCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATAtt                                                                                                                                              <  1:159775/250‑1 (MQ=255)
                                                                                                         tacacCTGACACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTAGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATGTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGAGATTTACAGCAATTTTAAt                                                                                                                               <  2:260954/249‑1 (MQ=255)
                                                                                                             ccTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGGGATGATATTTACAGCAATTTTAATCTAtt                                                                                                                          <  1:372579/251‑1 (MQ=255)
                                                                                                                                                    aCGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAAt                                                                                   <  2:116795/251‑1 (MQ=255)
                                                                                                                                                                              tctttatgtgaatgagcgcctccttCCATTAAAACTAACGTTTAGCCACCGTGAACGTGCTACAAAGATACATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACAGTCATATTTCCATAGAGACGTGATGAAATTTACAGCAATATTAATCTATTTATATGATTTCCTTATATTTAACTTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCg                                                          >  2:42721/24‑250 (MQ=255)
                                                                                                                                                                                tttATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGAcc                                                       >  1:168179/1‑251 (MQ=255)
                                                                                                                                                                                           gCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACaaa                                              <  2:200700/249‑1 (MQ=255)
                                                                                                                                                                                                                                     aGTGCTACAAAGAGAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAAAACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAAt  <  2:28554/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
TAGATTTCTGCAGCGTTTGCAGTACCTGCTACTAACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAAT  >  NC_000913/986906‑987385

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