Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 98,192 T→C I314T (ATT→ACT)  murC → UDP‑N‑acetylmuramate‑‑L‑alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927398,1920TC98.4% 194.0 / ‑3.9 61I314T (ATT→ACT) murCUDP‑N‑acetylmuramate‑‑L‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (36/24);  total (36/25)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-01

GAATGCGCCAGGTCGTCATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTATTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGCCAGTGAATGGTAAAA  >  NZ_CP009273/98109‑98278
                                                                                   |                                                                                      
gAATGCGCCAGGTCGTCATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCg                                                                                  >  2:989347/1‑90 (MQ=255)
gAATGCGCCAGGTCATCATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCg                                                                                  >  1:1331028/1‑90 (MQ=255)
     cgcCAGGTCGTCATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGgcgc                                                                             >  1:549693/1‑90 (MQ=255)
     cgcCAGGTCGTCATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGgcgc                                                                             >  1:143926/1‑90 (MQ=255)
      gcCAGGTCGTCATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCt                                                                            >  2:214778/1‑90 (MQ=255)
      gcCAGGTCGTCATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCt                                                                            >  2:459048/1‑90 (MQ=255)
      gcCAGGTCGTCATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCt                                                                            >  2:887509/1‑90 (MQ=255)
          ggTCGTCATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGaa                                                                        >  1:1331393/1‑90 (MQ=255)
          ggTCGTCATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGaa                                                                        >  1:505762/1‑90 (MQ=255)
          ggTCGTCATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGaa                                                                        >  1:438214/1‑90 (MQ=255)
           gtcgtcATACCGCGCCGAACGCCGCAGCTGCGGTTGCGGTTGGTGCGGAAGAGGGCATTGACGACGAGGGAACTGTGCGGGCGCTTGaaa                                                                       <  1:442451/90‑1 (MQ=255)
            tcgtcATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAg                                                                      <  1:1188366/90‑1 (MQ=255)
            tcgtcATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAg                                                                      <  1:1169684/90‑1 (MQ=255)
             cgtcATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTATTTTGCGGGCGCTTGAAAGc                                                                     <  1:399543/90‑1 (MQ=255)
              gtcATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCt                                                                    <  2:640226/90‑1 (MQ=255)
              gtcATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCt                                                                    <  1:564606/90‑1 (MQ=255)
              gtcATAACGCGCTGAACGCCGCAGATGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCt                                                                    <  2:1053971/90‑1 (MQ=255)
               tcATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCtt                                                                   >  2:1465970/1‑90 (MQ=255)
                   aaCGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAg                                                               <  2:657791/90‑1 (MQ=255)
                   aaCGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAg                                                               <  1:821223/90‑1 (MQ=255)
                    aCGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAgg                                                              >  1:1073404/1‑90 (MQ=255)
                    aCGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAgg                                                              >  1:948881/1‑90 (MQ=255)
                        gcTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGAc                                                          >  2:745164/1‑90 (MQ=255)
                        gcTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGAc                                                          >  1:206073/1‑90 (MQ=255)
                        gcTGAACGCCGCAGCTGCGGTTGCGGGTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGAc                                                          >  1:866987/1‑90 (MQ=255)
                         cTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACt                                                         >  1:659808/1‑90 (MQ=255)
                               gccgcAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGc                                                   <  1:1465970/90‑1 (MQ=255)
                                    aGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGttt                                              <  2:12945/90‑1 (MQ=255)
                                     gCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCa                                                                <  1:400797/71‑1 (MQ=255)
                                     gCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCa                                                                >  2:400797/1‑71 (MQ=255)
                                     gCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGtttt                                             >  2:657355/1‑90 (MQ=255)
                                         cggttgcggttgCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGAtt                                         <  1:946653/90‑1 (MQ=255)
                                         cggttgcggttgCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGAtt                                         <  2:21088/90‑1 (MQ=255)
                                          ggttgcggttgcTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGAttt                                        >  2:548341/1‑90 (MQ=255)
                                              gcggttgcTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTc                                    <  2:549693/90‑1 (MQ=255)
                                              gcggttgcTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTc                                    <  2:705499/90‑1 (MQ=255)
                                                ggttgcTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCgg                                  <  2:1367151/90‑1 (MQ=255)
                                                               gggCATTGACGACGAGGCTACTTTGCGGGCGCTTGaa                                                                        >  2:148185/1‑37 (MQ=38)
                                                               gggCATTGACGACGAGGCTACTTTGCGGGCGCTTGaa                                                                        <  1:148185/37‑1 (MQ=38)
                                                               gggCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTTCCGCTGGa                   >  1:367340/1‑90 (MQ=255)
                                                               gggCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGa                   >  1:1034042/1‑90 (MQ=255)
                                                               gggCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGa                   >  2:905144/1‑90 (MQ=255)
                                                               gggCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGa                   >  1:1100429/1‑90 (MQ=255)
                                                               gggCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGa                   >  1:977945/1‑90 (MQ=255)
                                                               gggCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGa                   >  2:478754/1‑90 (MQ=255)
                                                               gggCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGa                   >  1:920652/1‑90 (MQ=255)
                                                               gggCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGa                   >  1:56712/1‑90 (MQ=255)
                                                               gggCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGa                   >  1:291097/1‑90 (MQ=255)
                                                                  cATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGcc                <  1:1388358/90‑1 (MQ=255)
                                                                   aTTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGCCa               <  1:143489/90‑1 (MQ=255)
                                                                    ttGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGa                               <  1:417933/73‑1 (MQ=255)
                                                                    ttGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGa                               >  2:417933/1‑73 (MQ=255)
                                                                    ttGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGCCAg              >  2:1192839/1‑90 (MQ=255)
                                                                    ttGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGCCAg              >  2:1150065/1‑90 (MQ=255)
                                                                         gacgaGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGCCAGTGAAt         >  1:977832/1‑90 (MQ=255)
                                                                          acgaGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGc                 <  1:1323756/81‑1 (MQ=255)
                                                                          acgaGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGc                 >  2:1323756/1‑81 (MQ=255)
                                                                             aGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGCCAGTGAATGGTa     >  1:188985/1‑90 (MQ=255)
                                                                               gCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGCCAGTGAATGGTaaa   <  2:66745/90‑1 (MQ=255)
                                                                                cTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGCCAGTGAATGGTaaaa  <  2:626647/90‑1 (MQ=255)
                                                                                cTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGCCAGTGAATGGTaaaa  <  2:844346/90‑1 (MQ=255)
                                                                                   |                                                                                      
GAATGCGCCAGGTCGTCATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTATTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGCCAGTGAATGGTAAAA  >  NZ_CP009273/98109‑98278

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