Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,003,480 T→C F60S (TTC→TCC)  yedK → SOS response‑associated peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,003,4800TC92.0% 60.5 / ‑0.7 25F60S (TTC→TCC) yedKSOS response‑associated peptidase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (11/12);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

CCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATG  >  NZ_CP009273/2003396‑2003559
                                                                                    |                                                                               
ccATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTggg                                                                            <  1:467677/90‑1 (MQ=255)
 cATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTgggg                                                                           >  1:229635/1‑90 (MQ=255)
           tACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTccc                                                                 >  2:83717/1‑90 (MQ=255)
               aCGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAt                                                             <  1:331048/90‑1 (MQ=255)
                 gTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGAtgg                                                           <  1:362445/90‑1 (MQ=255)
                 gTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTTCCGGAtgg                                                           <  1:528143/90‑1 (MQ=255)
                  tCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggt                                                          <  1:411019/90‑1 (MQ=255)
                    gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGAtggtgg                                                        >  1:102196/1‑90 (MQ=255)
                    gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggtgg                                                        >  1:342490/1‑90 (MQ=255)
                                 aaGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACcgccgc                                           <  1:503871/90‑1 (MQ=255)
                                             tCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGGGGGATAAACCGCCGCTGATTAACGccc                               >  1:401706/1‑90 (MQ=255)
                                             ccAGTGAACGTGATGAACCCCTTCATCTGGATCCGGTTTCCTGGGGATATGCCCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                               <  2:681469/89‑1 (MQ=255)
                                                gTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCcgcg                            >  2:772367/1‑90 (MQ=255)
                                                gTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCcgcg                            >  2:39489/1‑90 (MQ=255)
                                                     cGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGaa                       >  1:652822/1‑90 (MQ=255)
                                                     cGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGaa                       <  1:39140/90‑1 (MQ=255)
                                                     cGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGaa                       >  1:350685/1‑90 (MQ=255)
                                                         atgaACCCCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTg                   <  2:187050/90‑1 (MQ=255)
                                                          tgaACCCCTTCATCTGGATCCGGTTTCCGGGGGATATGCCCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGc                  <  2:692081/90‑1 (MQ=255)
                                                          tgaACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGc                  >  2:605650/1‑90 (MQ=255)
                                                            aaCACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCgg                <  1:667597/90‑1 (MQ=255)
                                                                 cTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAAcc           <  1:605650/90‑1 (MQ=255)
                                                                       cTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGt     <  1:39489/90‑1 (MQ=255)
                                                                          gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATg  >  2:11929/1‑90 (MQ=255)
                                                                          gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATg  >  2:514693/1‑90 (MQ=255)
                                                                                    |                                                                               
CCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATG  >  NZ_CP009273/2003396‑2003559

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