Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 291,942 A→G K148R (AAG→AGG)  yagJ → CP4‑6 prophage; protein YagJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913291,9420AG95.2% 129.6 / ‑2.2 42K148R (AAG→AGG) yagJCP4‑6 prophage; protein YagJ
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (18/22);  total (19/23)
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.18e-01

CTCTGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAA  >  NC_000913/291807‑292072
                                                                                                                                       |                                                                                                                                  
ctctGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACCGGGAATCACTTATGATGAAGCCGAACAGTACAGGAtt                                                                                                                                >  1:350013/1‑140 (MQ=255)
   tGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGtt                                                                                                                             >  2:706377/1‑140 (MQ=255)
     ttGAACCATTGGGAAACCAAGGATTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTAAGTTCAAACTGAACAGTACAGGTTTGTTCg                                                                                                                           >  1:335796/1‑140 (MQ=255)
          ccATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCg                                                                                                                      <  1:170035/140‑1 (MQ=255)
          ccATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCg                                                                                                                      <  1:706377/140‑1 (MQ=255)
          ccATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCg                                                                                                                      <  1:611126/140‑1 (MQ=255)
          ccATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCCTATGTTGTAGCCGATCAGTACAGGATTGTTCGTCCCg                                                                                                                      >  1:398373/1‑140 (MQ=255)
            aTTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGAt                                                                                                                    <  1:563632/140‑1 (MQ=255)
             ttGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGAtt                                                                                                                   >  2:540894/1‑140 (MQ=255)
             ttGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGAtt                                                                                                                   >  2:549581/1‑140 (MQ=255)
               gggAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTc                                                                                                                 >  1:189743/1‑140 (MQ=255)
               gggAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTc                                                                                                                 >  1:760475/1‑140 (MQ=255)
                     ccAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCGTATGTTGAGGCTGAACCGTACAGGAGTGTTGGTCCCGATTTCCTGTTc                                                                                                           >  2:212329/1‑140 (MQ=255)
                          gTTTTTGAGATCGAAACCAAACGTTAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGc                                                                                                      >  2:648117/1‑140 (MQ=255)
                                           cAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCaaaa                                                                                     <  2:484872/140‑1 (MQ=255)
                                                gTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGtt                                                                                <  2:751895/140‑1 (MQ=255)
                                                gTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGtt                                                                                <  2:561677/140‑1 (MQ=255)
                                                 tGAAGGATCTGCTTTCTTGTATTGTAATCCAAACTCCACAGGGCAGTCATCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTtg                                                                               <  2:398373/140‑1 (MQ=255)
                                                     ggATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGa                                                                           <  2:189743/140‑1 (MQ=255)
                                                      gATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGAc                                                                          <  2:650337/140‑1 (MQ=255)
                                                          tGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTgg                                                                      <  2:762424/140‑1 (MQ=255)
                                                          tGCTTTCTGGTATCGTAATCCAAAGTACACAGGTCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTgg                                                                      <  2:350013/140‑1 (MQ=255)
                                                              ttCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGa                                                                  <  1:692602/140‑1 (MQ=255)
                                                                    tATCGTAATCCACAGTAAACAGTGCAGTCGTCACTGGGAGTCGCTTATATTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCcaca                                                            <  2:335796/140‑1 (MQ=255)
                                                                      tCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATa                                                          >  1:163309/1‑140 (MQ=255)
                                                                              cacaGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAGAATGGTTATGGGCCGGGTAGAGCCACATAGCCTACAt                                                  >  1:235725/1‑140 (MQ=255)
                                                                               acaGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATc                                                 <  1:549581/140‑1 (MQ=255)
                                                                                 aGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAGATGATTGCGGACCTGGGAGATCCACATAGCCTACATCTg                                               >  2:289538/1‑140 (MQ=255)
                                                                                   tACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGc                                             >  1:670399/1‑140 (MQ=255)
                                                                                     cacaGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTg                                           <  1:109026/140‑1 (MQ=255)
                                                                                      acaGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGa                                          >  1:500990/1‑140 (MQ=255)
                                                                                          ggCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCt                                      >  1:767948/1‑140 (MQ=255)
                                                                                             aGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTg                                   >  1:138015/1‑140 (MQ=255)
                                                                                               tcgtcACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCAGGTAGATTGACATAGCATACATTTGGCTGAGGCTTTGcc                                 >  2:394560/1‑140 (MQ=255)
                                                                                                       gggAATCGCTTATGTTGAAGCTCAATAGTGCAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTgg                         <  1:289538/140‑1 (MQ=255)
                                                                                                        ggAATCGCTTATGTTCAAGCTGCACCGTACAGGATTGATCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGa                        <  1:212329/140‑1 (MQ=255)
                                                                                                             cGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGAc                   <  2:760475/140‑1 (MQ=255)
                                                                                                             cGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGAc                   <  2:767948/140‑1 (MQ=255)
                                                                                                                ttATGTTGAAGCTGAACCGTACGGGATTGTTAGTCTCGACTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTg                <  2:235725/140‑1 (MQ=255)
                                                                                                                      tGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCtat          <  2:670399/140‑1 (MQ=255)
                                                                                                                       gAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCtata         >  2:94344/1‑140 (MQ=255)
                                                                                                                           cTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGcc     <  2:500990/140‑1 (MQ=255)
                                                                                                                              aaCAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGaa  >  2:593885/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                  
CTCTGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAA  >  NC_000913/291807‑292072

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