Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 294,762 A→G V13A (GTT→GCT)  yagM ← CP4‑6 prophage; protein YagM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913294,7620AG100.0% 152.7 / NA 42V13A (GTT→GCT) yagMCP4‑6 prophage; protein YagM
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/26);  total (16/26)

ATAATGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTA  >  NC_000913/294638‑294888
                                                                                                                            |                                                                                                                              
aTAATGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTc                                                                                                                        <  2:410049/133‑1 (MQ=255)
             aCTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAg                                                                                                           <  2:473749/133‑1 (MQ=255)
              cTCCTAATCAAGATTCATCTTTCAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGt                                                                                                          <  1:172247/133‑1 (MQ=255)
                ccTAATGACGATTGACTTATGAGCGTGACTCAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTaa                                                                                                        <  1:295634/133‑1 (MQ=255)
                ccTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTaa                                                                                                        <  2:398619/133‑1 (MQ=255)
                ccTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTaa                                                                                                        <  2:352284/133‑1 (MQ=255)
                ccTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTaa                                                                                                        <  1:371367/133‑1 (MQ=255)
                  tAATGAAGATTGATTTTTGAACGTGACTTAGTCTTTTTTTATGGGAGTAGTCTGTACTGCCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAaca                                                                                                      <  2:207092/133‑1 (MQ=255)
                   aaTGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAacac                                                                                                     >  2:42834/1‑133 (MQ=255)
                     tGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTa                                                                                                   >  1:557187/1‑133 (MQ=255)
                     tGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTATTCCTGGTGGAGAGCTTCTCATCTCTAAATGAGTAACACTa                                                                                                   >  2:164195/1‑133 (MQ=255)
                         gattgattTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGGTGCTCATGTCACAAGTAGTAAGACTATTAg                                                                                               >  1:199615/1‑133 (MQ=255)
                                      aCGTGACTGATTGTTTTTTTATGGCAGTTGTCCGTACTACCAGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGa                                                                                  <  2:199615/133‑1 (MQ=255)
                                      aCGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGa                                                                                  <  2:105860/133‑1 (MQ=255)
                                      aCGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGa                                                                                  <  1:71515/133‑1 (MQ=255)
                                      aCGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGa                                                                                  <  1:83304/133‑1 (MQ=255)
                                      aCGTGACTCAGTGTTTTTTTATTGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGa                                                                                  <  1:209823/133‑1 (MQ=255)
                                            cTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACctc                                                                            <  2:426408/133‑1 (MQ=255)
                                                      ttttATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTAtct                                                                  <  1:42834/133‑1 (MQ=255)
                                                           tGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAata                                                             <  1:333692/133‑1 (MQ=255)
                                                                          ctacCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTg                                              <  1:594343/133‑1 (MQ=255)
                                                                           tacCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGa                                             >  2:385515/1‑133 (MQ=255)
                                                                             ccTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGaaa                                           >  2:595076/1‑133 (MQ=255)
                                                                             ccTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGaaa                                           >  1:592350/1‑133 (MQ=255)
                                                                             ccTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGaaa                                           >  1:40059/1‑133 (MQ=255)
                                                                             ccTGTGACCCATTGTATTCGCCCTGACAAGCTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGaaa                                           <  1:164195/133‑1 (MQ=255)
                                                                                  gACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTtct                                      >  2:317606/1‑133 (MQ=255)
                                                                                  gACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTtct                                      <  1:348756/133‑1 (MQ=255)
                                                                                  gACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTtct                                      >  1:24164/1‑133 (MQ=255)
                                                                                  gACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTtct                                      >  1:585696/1‑133 (MQ=255)
                                                                                       aTTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAgtg                                 >  1:422436/1‑133 (MQ=255)
                                                                                               tGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCa                         <  2:316928/133‑1 (MQ=255)
                                                                                               tGCCTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCa                         >  2:436726/1‑133 (MQ=255)
                                                                                                 ccTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCAtt                       >  1:556289/1‑133 (MQ=255)
                                                                                                 ccTTGACAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACCATCTGAAATTTCACAGTGATGCCCCAtt                       >  2:279464/1‑133 (MQ=255)
                                                                                                      aCAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTtct                  <  1:436726/133‑1 (MQ=255)
                                                                                                      aCAAGGTTTTTTTCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTtct                  <  2:585696/133‑1 (MQ=255)
                                                                                                            tttttttCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGc            <  1:448791/133‑1 (MQ=255)
                                                                                                            tttttttCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGc            <  2:40059/133‑1 (MQ=255)
                                                                                                               ttttCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAAc         <  2:556289/133‑1 (MQ=255)
                                                                                                                 ttCCTGGTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACAt       <  1:23532/133‑1 (MQ=255)
                                                                                                                      ggTAGAGCGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCAATTTACTGGGGGCAGCATCGCta  >  1:166320/1‑133 (MQ=255)
                                                                                                                            |                                                                                                                              
ATAATGATCAATTACTCCTAATGAAGATTGATTTTTGAACGTGACTGAGTGTTTTTTTATGGGAGTAGTCTGTACTACCTGTGACTCATTGTTTTTGCCTTGACAAGGTTTTTTTCCTGGTAGAACGCTGCTCATCTCAAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTA  >  NC_000913/294638‑294888

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