Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,197,739 A→G S167P (TCC→CCC)  ymfE ← e14 prophage; putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,197,7390AG96.9% 209.3 / ‑2.2 65S167P (TCC→CCC) ymfEe14 prophage; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (35/28);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

CTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTT  >  NC_000913/1197643‑1197833
                                                                                                 |                                                                                              
ctGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAAt                                                                                             >  1:974103/1‑100 (MQ=255)
ctGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATACATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAAt                                                                                             >  1:118592/1‑100 (MQ=255)
 tGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAAtt                                                                                            >  2:1276269/1‑100 (MQ=255)
    gaATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAAt                                                                                         <  2:854865/100‑1 (MQ=255)
      ataaTAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGc                                                                                       <  2:427122/97‑1 (MQ=255)
      aTGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGc                                                                                       <  1:1052729/100‑1 (MQ=255)
       tGATAAAATTACAAAAAAGAGTAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCa                                                                                      >  2:969082/1‑100 (MQ=255)
        gATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCaa                                                                                     <  1:339872/100‑1 (MQ=255)
         aTAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAAt                                                                                    <  1:886281/100‑1 (MQ=255)
         aTAAAATAACAAAAAAGAATAAAGACAGAAAAGAAACTTCAACATGAAAAATATAAAAAAAAAAAAAAAAAAAGATTATTAGGAAGGGAATTAATGCAAt                                                                                    <  2:388986/100‑1 (MQ=255)
           aaaaTTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTa                                                                                  <  2:368687/100‑1 (MQ=255)
             aaTTACAAAAAAGAATAAATACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATc                                                                                <  2:567987/100‑1 (MQ=255)
              aTTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCg                                                                               >  1:1097617/1‑100 (MQ=255)
                 aCAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaaa                                                                            >  2:1396344/1‑100 (MQ=255)
                  cAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaaaa                                                                           <  1:1019712/100‑1 (MQ=255)
                  cAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaaaa                                                                           >  1:376387/1‑100 (MQ=255)
                   aaaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAt                                                                          <  1:210662/100‑1 (MQ=255)
                    aaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAtt                                                                         >  2:500498/1‑100 (MQ=255)
                     aaaaaaATAAAAACAGCAAAGAACCTTCAACATGAAAAATATCCAAATAATTAAAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAtt                                                                         <  2:1390520/99‑1 (MQ=255)
                     aaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTc                                                                        >  2:476505/1‑100 (MQ=255)
                     aaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTc                                                                        <  2:504287/100‑1 (MQ=255)
                     aaaaGAATAAAGACAGCAAAGAAACTTCAACATGCAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTc                                                                        <  2:834735/100‑1 (MQ=255)
                      aaaGAATAACGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGATAATTCa                                                                       <  1:456878/100‑1 (MQ=255)
                         gAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGTGAATTAATGCAATTATCGAAAATTCaaaa                                                                    >  2:70666/1‑100 (MQ=255)
                           aTAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGCGTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaaaaa                                                                  >  2:1499843/1‑100 (MQ=255)
                           aTAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaaaca                                                                  >  1:21107/1‑98 (MQ=255)
                           aTAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaaaaa                                                                  >  1:1299032/1‑100 (MQ=255)
                            tAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAat                                                                 >  1:583854/1‑100 (MQ=255)
                             aaaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAata                                                                <  1:227097/100‑1 (MQ=255)
                             aaaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAata                                                                >  2:1019301/1‑100 (MQ=255)
                                  cAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaa                                                           <  2:478830/100‑1 (MQ=255)
                                    gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaa                                                         >  1:250016/1‑100 (MQ=255)
                                     cAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaaa                                                        >  1:1535242/1‑100 (MQ=255)
                                      aaaGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAAt                                                       >  2:1272337/1‑100 (MQ=255)
                                       aaGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATa                                                      <  1:463453/100‑1 (MQ=255)
                                        aGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAg                                                     <  1:1341640/100‑1 (MQ=255)
                                         gAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGt                                                    <  2:144423/100‑1 (MQ=255)
                                                     tGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATACAAAAAATATCCAAAAATAGTATACTTTATTcc                                        <  2:271464/100‑1 (MQ=255)
                                                      gAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCa                                       >  2:872878/1‑100 (MQ=255)
                                                      gAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCa                                       >  1:473557/1‑100 (MQ=255)
                                                       aaaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCag                                      >  2:1029650/1‑100 (MQ=255)
                                                       aaaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCag                                      <  1:278755/100‑1 (MQ=255)
                                                          aaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCagaag                                   <  1:1243691/100‑1 (MQ=255)
                                                            tatCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCagaaga                                  >  2:348169/1‑99 (MQ=255)
                                                                cATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCa                             >  2:1366376/1‑100 (MQ=255)
                                                                cATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATCATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAAAAGAGTTCa                             >  1:1031298/1‑100 (MQ=255)
                                                                  tttgtttgCAAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCaa                            <  2:1172394/100‑1 (MQ=255)
                                                                  aaaaaaaaaaaaaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCaa                            <  2:1491840/90‑1 (MQ=255)
                                                                     gtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAAtata                        >  1:491196/1‑100 (MQ=255)
                                                                        tgCAAAAAAAGATTATTAGGTAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATg                     >  2:661594/1‑100 (MQ=255)
                                                                         gCAAAAAAAGATTATTAGGAAGGGAATTTATGCAATTATCGAATATTCAAACAATATCTTAAAATATTACACTTTATTCCAGAAGATTCCAATATAAt                      >  2:1492622/1‑98 (MQ=255)
                                                                         gCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                    >  2:266289/1‑100 (MQ=255)
                                                                         gCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                    >  2:861079/1‑100 (MQ=255)
                                                                          cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                   >  2:340954/1‑100 (MQ=255)
                                                                          cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                   >  1:1004581/1‑100 (MQ=255)
                                                                          cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                   >  1:1387669/1‑100 (MQ=255)
                                                                          cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                   <  2:1103572/100‑1 (MQ=255)
                                                                          cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                   >  2:469766/1‑100 (MQ=255)
                                                                          cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAAAGTTCAATATAATGtt                   >  2:112653/1‑100 (MQ=255)
                                                                             aaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTg                 <  1:1276274/100‑1 (MQ=255)
                                                                             aaaaaaGATTATTAAAAAGGGAATTAATGCAATTATCAAAAATAAAAAAAATATCAAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTg                 <  2:725409/100‑1 (MQ=255)
                                                                                    attattAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCaa          >  2:1215733/1‑100 (MQ=255)
                                                                                    attattAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCaa          >  1:871291/1‑100 (MQ=255)
                                                                                         taGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCATAAGAGTTCAATATAATGTTTGTCTTCAAttttt     <  2:486506/100‑1 (MQ=255)
                                                                                         taGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttttt     >  1:770194/1‑100 (MQ=255)
                                                                                          aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc    <  2:1229295/100‑1 (MQ=255)
                                                                                          aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc    >  2:688211/1‑100 (MQ=255)
                                                                                          aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc    <  1:688213/100‑1 (MQ=255)
                                                                                            gaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCtt  <  1:1274417/100‑1 (MQ=255)
                                                                                            gaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCtt  <  2:408907/100‑1 (MQ=255)
                                                                                            gaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCtt  <  2:408159/100‑1 (MQ=255)
                                                                                            gaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCtt  <  1:1362072/100‑1 (MQ=255)
                                                                                                 |                                                                                              
CTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTT  >  NC_000913/1197643‑1197833

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