Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 585,217 C→T W139* (TGG→TGA)  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,2170CT100.0% 153.2 / NA 45W139* (TGG→TGA) ompTDLP12 prophage, outer membrane protease VII, outer membrane protein 3b
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/22);  total (23/22)

GGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCC  >  NC_000913/585083‑585351
                                                                                                                                      |                                                                                                                                      
gggAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGagtcagc                                                                                                                                    <  1:1167456/139‑1 (MQ=255)
   aaGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCtt                                                                                                                                 <  2:541007/139‑1 (MQ=255)
     ggAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTg                                                                                                                               >  2:651637/1‑139 (MQ=255)
                    tctGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAAATGGGTTCGTTGAGGAGTGAGCCTTTGATATTCAGATCAAAt                                                                                                                >  1:1090136/1‑139 (MQ=255)
                          tCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTg                                                                                                          >  2:1485993/1‑139 (MQ=255)
                             ctcctcAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCa                                                                                                       >  2:432432/1‑139 (MQ=255)
                              tcctcAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCAt                                                                                                      >  2:60599/1‑139 (MQ=255)
                                  cAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAAtt                                                                                                  <  1:60599/139‑1 (MQ=255)
                                  cAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAAtt                                                                                                  <  1:651637/139‑1 (MQ=255)
                                     aaCTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTgg                                                                                                         <  1:200687/129‑1 (MQ=255)
                                     aaCTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTgg                                                                                                         >  2:200687/1‑129 (MQ=255)
                                       cTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGtt                                                                                             <  2:1090136/139‑1 (MQ=255)
                                          tAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGagtcagc                                                                                                                                    >  1:1282868/1‑97 (MQ=255)
                                          tAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGagtcagc                                                                                                                                    <  2:1282868/97‑1 (MQ=255)
                                          tAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTtgtg                                                                                          >  1:1140415/1‑139 (MQ=255)
                                              tataGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATc                                                                                      <  1:250760/139‑1 (MQ=255)
                                                  ggAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGagtcagc                                                                                                                                    >  2:1061676/1‑89 (MQ=255)
                                                  ggAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGagtcagc                                                                                                                                    <  1:1061676/89‑1 (MQ=255)
                                                      ccaccTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGtt                                                                                             >  1:569898/1‑124 (MQ=255)
                                                      ccaccTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGtt                                                                                             <  2:569898/124‑1 (MQ=255)
                                                        accTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTa                                                                            >  2:334645/1‑139 (MQ=255)
                                                              ggCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTGTCa                                                                      <  1:432432/139‑1 (MQ=255)
                                                               gCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCAAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCAt                                                                     <  1:502511/139‑1 (MQ=255)
                                                                   tAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGt                                                                 >  1:1178707/1‑139 (MQ=255)
                                                                             aCGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATaa                                                                                                    <  2:1055382/94‑1 (MQ=255)
                                                                             aCGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATaa                                                                                                    >  1:1055382/1‑94 (MQ=255)
                                                                             aCGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGtt                                                                                             <  1:459556/101‑1 (MQ=255)
                                                                             aCGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGtt                                                                                             >  2:459556/1‑101 (MQ=255)
                                                                                  tttCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGtt                                                                                             >  1:512729/1‑96 (MQ=255)
                                                                                  tttCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGtt                                                                                             <  2:512729/96‑1 (MQ=255)
                                                                                  tttCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATc                                                                                      >  2:122365/1‑103 (MQ=255)
                                                                                  tttCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATc                                                                                      <  1:122365/103‑1 (MQ=255)
                                                                                     ccTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACt                                               >  2:266642/1‑139 (MQ=255)
                                                                                          tatCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAat                                          <  2:1140415/139‑1 (MQ=255)
                                                                                          tatCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAat                                          <  1:1485993/139‑1 (MQ=255)
                                                                                               ggCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAatccatc                                     <  1:334645/139‑1 (MQ=255)
                                                                                                gCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTgg                                                                                                         <  2:1182084/70‑1 (MQ=255)
                                                                                                gCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTgg                                                                                                         >  1:1182084/1‑70 (MQ=255)
                                                                                                    tGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGagtca                                                                                                                                      >  2:312055/1‑37 (MQ=255)
                                                                                                                 gTAATTGGGTTCGTTGAGGAGTCAGCCTTTGATATTGAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCata                   >  1:1084989/1‑139 (MQ=255)
                                                                                                                  tAATTGGGTTCGTTGAGGAGTCAGCCTTTGATAtt                                                                                                                          <  2:700899/35‑1 (MQ=255)
                                                                                                                  tAATTGGGTTCGTTGAGGAGTCAGCCTTTGATAtt                                                                                                                          >  1:700899/1‑35 (MQ=255)
                                                                                                                        ggTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGtt                                                                                             >  2:84595/1‑58 (MQ=255)
                                                                                                                        ggTTCGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGtt                                                                                             <  1:84595/58‑1 (MQ=255)
                                                                                                                            cGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACttt                                                                        >  2:647907/1‑75 (MQ=255)
                                                                                                                            cGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACttt                                                                        <  1:647907/75‑1 (MQ=255)
                                                                                                                            cGTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCg        >  2:1346426/1‑139 (MQ=255)
                                                                                                                             gTTGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCgg       <  1:1346426/139‑1 (MQ=255)
                                                                                                                              ttGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGtt                                                                                             >  2:1505515/1‑52 (MQ=255)
                                                                                                                              ttGAGGAGTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGtt                                                                                             <  1:1505515/52‑1 (MQ=255)
                                                                                                                                  ggagTCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGcc  >  2:1005528/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                      
GGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCC  >  NC_000913/585083‑585351

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