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

CTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTT  >  NC_000913/1197643‑1197833
                                                                                                 |                                                                                              
ctGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAAt                                                                                             >  1:974103/1‑100 (MQ=255)
ctGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATACATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAAt                                                                                             >  1:118592/1‑100 (MQ=255)
 tGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAAtt                                                                                            >  2:1276269/1‑100 (MQ=255)
    gaATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAAt                                                                                         <  2:854865/100‑1 (MQ=255)
      ataaTAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGc                                                                                       <  2:427122/97‑1 (MQ=255)
      aTGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGc                                                                                       <  1:1052729/100‑1 (MQ=255)
       tGATAAAATTACAAAAAAGAGTAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCa                                                                                      >  2:969082/1‑100 (MQ=255)
        gATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCaa                                                                                     <  1:339872/100‑1 (MQ=255)
         aTAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAAt                                                                                    <  1:886281/100‑1 (MQ=255)
         aTAAAATAACAAAAAAGAATAAAGACAGAAAAGAAACTTCAACATGAAAAATATAAAAAAAAAAAA‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAAt                                                                                    <  2:388986/100‑1 (MQ=255)
           aaaaTTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTa                                                                                  <  2:368687/100‑1 (MQ=255)
             aaTTACAAAAAAGAATAAATACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATc                                                                                <  2:567987/100‑1 (MQ=255)
              aTTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCg                                                                               >  1:1097617/1‑100 (MQ=255)
                 aCAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaaa                                                                            >  2:1396344/1‑100 (MQ=255)
                  cAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaaaa                                                                           <  1:1019712/100‑1 (MQ=255)
                  cAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaaaa                                                                           >  1:376387/1‑100 (MQ=255)
                   aaaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAt                                                                          <  1:210662/100‑1 (MQ=255)
                    aaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAtt                                                                         >  2:500498/1‑100 (MQ=255)
                     aaaaaaATAAAAACAGCAAAGAACCTTCAACATGAAAAATATCCAAATAATTAA‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAtt                                                                         <  2:1390520/99‑1 (MQ=255)
                     aaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTc                                                                        >  2:476505/1‑100 (MQ=255)
                     aaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTc                                                                        <  2:504287/100‑1 (MQ=255)
                     aaaaGAATAAAGACAGCAAAGAAACTTCAACATGCAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTc                                                                        <  2:834735/100‑1 (MQ=255)
                      aaaGAATAACGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGATAATTCa                                                                       <  1:456878/100‑1 (MQ=255)
                         gAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGTGAATTAATGCAATTATCGAAAATTCaaaa                                                                    >  2:70666/1‑100 (MQ=255)
                           aTAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGCGTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaaaaa                                                                  >  2:1499843/1‑100 (MQ=255)
                           aTAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaaaca                                                                  >  1:21107/1‑98 (MQ=255)
                           aTAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaaaaa                                                                  >  1:1299032/1‑100 (MQ=255)
                            tAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAat                                                                 >  1:583854/1‑100 (MQ=255)
                             aaaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAata                                                                <  1:227097/100‑1 (MQ=255)
                             aaaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAata                                                                >  2:1019301/1‑100 (MQ=255)
                                  cAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaa                                                           <  2:478830/100‑1 (MQ=255)
                                    gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaa                                                         >  1:250016/1‑100 (MQ=255)
                                     cAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaaa                                                        >  1:1535242/1‑100 (MQ=255)
                                      aaaGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAAt                                                       >  2:1272337/1‑100 (MQ=255)
                                       aaGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATa                                                      <  1:463453/100‑1 (MQ=255)
                                        aGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAg                                                     <  1:1341640/100‑1 (MQ=255)
                                         gAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGt                                                    <  2:144423/100‑1 (MQ=255)
                                                     tGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATACAAAAAATATCCAAAAATAGTATACTTTATTcc                                        <  2:271464/100‑1 (MQ=255)
                                                      gAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCa                                       >  2:872878/1‑100 (MQ=255)
                                                      gAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCa                                       >  1:473557/1‑100 (MQ=255)
                                                       aaaaaTATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCag                                      >  2:1029650/1‑100 (MQ=255)
                                                       aaaaaTATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCag                                      <  1:278755/100‑1 (MQ=255)
                                                          aaTATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCagaag                                   <  1:1243691/100‑1 (MQ=255)
                                                            tatCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCagaaga                                  >  2:348169/1‑99 (MQ=255)
                                                                cATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCa                             >  2:1366376/1‑100 (MQ=255)
                                                                cATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATCATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAAAAGAGTTCa                             >  1:1031298/1‑100 (MQ=255)
                                                                  tttgtttgCAAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCaa                            <  2:1172394/100‑1 (MQ=255)
                                                                  aaaaaaaaaaaaaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCaa                            <  2:1491840/90‑1 (MQ=255)
                                                                     gtttgC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAAtata                        >  1:491196/1‑100 (MQ=255)
                                                                        tgC‑AAAAAAAGATTATTAGGTAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATg                     >  2:661594/1‑100 (MQ=255)
                                                                         gC‑AAAAAAAGATTATTAGGAAGGGAATTTATGCAATTATCGAATATTCAAACAATATCTTAAAATATTACACTTTATTCCAGAAGATTCCAATATAAt                      >  2:1492622/1‑98 (MQ=255)
                                                                         gC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                    >  2:266289/1‑100 (MQ=255)
                                                                         gC‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                    >  2:861079/1‑100 (MQ=255)
                                                                          c‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                   >  2:340954/1‑100 (MQ=255)
                                                                          c‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                   >  1:1004581/1‑100 (MQ=255)
                                                                          c‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                   >  1:1387669/1‑100 (MQ=255)
                                                                          c‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                   <  2:1103572/100‑1 (MQ=255)
                                                                          c‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                   >  2:469766/1‑100 (MQ=255)
                                                                          c‑AAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAAAGTTCAATATAATGtt                   >  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)
                                                                                                 |                                                                                              
CTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGC‑AAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTT  >  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: ‑