Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,207,981 G→A A161T (GCC→ACC)  stfP → e14 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,207,9810GA95.2% 188.0 / 0.8 64A161T (GCC→ACC) stfPe14 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (26/34);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01

AATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGG  >  NC_000913/1207886‑1208076
                                                                                               |                                                                                               
aaTGAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCtg                                                                                             <  2:93694/100‑1 (MQ=255)
    aGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGtt                                                                                         <  1:1207227/100‑1 (MQ=255)
     gTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGCATTCCCGGTTGCGTTTAATAATGCCTGTGTTg                                                                                        >  2:1134434/1‑100 (MQ=255)
      tAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGc                                                                                       <  2:1509192/100‑1 (MQ=255)
         tCATTCAGTGGGGGTGGAATGCCGCAAAAGCTTCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCgt                                                                                    <  2:621296/100‑1 (MQ=255)
           aTTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCgttg                                                                                  <  2:1330302/100‑1 (MQ=255)
             tCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCgttgtt                                                                                >  2:466065/1‑100 (MQ=255)
              cAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGttgttg                                                                               <  1:233461/100‑1 (MQ=255)
              cAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGCATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGttgttg                                                                               <  1:1134439/100‑1 (MQ=255)
               agggggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGc                                                                              <  2:488007/97‑1 (MQ=255)
               aGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGc                                                                              <  2:1154118/100‑1 (MQ=255)
               aGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGc                                                                              >  1:367805/1‑100 (MQ=255)
                  gggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCaaa                                                                           <  1:1183429/100‑1 (MQ=255)
                   ggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAAt                                                                          >  1:999762/1‑100 (MQ=255)
                   ggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAAt                                                                          <  1:327245/100‑1 (MQ=255)
                    gggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAAtt                                                                         <  2:1220340/100‑1 (MQ=255)
                     ggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTa                                                                        <  1:1033697/100‑1 (MQ=255)
                      gTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATTGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATTCCTGTGTTGCCGTTGTTGCAAATTAt                                                                       >  1:117311/1‑100 (MQ=255)
                      gTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTAt                                                                       >  1:1321666/1‑100 (MQ=255)
                      gTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTAt                                                                       <  2:1042379/100‑1 (MQ=255)
                      gTGGAATGCCGCAAAAGCAGCTGCCTATGGGGGGGAGGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTAt                                                                       <  1:1109576/100‑1 (MQ=255)
                       tGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATg                                                                      <  2:134481/100‑1 (MQ=255)
                          aaTGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACa                                                                   >  2:468802/1‑100 (MQ=255)
                              ccGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGACGTTGTTGAAAATTATGACaa                                                                  >  2:649919/1‑97 (MQ=255)
                                gCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCa                                                             >  1:619359/1‑100 (MQ=255)
                                 cAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTACCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCa                                                             >  2:610623/1‑99 (MQ=255)
                                  aaaaGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgc                                                           <  2:467618/100‑1 (MQ=255)
                                   aaaGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgcg                                                          <  1:1568841/100‑1 (MQ=255)
                                         tCTGCCTCTGGGGGGGATGGAAATACAGTTGTCTTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGAACCTa                                                    >  2:1200966/1‑100 (MQ=255)
                                            gCCTCTGGGGGGGATGGAAATACCGTTGTATTCCCGGTTGCGTGTAATAATACCTGAGTTGCTGTTTTTGCAATTTATGACAAGCTCAgcg                                                          >  2:1569964/1‑91 (MQ=255)
                                             ccTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCaa                                                >  2:184838/1‑100 (MQ=255)
                                              ctctGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAAt                                               <  2:999755/100‑1 (MQ=255)
                                               tctGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATg                                              >  2:445870/1‑100 (MQ=255)
                                               tctGGGGGAGATGGAAATACAGTTGTATTCCTGGTTGCGTTTAAGAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTTAGCGCACCTATCAATg                                              <  2:540874/100‑1 (MQ=255)
                                                  gggggggATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAg                                           >  1:64352/1‑100 (MQ=255)
                                                    gggggATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTg                                         >  1:327113/1‑100 (MQ=255)
                                                        gATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCaa                                     <  1:466066/100‑1 (MQ=255)
                                                              aaTACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACggggg                               <  1:1318274/100‑1 (MQ=255)
                                                               aTACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACgggggg                              <  2:1482569/100‑1 (MQ=255)
                                                               aTACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAAGGTCAGCGCACCTATCAATGCAGTGGCAACgggggg                              <  2:1380301/100‑1 (MQ=255)
                                                                  cAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGata                           >  1:791039/1‑100 (MQ=255)
                                                                  cAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGata                           <  1:845112/100‑1 (MQ=255)
                                                                    nttGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGAtata                         >  2:453689/2‑100 (MQ=255)
                                                                    gTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCACATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGAtata                         >  2:1560826/1‑100 (MQ=255)
                                                                    gTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGAtata                         >  1:652175/1‑100 (MQ=255)
                                                                      tGTATTCGCGGTTGCGTTTAATAATACCTGTTTTGCCGTTGTTGAAAATTATGACAATGTCAGCGCACCTATTTATTTAGTTGCAACGGAGGGATATACa                       <  2:161832/100‑1 (MQ=255)
                                                                               ggTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCg              <  2:867079/100‑1 (MQ=255)
                                                                                gTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGt             >  1:1330430/1‑100 (MQ=255)
                                                                                gTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGt             >  2:1044301/1‑100 (MQ=255)
                                                                                gTTGCGTTTAATAATACCTGTGTTGCCATTGTTGCAAATTATGAGAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGt             <  2:1358128/100‑1 (MQ=255)
                                                                                  tGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGttt           <  2:1069646/100‑1 (MQ=255)
                                                                                   gCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGtttt          >  1:497435/1‑100 (MQ=255)
                                                                                    cGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTCTCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGttttt         >  2:710166/1‑100 (MQ=255)
                                                                                    cGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGttttt         >  2:636764/1‑100 (MQ=255)
                                                                                     gTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTtta        >  1:99966/1‑100 (MQ=255)
                                                                                     gTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTtta        >  1:1426190/1‑100 (MQ=255)
                                                                                      tttAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTttat       <  2:478952/100‑1 (MQ=255)
                                                                                      tttAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTttat       <  1:1276606/100‑1 (MQ=255)
                                                                                      tttAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTttat       <  2:1321659/100‑1 (MQ=255)
                                                                                        taataaTGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTttatta     <  2:682977/100‑1 (MQ=255)
                                                                                        taataaTACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTttatta     >  1:912030/1‑100 (MQ=255)
                                                                                           taatACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACgg  <  1:1207030/100‑1 (MQ=255)
                                                                                           taatACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACgg  <  1:636766/100‑1 (MQ=255)
                                                                                           taatACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACgg  <  1:1002245/100‑1 (MQ=255)
                                                                                               |                                                                                               
AATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGG  >  NC_000913/1207886‑1208076

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: