Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 585,420 Δ1 bp coding (214/954 nt) ompT ← DLP12 prophage, outer membrane protease VII, outer membrane protein 3b

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913585,4200T.100.0% 61.0 / NA 15coding (214/954 nt)ompTDLP12 prophage, outer membrane protease VII, outer membrane protein 3b
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (7/8);  total (7/8)

GGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTTTGTTTTTCCGCTCAGAGTTCCAAGACTAATGTCCGCATTTATGTTGTCAGGAGTA  >  NC_000913/585299‑585553
                                                                                                                         |                                                                                                                                     
gggTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAA‑AATTGCAGCGTTATTGaa                                                                                                                     >  1:165794/1‑139 (MQ=255)
               atccaGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAA‑AATTGCAGCGTTATTGAATTTCCAGTCGAGTtg                                                                                                      >  2:402729/1‑139 (MQ=255)
                                     tGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAA‑AATTGCAGCGTTATTGAATTTCCAGTCGAg                                                                                                         <  1:321837/114‑1 (MQ=255)
                                     tGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAA‑AATTGCAGCGTTATTGAATTTCCAGTCGAg                                                                                                         >  2:321837/1‑114 (MQ=255)
                                              ggCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAA‑AATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGa                                                                       >  2:49521/1‑139 (MQ=255)
                                                              ccagccagCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAA‑AATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTTTg                                                       <  2:393157/139‑1 (MQ=255)
                                                                    agcagcCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAA‑AATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCgg                                                                                         <  1:368566/99‑1 (MQ=255)
                                                                    agcagcCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAA‑AATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCgg                                                                                         >  2:368566/1‑99 (MQ=255)
                                                                         ccccGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAA‑AATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTTTGTTTTTCCGCTc                                            <  1:402729/139‑1 (MQ=255)
                                                                                       tGGGGCATCAAATCCCAATTAATTGCACCTTTAA‑AATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTTTGTTTTTCCGCTCAGAGTTCCAAGACt                              <  1:49521/139‑1 (MQ=255)
                                                                                            cATCAAATCCCAATTAATTGCACCTTTAA‑AATTGCAGTGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACTCGCTCTTTTGTTTTTCCGCTCAGAGTTCCAAGACTAATGt                         >  1:305143/1‑139 (MQ=255)
                                                                                            cATCAAATCCCAATTAATTGCACCTTTAA‑AATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTTTGTTTTTCCGCTCAGAGTTCCa                                   >  1:178444/1‑129 (MQ=255)
                                                                                            cATCAAATCCCAATTAATTGCACCTTTAA‑AATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTTTGTTTTTCCGCTCAGAGTTCCa                                   <  2:178444/129‑1 (MQ=255)
                                                                                             aTCAAATCCCAATTAATTGCACCTTTAA‑AATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACCTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTTTGTTTTTCCGCTCAGAGTTCCAAGACTAATGTc                        <  1:17196/139‑1 (MQ=255)
                                                                                                                   cTTTAA‑AATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTTTGTTTTTCCGCTCAGAGTTCCAAGACTAATGTCCGCATTTATGTTGTCAGGAGTa  <  2:305143/139‑1 (MQ=255)
                                                                                                                         |                                                                                                                                     
GGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTTTGTTTTTCCGCTCAGAGTTCCAAGACTAATGTCCGCATTTATGTTGTCAGGAGTA  >  NC_000913/585299‑585553

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