Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,907,320 A→G F585S (TTC→TCC)  prc ← carboxy terminal‑processing peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,907,3200AG100.0% 92.2 / NA 27F585S (TTC→TCC) prccarboxy terminal‑processing peptidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/11);  total (16/11)

TCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCA  >  NZ_CP009273/1907234‑1907398
                                                                                      |                                                                              
tctcACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGGACt                                                                             >  1:266308/1‑90 (MQ=255)
  tcACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGGACTCa                                                                           >  1:457553/1‑90 (MQ=255)
       cacaGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGGACTCAGGATc                                                                      >  2:709138/1‑90 (MQ=255)
         caGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGGACTCAGGATCtt                                                                    >  1:809048/1‑90 (MQ=255)
           gCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGGACTCAGGATCTTTc                                                                  >  1:1054552/1‑90 (MQ=255)
           gCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGGACTCAGGATCTTTc                                                                  >  2:851784/1‑90 (MQ=255)
           gCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGGACTCAGGATCTTTc                                                                  >  2:339743/1‑90 (MQ=255)
                   cAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGGACTCAGGATCTTTCGCGATAcg                                                          <  1:325529/90‑1 (MQ=255)
                   cAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGGACTCAGGATCTTTCGCGATAcg                                                          >  2:80143/1‑90 (MQ=255)
                     gagaAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGGACTCAGGATCTTTCGCGATAcgcg                                                        <  2:578350/90‑1 (MQ=255)
                     gagaAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGGACTCAGGATCTTTCGCGATAcgcg                                                        >  2:773611/1‑90 (MQ=255)
                     gagaAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGACATCCTTCATGATGTTCTGGGACTCAGGATCTTTCGCGATAcgcg                                                        >  1:358164/1‑90 (MQ=255)
                      agaAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGGACTCAGGATCTTTCGCGATACgcgc                                                       <  2:809048/90‑1 (MQ=255)
                        aaaCGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGGACTCAGGATCTTTCGCGATACGCGCAt                                                     <  1:709138/90‑1 (MQ=255)
                                       tGTCCTTCATAGCGTTGAAGCTCTCGATATCCTTCATGATGTTCTGGGACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTcagca                                      >  1:848529/1‑90 (MQ=255)
                                            ttCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGGACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTcc                                 <  2:790456/90‑1 (MQ=255)
                                            ttCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGGACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTcc                                 <  2:787469/90‑1 (MQ=255)
                                                    gTTGAAGCGCGCGATATCCTTCATGATGTTCTGGGACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCaaa                         <  1:339743/90‑1 (MQ=255)
                                                    gTTGAAGCGCGCGATATCCTTCATGATGTTCTGGGACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCaaa                         <  1:851784/90‑1 (MQ=255)
                                                            gcgcgATATCCTTCATGATGTTCTGGGACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTa                 >  1:103305/1‑90 (MQ=255)
                                                            gcgcgATATCCTTCATGATGTTCTGGGACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTa                 >  2:804543/1‑90 (MQ=255)
                                                                   atCCTTCATGATGTTCTGGGACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTcc          <  1:986942/90‑1 (MQ=255)
                                                                       ttCATGATGTTCTGGGACTCAGGATCTTTCGCGATAcg                                                          >  1:159637/1‑38 (MQ=38)
                                                                       ttCATGATGTTCTGGGACTCAGGATCTTTCGCGATAcg                                                          <  2:159637/38‑1 (MQ=38)
                                                                           tgatgTTCTGGGACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTca  >  1:290929/1‑90 (MQ=255)
                                                                              tgTTCTGGGACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGAtt     <  1:276145/84‑1 (MQ=255)
                                                                              tgTTCTGGGACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGAtt     >  2:276145/1‑84 (MQ=255)
                                                                                      |                                                                              
TCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCA  >  NZ_CP009273/1907234‑1907398

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