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,4800TC96.0% 72.9 / NA 25F60S (TTC→TCC) yedKSOS response‑associated peptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (8/16);  minor base A (1/0);  total (9/16)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAAC  >  NZ_CP009273/2003393‑2003566
                                                                                       |                                                                                      
aaCCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCt                                                                                      <  2:133799/90‑1 (MQ=255)
    cATTGGCAGATACAACGTCGCGCCCGCGACTAAAGTCCTACTGCTCAATACACGTGATGAACACCTTCCTCTGACTCCGGTTTACTgctg                                                                                  >  2:335011/1‑87 (MQ=255)
      ttGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGat                                                                                >  1:468760/1‑90 (MQ=255)
       tGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGata                                                                               <  2:630281/90‑1 (MQ=255)
         gCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATg                                                                             <  2:165714/90‑1 (MQ=255)
              tACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTccc                                                                        <  2:115519/90‑1 (MQ=255)
                    gTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtgg                                                                  <  2:272461/90‑1 (MQ=255)
                       gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGatat                                                                              >  2:75781/1‑75 (MQ=255)
                       gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGatat                                                                              <  1:75781/75‑1 (MQ=255)
                       gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggtgg                                                               >  1:250909/1‑90 (MQ=255)
                                    aaGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACcgccgc                                                  <  2:140275/90‑1 (MQ=255)
                                             tGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGACAAACCGCCGCTGATTAACg                                         >  2:61440/1‑90 (MQ=255)
                                              gCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGc                                        >  1:568653/1‑90 (MQ=255)
                                                tCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGatat                                                                              <  1:543693/50‑1 (MQ=255)
                                                tCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGatat                                                                              >  2:543693/1‑50 (MQ=255)
                                                tCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                                      <  1:82083/90‑1 (MQ=255)
                                                                  acCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCaa                    >  2:365258/1‑90 (MQ=255)
                                                                    cTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAGACTGCGGCAAcc                  <  1:278566/90‑1 (MQ=255)
                                                                    cTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAAcc                  <  1:558745/90‑1 (MQ=255)
                                                                    cTTCATCTGGATCCGGTTTCCTGGGGATATGCCCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAAcc                  <  2:477458/90‑1 (MQ=255)
                                                                          cTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGt            <  1:515513/90‑1 (MQ=255)
                                                                           tGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTa           <  1:577745/90‑1 (MQ=255)
                                                                           tGGATCCGGTTTCCTGGGGATATCCCCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTa           <  2:609560/90‑1 (MQ=255)
                                                                             gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGAt                                                            <  2:45945/39‑1 (MQ=255)
                                                                             gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGAt                                                            >  1:45945/1‑39 (MQ=39)
                                                                                    tttcCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAAc  <  2:568653/90‑1 (MQ=255)
                                                                                       |                                                                                      
AACCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAAC  >  NZ_CP009273/2003393‑2003566

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