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,7390AG93.0% 171.8 / 3.6 57S167P (TCC→CCC) ymfEe14 prophage; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base G (32/21);  total (33/24)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01

TGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTA  >  NC_000913/1197644‑1197834
                                                                                               |                                                                                               
tGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAAtt                                                                                             >  1:84070/1‑100 (MQ=255)
 gagaATGATAAAA‑TACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTaa                                                                                           >  1:482028/1‑100 (MQ=255)
   gaATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAAt                                                                                          <  1:657915/100‑1 (MQ=255)
    aaTGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATACATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATg                                                                                         >  2:1063938/1‑100 (MQ=255)
         tAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAAtt                                                                                    >  2:495232/1‑100 (MQ=255)
          aaaaTTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTa                                                                                   <  2:374056/100‑1 (MQ=255)
            aaTTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATc                                                                                 <  1:533600/100‑1 (MQ=255)
             aTTACAAAAAAGAATAAAGACAGCAAAGAACATTCAACATGAAAAATATCCATTTAATAGAAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCg                                                                                <  2:876881/100‑1 (MQ=255)
                aCAAAAAAGAATAAAGACATCAAAGAACCTTCAACATGAAAAATATCCATTTTTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaa                                                                              <  1:1088335/99‑1 (MQ=255)
                  aaaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAt                                                                           <  1:127176/100‑1 (MQ=255)
                   aaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAACATATCCATTTGTTTGCAAAAAAACATAATTAGGAAGGGAATAAATGCACTTATCGAAAAtt                                                                          >  2:1107692/1‑100 (MQ=255)
                   aaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAtt                                                                          >  2:808065/1‑100 (MQ=255)
                    aaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTc                                                                         >  1:410257/1‑100 (MQ=255)
                      aaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaa                                                                       >  1:836290/1‑100 (MQ=255)
                       aGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaa                                                                      <  2:506137/100‑1 (MQ=255)
                         aaTAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaaaa                                                                    >  1:778479/1‑100 (MQ=255)
                         aaTAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaaaa                                                                    <  1:639174/100‑1 (MQ=255)
                          aTAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaaaaa                                                                   <  1:105143/100‑1 (MQ=255)
                           tAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAat                                                                  <  2:684477/100‑1 (MQ=255)
                           tAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAat                                                                  <  1:1038999/100‑1 (MQ=255)
                             aaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAatat                                                                >  2:20195/1‑100 (MQ=255)
                             aaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAatat                                                                >  2:742699/1‑100 (MQ=255)
                             aaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAatat                                                                >  1:1049541/1‑100 (MQ=255)
                              aGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATc                                                               >  1:928602/1‑100 (MQ=255)
                              aGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATc                                                               >  1:356041/1‑100 (MQ=255)
                                acaaaaaaGAAACTTAAACATGAAAAATATACAAAAAAAAAAAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATcc                                                              <  2:1014122/94‑1 (MQ=255)
                                aaaaaaaaGAACCTTCAACATGAAAAATATCCATTTGTTTACAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTAAAAAAAAAAAAAAAAATATcc                                                              <  2:1002604/94‑1 (MQ=255)
                                 cAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaa                                                            >  1:336149/1‑100 (MQ=255)
                                 cAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaa                                                            >  2:1097913/1‑100 (MQ=255)
                                   gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaa                                                          >  1:1034817/1‑100 (MQ=255)
                                     aaaGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAAt                                                        >  1:801795/1‑100 (MQ=255)
                                       aGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAg                                                      <  1:1219331/100‑1 (MQ=255)
                                           ccTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGtata                                                  >  1:1156250/1‑100 (MQ=255)
                                           ccTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAAGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGtat                                                   <  1:228234/99‑1 (MQ=255)
                                             ttCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACt                                                <  1:113271/100‑1 (MQ=255)
                                                aaCATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTa                                             >  2:895165/1‑100 (MQ=255)
                                                   aTGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTc                                          >  2:787568/1‑100 (MQ=255)
                                                    tGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTcc                                         >  2:786812/1‑100 (MQ=255)
                                                     gAAAAATATCCATTTGTTTGCAAAAAAAGAGTACTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATCCTTTATTCCa                                        >  2:858445/1‑100 (MQ=255)
                                                       aaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCaga                                      >  2:720640/1‑100 (MQ=255)
                                                       aaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTGATTCCaga                                      >  1:567840/1‑100 (MQ=255)
                                                         aaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCagaag                                    <  1:20195/100‑1 (MQ=255)
                                                           tatCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAagag                                  >  2:973414/1‑100 (MQ=255)
                                                            atacATTTGTTTGAAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAagagn                                 <  2:95991/97‑2 (MQ=255)
                                                            agccAGTTGTTTGCAAAAAAAGGTTATTAGGAAGGGAATTAATGCAATTATCGAAAATTGAAAAAATATCCAAAAATAGTATACTTTATTTCAGAAGAGt                                 >  2:1089260/3‑100 (MQ=255)
                                                                 tttgtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAat                            <  1:331379/100‑1 (MQ=255)
                                                                   tgtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAatat                          >  1:606458/1‑100 (MQ=255)
                                                                        aaaaaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAACATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                     <  2:695800/98‑1 (MQ=255)
                                                                         cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                    >  2:331013/1‑100 (MQ=255)
                                                                         cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                    >  1:492285/1‑100 (MQ=255)
                                                                         cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                    <  1:554991/100‑1 (MQ=255)
                                                                          aaaaaaaGATTATTAGGAAGGGAATTAATGCAATTAAAAAAAAAAAAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                     <  2:265167/98‑1 (MQ=255)
                                                                           aaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTg                  <  1:1063948/100‑1 (MQ=255)
                                                                                aGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTc             >  1:929412/1‑100 (MQ=255)
                                                                                aGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTc             <  1:1022097/100‑1 (MQ=255)
                                                                                     attaGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttt        >  2:80555/1‑100 (MQ=255)
                                                                                       taGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttttt      >  1:841761/1‑100 (MQ=255)
                                                                                        aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc     <  1:1145694/100‑1 (MQ=255)
                                                                                        aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc     >  2:583431/1‑100 (MQ=255)
                                                                                        aggaaggaAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc     <  1:1141753/100‑1 (MQ=255)
                                                                                        aggaaggaAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc     >  2:1141743/1‑100 (MQ=255)
                                                                                          gaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCtt   <  1:732757/100‑1 (MQ=255)
                                                                                          gaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCtt   <  2:481846/100‑1 (MQ=255)
                                                                                          gaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCtt   >  1:41292/1‑100 (MQ=255)
                                                                                          gaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCtt   >  1:1196838/1‑100 (MQ=255)
                                                                                           aagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTa  <  2:69284/100‑1 (MQ=255)
                                                                                               |                                                                                               
TGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTA  >  NC_000913/1197644‑1197834

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