Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 281,662 Δ1 bp coding (322/1155 nt) yagA ← CP4‑6 prophage, putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913281,6600G.98.0% 194.0 / ‑2.7 49coding (324/1155 nt)yagACP4‑6 prophage, putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (21/27);  minor base C (0/1);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGA  >  NC_000913/281525‑281779
                                                                                                                                       |                                                                                                                       
aaaGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggca                                                                                                                      <  1:1464319/139‑3 (MQ=255)
      gCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGgtgt                                                                                                               >  2:360348/1‑139 (MQ=255)
      gCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGgtgt                                                                                                               >  1:1395774/1‑139 (MQ=255)
         cTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGccc                                                                                                           >  2:829891/1‑139 (MQ=255)
           tAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGcccc                                                                                                          <  2:715110/139‑1 (MQ=255)
                         cAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCa                                                                                            <  1:550226/139‑1 (MQ=255)
                           gcggCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCCGGGCATCACGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAgc                                                                                          >  2:1474552/3‑139 (MQ=255)
                             ggCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAgcgc                                                                                        <  2:1158697/139‑1 (MQ=255)
                              gCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCt                                                                                       <  2:199086/139‑1 (MQ=255)
                                                 tCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCa                                                                    <  2:1000649/139‑1 (MQ=255)
                                                          cccGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATg                                                           >  1:1244526/1‑139 (MQ=255)
                                                             gTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACg                                                        <  1:360348/139‑1 (MQ=255)
                                                                      aTGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCggg                                                                           <  1:459169/111‑1 (MQ=255)
                                                                      aTGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCggg                                                                           >  2:459169/1‑111 (MQ=255)
                                                                          ccGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATa                                           <  2:820475/139‑1 (MQ=255)
                                                                                   cgcCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCTCAGCGTTCATGACGGTCATGGGCCATACg                                         <  2:871856/132‑1 (MQ=255)
                                                                                   cgcCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACg                                         >  1:1089704/1‑132 (MQ=255)
                                                                                   cgcCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACg                                         <  2:1089704/132‑1 (MQ=255)
                                                                                   cgcCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACg                                         >  1:871856/1‑132 (MQ=255)
                                                                                        ggCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa                                                          <  1:821637/110‑1 (MQ=255)
                                                                                        ggCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa                                                          >  2:821637/1‑110 (MQ=255)
                                                                                         gcagcaGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTg                            <  1:829891/139‑1 (MQ=255)
                                                                                         gcagcaGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTg                            <  2:915950/139‑1 (MQ=255)
                                                                                               ggCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCtt                                                                                >  1:1494913/1‑81 (MQ=255)
                                                                                               ggCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCtt                                                                                <  2:1494913/81‑1 (MQ=255)
                                                                                                    tGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATggg                                                >  1:1241969/1‑108 (MQ=255)
                                                                                                    tGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATggg                                                <  2:1241969/108‑1 (MQ=255)
                                                                                                        gggCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg             <  2:758504/139‑1 (MQ=255)
                                                                                                          gCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCGTACgcagca                                    <  1:1497177/114‑1 (MQ=255)
                                                                                                          gCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCGTACgcagca                                    >  2:1497177/1‑114 (MQ=255)
                                                                                                          gCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGc            >  2:1082070/1‑138 (MQ=255)
                                                                                                          gCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGc            <  1:1082070/138‑1 (MQ=255)
                                                                                                                   ttATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa  <  1:1474552/139‑1 (MQ=255)
                                                                                                                    tATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTcccc                                                                     >  2:58930/1‑71 (MQ=255)
                                                                                                                    tATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTcccc                                                                     <  1:58930/71‑1 (MQ=255)
                                                                                                                       ggACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGcc                               >  2:894413/1‑106 (MQ=255)
                                                                                                                       ggACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGcc                               <  1:894413/106‑1 (MQ=255)
                                                                                                                            gTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCg      <  2:179035/126‑1 (MQ=255)
                                                                                                                            gTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCg      >  1:179035/1‑126 (MQ=255)
                                                                                                                            gTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa  >  2:1206972/1‑130 (MQ=255)
                                                                                                                            gTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa  <  1:1206972/130‑1 (MQ=255)
                                                                                                                             tGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATggg                                                >  1:859363/1‑83 (MQ=255)
                                                                                                                             tGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATggg                                                <  2:859363/83‑1 (MQ=255)
                                                                                                                             tGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATgg                                                 >  2:923992/1‑82 (MQ=255)
                                                                                                                             tGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg             >  1:59511/1‑118 (MQ=255)
                                                                                                                             tGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg             <  1:1276991/118‑1 (MQ=255)
                                                                                                                             tGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg             <  2:59511/118‑1 (MQ=255)
                                                                                                                             tGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg             >  2:1276991/1‑118 (MQ=255)
                                                                                                                             tGCTGAAGCCGGCATGGTGTGCCCCTGGTCTTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATgg                                                 <  1:923992/82‑1 (MQ=255)
                                                                                                                                tGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa  <  2:1194712/126‑1 (MQ=255)
                                                                                                                                tGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa  >  1:1194712/1‑126 (MQ=255)
                                                                                                                                       |                                                                                                                       
AAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGA  >  NC_000913/281525‑281779

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: