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)
      gCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGgtgt                                                                                                               >  2:360348/1‑139 (MQ=255)
      gCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGgtgt                                                                                                               >  1:1395774/1‑139 (MQ=255)
         cTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCG‑CCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGccc                                                                                                           >  2:829891/1‑139 (MQ=255)
           tAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGcccc                                                                                                          <  2:715110/139‑1 (MQ=255)
                         cAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCa                                                                                            <  1:550226/139‑1 (MQ=255)
                           gcggCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCCGGGCATCACGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAgc                                                                                          >  2:1474552/3‑139 (MQ=255)
                             ggCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAgcgc                                                                                        <  2:1158697/139‑1 (MQ=255)
                              gCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCt                                                                                       <  2:199086/139‑1 (MQ=255)
                                                 tCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCa                                                                    <  2:1000649/139‑1 (MQ=255)
                                                          cccGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATg                                                           >  1:1244526/1‑139 (MQ=255)
                                                             gTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACg                                                        <  1:360348/139‑1 (MQ=255)
                                                                      aTGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCggg                                                                           <  1:459169/111‑1 (MQ=255)
                                                                      aTGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCggg                                                                           >  2:459169/1‑111 (MQ=255)
                                                                          ccGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATa                                           <  2:820475/139‑1 (MQ=255)
                                                                                   cgcCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCTCAGCGTTCATGACGGTCATGGGCCATACg                                         <  2:871856/132‑1 (MQ=255)
                                                                                   cgcCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACg                                         >  1:1089704/1‑132 (MQ=255)
                                                                                   cgcCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACg                                         <  2:1089704/132‑1 (MQ=255)
                                                                                   cgcCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACg                                         >  1:871856/1‑132 (MQ=255)
                                                                                        ggCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa                                                          <  1:821637/110‑1 (MQ=255)
                                                                                        ggCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa                                                          >  2:821637/1‑110 (MQ=255)
                                                                                         gcagcaGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTg                            <  1:829891/139‑1 (MQ=255)
                                                                                         gcagcaGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTg                            <  2:915950/139‑1 (MQ=255)
                                                                                               ggCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCtt                                                                                >  1:1494913/1‑81 (MQ=255)
                                                                                               ggCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCtt                                                                                <  2:1494913/81‑1 (MQ=255)
                                                                                                    tGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATggg                                                >  1:1241969/1‑108 (MQ=255)
                                                                                                    tGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATggg                                                <  2:1241969/108‑1 (MQ=255)
                                                                                                        gggCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg             <  2:758504/139‑1 (MQ=255)
                                                                                                          gCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCGTACgcagca                                    <  1:1497177/114‑1 (MQ=255)
                                                                                                          gCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCGTACgcagca                                    >  2:1497177/1‑114 (MQ=255)
                                                                                                          gCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGc            >  2:1082070/1‑138 (MQ=255)
                                                                                                          gCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGc            <  1:1082070/138‑1 (MQ=255)
                                                                                                                   ttATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa  <  1:1474552/139‑1 (MQ=255)
                                                                                                                    tATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTcccc                                                                     >  2:58930/1‑71 (MQ=255)
                                                                                                                    tATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTcccc                                                                     <  1:58930/71‑1 (MQ=255)
                                                                                                                       ggACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGcc                               >  2:894413/1‑106 (MQ=255)
                                                                                                                       ggACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGcc                               <  1:894413/106‑1 (MQ=255)
                                                                                                                            gTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCg      <  2:179035/126‑1 (MQ=255)
                                                                                                                            gTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCg      >  1:179035/1‑126 (MQ=255)
                                                                                                                            gTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa  >  2:1206972/1‑130 (MQ=255)
                                                                                                                            gTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa  <  1:1206972/130‑1 (MQ=255)
                                                                                                                             tGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATggg                                                >  1:859363/1‑83 (MQ=255)
                                                                                                                             tGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATggg                                                <  2:859363/83‑1 (MQ=255)
                                                                                                                             tGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATgg                                                 >  2:923992/1‑82 (MQ=255)
                                                                                                                             tGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg             >  1:59511/1‑118 (MQ=255)
                                                                                                                             tGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg             <  1:1276991/118‑1 (MQ=255)
                                                                                                                             tGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg             <  2:59511/118‑1 (MQ=255)
                                                                                                                             tGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg             >  2:1276991/1‑118 (MQ=255)
                                                                                                                             tGCTGAA‑GCCGGCATGGTGTGCCCCTGGTCTTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATgg                                                 <  1:923992/82‑1 (MQ=255)
                                                                                                                                tGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa  <  2:1194712/126‑1 (MQ=255)
                                                                                                                                tGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa  >  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: ‑