Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 987,537 A→C intergenic (‑555/+48) ompF ← / ← asnS outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,5370AC100.0% 151.9 / NA 47intergenic (‑555/+48)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (19/28);  total (19/28)

GCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCG  >  NC_000913/987389‑987684
                                                                                                                                                     |                                                                                                                                                   
gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGcaa                                                                                                                                                   <  2:874849/151‑1 (MQ=255)
   aTTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGcaaa                                                                                                                                                  <  1:160187/149‑1 (MQ=255)
          aaGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATaaaaa                                                                                                                                          <  2:523490/150‑1 (MQ=255)
              aaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATAGAATTGTTGACTTTAAAAACAGATAAATTTACAACAaaaaaaaaaaaaaa                                                                                                                                      >  1:846632/1‑148 (MQ=255)
               aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGc                                                                                                                                    <  2:356561/151‑1 (MQ=255)
               aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTTGTTGGATTTAAAAAAAAAAAAATTTACAGCAAAGAaaaaaaaaaa                                                                                                                                      >  2:134933/1‑147 (MQ=255)
                 aaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTTGTTGACTTAACAAAAAAAAAACTTGACAGCAAAGATaaaaaaaaa                                                                                                                                      >  2:531106/1‑145 (MQ=255)
                     aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGc                                                                                                                                <  1:420992/149‑1 (MQ=255)
                     aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGc                                                                                                                                <  1:29108/149‑1 (MQ=255)
                     aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAAAAGATTGGTTGACTTAACAAAAAAAAAAATTTACAGCAaaaaaaaaaaaaaa                                                                                                                                      >  1:908235/1‑141 (MQ=255)
                     aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATAAATTTCTTGTCTTAACAAAAAAATAAATTGACAGCAAAGATAAAAAAAAAAACAGCgg                                                                                                                              >  1:48641/1‑151 (MQ=255)
                         cGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGCAAAAAATTTCTTTACTTTACAAACAGAAAAAATTAAAGCAAAAAAAAAAAAAAAAACAGCggg                                                                                                                             >  1:276456/1‑148 (MQ=255)
                          gCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCt                                                                                                                         <  1:234/151‑1 (MQ=255)
                          gCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATGCTTTAATTTACAAAAAAATAAATTTACAACAaaaaaaaaaaaaaaacccacggg                                                                                                                             >  2:811646/1‑136 (MQ=255)
                           cAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTg                                                                                                                        <  2:96366/151‑1 (MQ=255)
                           cAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAGATTTCTTTGCTTTACAAACAAAAAACTTGACAGCAAAGATAAAAAAAAAAACAGCggggg                                                                                                                           >  1:30873/1‑146 (MQ=255)
                                            aGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTa                                                                                                        <  2:393019/150‑1 (MQ=255)
                                                 ttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTa                                                                                                   <  1:24783/150‑1 (MQ=255)
                                                     gcgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATGCTTGACTTAAAAAAAAAAAAAATTTACAGCAAAAAAAAAAAAAAAAACAACGGGGGCGGGGTCTTTGAGAAGCtaaaataaaa                                                                                                >  1:723040/1‑149 (MQ=255)
                                                       gTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGACATAAAATTTTTTCATTTACAAACAGAGAAATTGACAACACAGATAAAAAAAAAGCAAAGGGGGGCTGGGGCTTCTAGTCGCTTAATTAGAAACtc                                                                                            >  1:495146/1‑150 (MQ=255)
                                                         aaaTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTgg                                                                                           <  1:186054/150‑1 (MQ=255)
                                                          aaTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCg                                                                                         <  1:101748/151‑1 (MQ=255)
                                                             gTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTa                                                                                      <  2:398042/151‑1 (MQ=255)
                                                                      tAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGTAAAAAAAATTGTTTGCTTTAAAAAAAAATAAAATGATAACAAAAATAAAAAAAAAAACAACGGGCGCCGGGGCTTTTAGAAGCGAAAATAAAAAGAGGCGGTACGCgg                                                                                 >  2:396102/1‑147 (MQ=255)
                                                                        atatTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACg                                                                           <  1:664921/151‑1 (MQ=255)
                                                                          atTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTg                                                                         <  1:26595/151‑1 (MQ=255)
                                                                              cGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGaa                                                                     <  2:880756/151‑1 (MQ=255)
                                                                                 aaCAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACgg                                                                  <  1:813138/151‑1 (MQ=255)
                                                                                      gatgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatca                                                             <  1:795680/151‑1 (MQ=255)
                                                                                          atCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatc                                                         <  2:195724/151‑1 (MQ=255)
                                                                                                aTGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtac                                                   <  1:37493/151‑1 (MQ=255)
                                                                                                  ggCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtacgt                                                 <  2:90051/151‑1 (MQ=255)
                                                                                                   gCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGtt                                                <  1:148463/151‑1 (MQ=255)
                                                                                                          tGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc                                          <  1:905970/150‑1 (MQ=255)
                                                                                                              aaaTAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGcc                                       >  2:44631/1‑150 (MQ=255)
                                                                                                                aTAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGt                                    <  2:81657/151‑1 (MQ=255)
                                                                                                                 tAGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGtt                                   >  1:868247/1‑151 (MQ=255)
                                                                                                                  aGATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTa                                  <  2:18646/151‑1 (MQ=255)
                                                                                                                   gATTGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTAc                                 <  1:606520/151‑1 (MQ=255)
                                                                                                                      tGCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGt                               <  2:868797/150‑1 (MQ=255)
                                                                                                                       gCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTaa                             <  2:592822/151‑1 (MQ=255)
                                                                                                                       gCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTa                              >  1:556381/1‑150 (MQ=255)
                                                                                                                       gCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTa                              >  1:556743/1‑150 (MQ=255)
                                                                                                                       gCTTGACTTAACAAACAGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTa                              >  2:869719/1‑150 (MQ=255)
                                                                                                                                    aacaGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGtt                <  1:120871/151‑1 (MQ=255)
                                                                                                                                      caGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTc               >  1:65252/1‑150 (MQ=255)
                                                                                                                                         aTAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCaaa            >  2:908478/1‑150 (MQ=255)
                                                                                                                                           aaCTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAAcc         >  1:318056/1‑151 (MQ=255)
                                                                                                                                                  cAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCg  <  1:1090/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCG  >  NC_000913/987389‑987684

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