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,4800TC93.5% 124.5 / ‑0.5 46F60S (TTC→TCC) yedKSOS response‑associated peptidase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (24/19);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-01

CCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAAC  >  NZ_CP009273/2003395‑2003566
                                                                                     |                                                                                       
cccATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTgg                                                                                    <  2:916731/90‑1 (MQ=255)
   aTTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGa                                                                                 >  2:311500/1‑90 (MQ=255)
   aTTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGa                                                                                 >  2:1009220/1‑90 (MQ=255)
        cAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGc                                                                            <  2:1407156/90‑1 (MQ=255)
          gATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTc                                                                          >  2:548618/1‑90 (MQ=255)
              cAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCgg                                                                      <  2:989618/90‑1 (MQ=255)
                  gTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtgg                                                                  <  2:1255011/90‑1 (MQ=255)
                  gTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtgg                                                                  <  1:1557731/90‑1 (MQ=255)
                   tCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggt                                                                 <  1:480790/90‑1 (MQ=255)
                     gcgtCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggtgg                                                               >  1:1211111/1‑90 (MQ=255)
                     gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGatat                                                                              >  2:1104705/1‑75 (MQ=255)
                     gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGatat                                                                              <  1:1104705/75‑1 (MQ=255)
                     gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggtgg                                                               >  2:1191523/1‑90 (MQ=255)
                     gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggtgg                                                               >  1:1185149/1‑90 (MQ=255)
                            gAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAAc                                                        >  2:1535859/1‑90 (MQ=255)
                                  aaGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACcgccgc                                                  <  2:1211111/90‑1 (MQ=255)
                                            gCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGc                                        >  2:305899/1‑90 (MQ=255)
                                            gCTCAGTGAACGTGATGAACACCTTCATATGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGc                                        >  1:191453/1‑90 (MQ=255)
                                              tCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                                      <  1:305899/90‑1 (MQ=255)
                                              tCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                                      >  1:815112/1‑90 (MQ=255)
                                              tCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGATGGGATAAACCGCCGCTGATTAACGccc                                      <  1:311500/90‑1 (MQ=255)
                                                 gTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGAtgg                                                                  <  2:1129650/59‑1 (MQ=255)
                                                 gTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGAtgg                                                                  >  1:1129650/1‑59 (MQ=255)
                                                 gTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCcgcg                                   >  2:1055500/1‑90 (MQ=255)
                                                    aaCGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAg                                >  2:1217078/1‑90 (MQ=255)
                                                      cGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                                      <  1:1164993/82‑1 (MQ=255)
                                                      cGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                                      >  2:1164993/1‑82 (MQ=255)
                                                      cGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGaa                              >  1:570398/1‑90 (MQ=255)
                                                           tgaACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGc                         >  2:661067/1‑90 (MQ=255)
                                                           tgaACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGc                         >  2:73443/1‑90 (MQ=255)
                                                           tgaACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGc                         >  1:269548/1‑90 (MQ=255)
                                                                 ccTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATaa                                                          <  2:653020/51‑1 (MQ=255)
                                                                 ccTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATaa                                                          >  1:653020/1‑51 (MQ=255)
                                                                 ccTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAAc                   >  1:1032259/1‑90 (MQ=255)
                                                                  cTTCATCTGGATCCGTTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAAcc                  <  1:23311/90‑1 (MQ=255)
                                                                  cTTCATCTGGATCCGGTTTCCTGGGGATATGCCCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAAcc                  <  2:1165523/90‑1 (MQ=255)
                                                                    tCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAg                <  1:890816/90‑1 (MQ=255)
                                                                      aTCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTc              >  2:462394/1‑90 (MQ=255)
                                                                        cTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGt            <  1:311442/90‑1 (MQ=255)
                                                                        cTGGATCCGGTTTCCTGGGGATATGCCCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGt            <  1:1196727/90‑1 (MQ=255)
                                                                           gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                                      <  2:1528723/61‑1 (MQ=255)
                                                                           gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                                      >  1:1528723/1‑61 (MQ=255)
                                                                           gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATg         >  2:719202/1‑90 (MQ=255)
                                                                            gaccGGTTTTCCGGGGGATATCCCCCGGGATGGTGGGAAAACCCGCCGCTGATTAACCCCCGCGTAGAAACTGCGGCAACCAGTCGTATg         <  1:830401/88‑1 (MQ=255)
                                                                              ccGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                                      >  2:1103370/1‑58 (MQ=255)
                                                                              ccGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                                      <  1:1103370/58‑1 (MQ=255)
                                                                                  tttcCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAAc  <  1:1535859/90‑1 (MQ=255)
                                                                                     |                                                                                       
CCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAAC  >  NZ_CP009273/2003395‑2003566

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