Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,419,344 G→A A38V (GCT→GTT)  ydaG ← Rac prophage, uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,419,3440GA96.5% 173.6 / ‑1.7 58A38V (GCT→GTT) ydaGRac prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (26/30);  total (27/31)
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.30e-01

TTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCG  >  NC_000913/1419214‑1419478
                                                                                                                                  |                                                                                                                                      
ttGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTct                                                                                                                                <  1:1366100/139‑1 (MQ=255)
 tGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctg                                                                                                                               <  1:7519/139‑1 (MQ=255)
     tttGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctgcctg                                                                                                                           <  2:213683/139‑1 (MQ=255)
              cGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATc                                                                                                                  <  2:781059/139‑1 (MQ=255)
              cGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATc                                                                                                                  >  1:333274/1‑139 (MQ=255)
                   tGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGt                                                                                                             >  1:1408817/1‑139 (MQ=255)
                    ggCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctg                                                                                                                               <  1:1036436/120‑1 (MQ=255)
                    ggCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctg                                                                                                                               >  2:1036436/1‑120 (MQ=255)
                    ggCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTc                                                                                                            >  2:1193676/1‑139 (MQ=255)
                           aaCGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgctg                                                                                                     >  1:572606/1‑139 (MQ=255)
                              gTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGt                                                                                                  <  1:1112339/139‑1 (MQ=255)
                              gTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGt                                                                                                  <  2:1408817/139‑1 (MQ=255)
                                tCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTaa                                                                                                >  2:836838/1‑139 (MQ=255)
                                      cTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTg                                                                                          <  1:836838/139‑1 (MQ=255)
                                          ttGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGcaca                                                                                      >  2:59087/1‑139 (MQ=255)
                                          ttGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGcaca                                                                                      >  1:1430260/1‑139 (MQ=255)
                                            ggggAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGc                                                                                    <  1:512646/139‑1 (MQ=255)
                                            ggggAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACAGCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGc                                                                                    <  1:1181300/139‑1 (MQ=255)
                                              ggAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCa                                                                                                                        <  1:1185444/101‑1 (MQ=255)
                                              ggAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCa                                                                                                                        >  2:1185444/1‑101 (MQ=255)
                                                  cacGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGtt                                                                              <  2:986744/139‑1 (MQ=255)
                                                    cgccgcACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGt                                                                                                  <  1:312189/117‑1 (MQ=255)
                                                    cgccgcACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGt                                                                                                  >  2:312189/1‑117 (MQ=255)
                                                        gcACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCtgcctgcc                                                                                                                         <  1:233325/90‑1 (MQ=255)
                                                        gcACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCtgcctgcc                                                                                                                         >  2:233325/1‑90 (MQ=255)
                                                        gcACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCg                                                                                                           >  2:496454/1‑104 (MQ=255)
                                                        gcACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCg                                                                                                           <  1:496454/104‑1 (MQ=255)
                                                                 ggATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAg                                                                                              >  1:1615154/1‑108 (MQ=255)
                                                                 ggATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAg                                                                                              <  2:1615154/108‑1 (MQ=255)
                                                                      cGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAAt                                                          >  2:1385465/1‑139 (MQ=255)
                                                                      cGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAAt                                                          >  1:1596702/1‑139 (MQ=255)
                                                                      cGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAAt                                                          >  2:1276627/1‑139 (MQ=255)
                                                                      cGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAAt                                                          <  1:1276627/139‑1 (MQ=255)
                                                                      agttgtaGCTGAGGTCGATTTTTTGAATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCCGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAAt                                                          <  2:333274/138‑1 (MQ=255)
                                                                        ttgttgCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATcc                                                        <  1:59087/139‑1 (MQ=255)
                                                                          gttgCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAg                                                                                              <  2:978069/99‑1 (MQ=255)
                                                                          gttgCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAg                                                                                              >  1:978069/1‑99 (MQ=255)
                                                                           ttgCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATc                                                     >  1:704802/1‑139 (MQ=255)
                                                                             gCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTgg                                                                                >  1:817236/1‑110 (MQ=255)
                                                                             gCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTgg                                                                                <  2:817236/110‑1 (MQ=255)
                                                                                     cGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGaa                                           <  1:1385465/139‑1 (MQ=255)
                                                                                      gATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAac                                          <  2:704802/139‑1 (MQ=255)
                                                                                              gCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGt                                  <  2:1430260/139‑1 (MQ=255)
                                                                                                  tttGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATg                              >  2:819880/1‑139 (MQ=255)
                                                                                                   ttGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGc                                                                                    <  2:1026173/84‑1 (MQ=255)
                                                                                                   ttGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGc                                                                                    >  1:1026173/1‑84 (MQ=255)
                                                                                                   ttGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGc                             <  2:572606/139‑1 (MQ=255)
                                                                                                     gCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCAc                           <  1:1193676/139‑1 (MQ=255)
                                                                                                          tctcACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCt                      <  2:1596702/139‑1 (MQ=255)
                                                                                                              acaTCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAacta                                        <  1:999154/117‑1 (MQ=255)
                                                                                                              acaTCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAacta                                        >  2:999154/1‑117 (MQ=255)
                                                                                                              acaTCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCAc                           <  2:727330/130‑1 (MQ=255)
                                                                                                              acaTCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCAc                           >  1:727330/1‑130 (MQ=255)
                                                                                                                aTCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGca                                                                                        <  2:80052/67‑1 (MQ=255)
                                                                                                                aTCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGca                                                                                        >  1:80052/1‑67 (MQ=255)
                                                                                                                             gTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGcc   >  2:1200081/1‑139 (MQ=255)
                                                                                                                              tataACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCg  >  2:512287/1‑139 (MQ=255)
                                                                                                                              tataACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCg  >  1:1463662/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                      
TTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCG  >  NC_000913/1419214‑1419478

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