Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,432,273 G→A A1075A (GCG→GCA stfR → Rac prophage, putative tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,432,2730GA96.9% 196.6 / ‑0.9 64coding (3225/3363 nt)stfRRac prophage, putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/37);  minor base . (1/1);  total (26/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CGGCTAGTGCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGTGAGGC  >  NC_000913/1432137‑1432403
                                                                                                                                        |                                                                                                                                  
cGGCTAGTGCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGcaca                                                                                                                                  <  2:390515/139‑1 (MQ=255)
cGGCTAGTGCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGcaca                                                                                                                                  <  1:597836/139‑1 (MQ=255)
  gCTAGTGCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATa                                                                                                                                <  2:95896/139‑1 (MQ=255)
  gCTAGTGCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATa                                                                                                                                <  1:521659/139‑1 (MQ=255)
   cTAGTGCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATAc                                                                                                                               >  2:857596/1‑139 (MQ=255)
               ccGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTgg                                                                                                                   <  2:702991/139‑1 (MQ=255)
                 ggTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTg                                                                                                                 >  1:959384/1‑139 (MQ=255)
                  gTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGc                                                                                                                >  1:1033860/1‑139 (MQ=255)
                   tGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCGACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACGGTCGGTATTGGTGCt                                                                                                               <  1:1075578/139‑1 (MQ=255)
                    gCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTc                                                                                                              <  2:303300/139‑1 (MQ=255)
                      tGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCAt                                                                                                            <  2:1033001/139‑1 (MQ=255)
                          aGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGc                                                                                                        <  2:1581004/139‑1 (MQ=255)
                           gTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCa                                                                                                       <  1:90181/139‑1 (MQ=255)
                                  aacaaGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGtt                                                                                                >  1:634680/1‑139 (MQ=255)
                                   acaaGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTg                                                                                               >  2:1263670/1‑139 (MQ=255)
                                          tGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTgg                                                                                        <  1:1263670/139‑1 (MQ=255)
                                            tCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTAt                                                                                                                      <  1:1054206/107‑1 (MQ=255)
                                            tCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCACACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTAt                                                                                                                      >  2:1054206/1‑107 (MQ=255)
                                                 tgtgCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGtt                                                                                      >  1:1435886/1‑134 (MQ=255)
                                                 tgtgCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGtt                                                                                      <  2:1435886/134‑1 (MQ=255)
                                                      aTAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGaca                                                                            >  2:1255503/1‑139 (MQ=255)
                                                        aaTCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGacaca                                                                          <  2:45495/139‑1 (MQ=255)
                                                             aaaCTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATc                                                                     >  1:1508790/1‑139 (MQ=255)
                                                               aCTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCAc                                                                   >  1:1474992/1‑139 (MQ=255)
                                                                  aTGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGt                                                                >  1:1214617/1‑139 (MQ=255)
                                                                   tGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCTGTTGCGATTGGTTCACATGGa                                                                              >  1:1259764/1‑124 (MQ=255)
                                                                   tGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCTGTTGCGATTGGTTCACATGGa                                                                              <  2:1259764/124‑1 (MQ=255)
                                                                    gCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCGGGTGCACACGCCATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGaca                                                                            <  1:694656/124‑1 (MQ=255)
                                                                    gCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCGGGTGCACACGCCATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGaca                                                                            >  2:694656/1‑124 (MQ=255)
                                                                    gCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGAt                                                                                           >  2:449696/1‑110 (MQ=255)
                                                                    gCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGAt                                                                                           <  1:449696/110‑1 (MQ=255)
                                                                           catcTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACgctgc                                                       >  1:825417/1‑139 (MQ=255)
                                                                                gctgGCGCATATACCCACTCACTGTCCGGCACTGTTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTa                                                  <  2:1033860/139‑1 (MQ=255)
                                                                                gctgGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTa                                                  <  2:1214617/139‑1 (MQ=255)
                                                                                  tgGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAAc                                                >  2:1430178/1‑139 (MQ=255)
                                                                                   ggCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAAcg                                               <  1:857596/139‑1 (MQ=255)
                                                                                    gcgcACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGaca                                                                            >  2:782826/1‑109 (MQ=255)
                                                                                    gcgcACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGaca                                                                            <  1:782826/109‑1 (MQ=255)
                                                                                     cgcACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCt                                                                                                               >  2:693818/1‑73 (MQ=255)
                                                                                     cgcACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCt                                                                                                               <  1:693818/73‑1 (MQ=255)
                                                                                            aCCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTTTCGGTATTGGTGCTCATACGCACTCCGTTTCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAaca                                      <  2:1100466/139‑1 (MQ=255)
                                                                                               cactcactGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCg                                   >  2:1265473/1‑139 (MQ=255)
                                                                                                    actGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCa                                                                       >  2:1068387/1‑98 (MQ=255)
                                                                                                    actGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCa                                                                       <  1:1068387/98‑1 (MQ=255)
                                                                                                      tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGAt                                                                                           >  1:652584/1‑76 (MQ=255)
                                                                                                      tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGAt                                                                                           <  2:652584/76‑1 (MQ=255)
                                                                                                      tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGa                                                                              <  2:1308748/89‑1 (MQ=255)
                                                                                                      tGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGa                                                                              >  1:1308748/1‑89 (MQ=255)
                                                                                                        tCCGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACa                          >  1:846248/1‑139 (MQ=37)
                                                                                                         ccGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACAt                         <  1:1430178/139‑1 (MQ=21)
                                                                                                         ccGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACAt                         <  1:247832/139‑1 (MQ=21)
                                                                                                          cGGCACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATc                        <  2:634680/139‑1 (MQ=21)
                                                                                                             cACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                     <  2:1508790/139‑1 (MQ=37)
                                                                                                             cACTGCTGCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCa                     <  2:825417/139‑1 (MQ=37)
                                                                                                               ctgctgCAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCAtt                   <  1:1255503/139‑1 (MQ=37)
                                                                                                                  ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTca                                                                                    >  2:216422/1‑71 (MQ=255)
                                                                                                                  ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTca                                                                                    <  1:216422/71‑1 (MQ=255)
                                                                                                                  ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACAcc                                                                        <  1:1352973/83‑1 (MQ=255)
                                                                                                                  ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACAcc                                                                        >  2:1352973/1‑83 (MQ=255)
                                                                                                                  ctgcAAGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTaa                <  2:959384/139‑1 (MQ=21)
                                                                                                                      aaGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACtat            <  1:495885/139‑1 (MQ=21)
                                                                                                                       aGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATGCGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAAc                                                >  2:1052108/1‑102 (MQ=255)
                                                                                                                       aGCGCAGGTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGATGGTAAc                                                <  1:1052108/102‑1 (MQ=255)
                                                                                                                              gTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGTGAg    <  2:768219/139‑1 (MQ=18)
                                                                                                                              gTGCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGTGAg    <  1:1265473/139‑1 (MQ=18)
                                                                                                                                gCACACGCACATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGTGAGGc  <  2:846248/139‑1 (MQ=14)
                                                                                                                                        |                                                                                                                                  
CGGCTAGTGCTAACTCCGGTGCTGGTAGTGCATCAACAAGATTGTCTGTTGTGCATAATCAAAACTATGCAACATCATCTGCTGGCGCACATACCCACTCACTGTCCGGCACTGCTGCAAGCGCAGGTGCACACGCGCATACTGTCGGTATTGGTGCTCATACGCACTCCGTTGCGATTGGTTCACATGGACACACCATCACCGTTAACGCTGCTGGTAACGCGGAAAACACCGTCAAAAACATCGCATTTAACTATATTGTGAGGC  >  NC_000913/1432137‑1432403

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