Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,763,869 (T)7→8 coding (1285/1617 nt) abpA ← CP4‑57 prophage; anti‑bacteriophage protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,763,8621.T100.0% 64.0 / NA 19I431N (ATC→AAC) abpACP4‑57 prophage; anti‑bacteriophage protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (9/10);  total (9/10)

CGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCT  >  NC_000913/2763788‑2763932
                                                                           |                                                                      
cGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGAtttttt                                                                  >  1:445153/1‑81 (MQ=255)
 gAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGAttttttt                                                                 >  1:445160/1‑81 (MQ=255)
     ggTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGAttttttttgcg                                                             <  1:445111/81‑5 (MQ=255)
           cTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGt                                                        <  1:445112/81‑1 (MQ=255)
           cTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGt                                                        <  2:445113/81‑1 (MQ=255)
               cATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCg                                                    <  1:445114/81‑1 (MQ=255)
                       aTTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATc                                            <  2:445115/81‑1 (MQ=255)
                         ttCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATcc                                           >  2:445161/1‑80 (MQ=255)
                          tCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCg                                          <  2:445116/80‑1 (MQ=255)
                                 tGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTcc                                  <  2:445117/81‑1 (MQ=255)
                                     gTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCg                               >  2:445159/1‑80 (MQ=255)
                                     gTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCg                               >  2:445166/1‑80 (MQ=255)
                                          gAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAAc                         >  1:445155/1‑81 (MQ=255)
                                          gAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAAc                         >  1:445167/1‑81 (MQ=255)
                                          gAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAAc                         >  2:445162/1‑81 (MQ=255)
                                              cAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGAt                     <  2:445118/81‑1 (MQ=255)
                                                  tAGTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCg                 >  1:445172/1‑81 (MQ=255)
                                                    gTTATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTa               <  1:445119/81‑1 (MQ=255)
                                                     ttATCGTCTTTAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTAt              <  2:445120/81‑1 (MQ=255)
                                                               tAATGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAAt    >  2:445156/1‑81 (MQ=255)
                                                                 aTGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATc   <  1:445121/80‑1 (MQ=255)
                                                                 aTGTCGAGGATTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCt  >  2:445154/1‑81 (MQ=255)
                                                                           |                                                                      
CGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCT  >  NC_000913/2763788‑2763932

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