Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,328,749 A→G S459P (TCC→CCC)  yfaA ← DUF2138 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,328,7490AG95.7% 147.8 / ‑0.7 46S459P (TCC→CCC) yfaADUF2138 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (28/16);  total (29/17)
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.87e-01

CTATTCCTTGTGGATTGATATAGAGCGGAACGATGCCATCAGTGGGTATTACATCCACCATTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGAGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGCGCCGCCTGCGCTTTCGGATACTGGCCGTATCGGGAACTCACTTCGCGACGCCAGATTTGCG  >  NZ_CP009273/2328629‑2328883
                                                                                                                        |                                                                                                                                      
cTATTCCTTGTGGATTGATATAGAGCGGAACGATGCCATCAGTGGGTATTACATCCACCATTGCCGGGC‑GGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCgttttgttt                                                                                                                     >  2:668573/1‑139 (MQ=255)
    tCCTTGTGGATTGATATAGAGCGGAACGATGCCATCAGTGGGTATTACATCCACCATTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGttttgttttg                                                                                                                   >  1:357260/1‑138 (MQ=255)
      cTTGTGGATTGATATAGAGCGGAACGATGCCATCAGTGGGTATTACATCCACCATTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCa                                                                                                                 <  2:981703/138‑1 (MQ=255)
        tgtgGATTGATATAGAGCGGAACGATGCCATCAGTGGGTATTACATCCACCATTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAagca                                                                                                                                <  1:1067706/121‑1 (MQ=255)
        tgtgGATTGATATAGAGCGGAACGATGCCATCAGTGGGTATTACATCCACCATTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAagca                                                                                                                                >  2:1067706/1‑121 (MQ=255)
             aTTGATATAGAGCGGAACGATGCCATCAGTGGGTATTACATCCACCATTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGc                                                                                                         >  2:572869/1‑139 (MQ=255)
                 gtatcGAGGGGTACGATGCCATCAGTGGGTATTACATCCACCATTGCCGGGCGGGTTTTATTCAGTGCTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGaca                                                                                                     >  2:286934/2‑139 (MQ=255)
                         cGGAACGATGCCATCAGTGGGTATTACATCCACCATTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGaa                                                                                              >  1:745041/1‑138 (MQ=255)
                         cGGAACGATGCCATCAGTGGGTATTACATCCACCATTGCCGGGC‑GGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGaaaa                                                                                            <  1:668573/139‑1 (MQ=255)
                            aaCGATGCCATCAGTGGGTATTACATCCACCATTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGaaaaaa                                                                                          >  1:1048961/1‑139 (MQ=255)
                            aaCGATGCCATCAGTGGGTATTACATCCACCATTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGaaaaaa                                                                                          >  1:330811/1‑139 (MQ=255)
                            aaCGATGCCATCAGTGGGTATTACATCCACCATTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGaaaaaa                                                                                          <  2:1070950/139‑1 (MQ=255)
                                aTGCCATCAGTGGGTATTACATCCACCATTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAaataat                                                                                      <  1:360595/139‑1 (MQ=255)
                                          tGGGTATTACATCCACCATTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTaa                                                                            >  2:229597/1‑139 (MQ=255)
                                               aTTACATCCACCATTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGAt                                                                       >  2:668132/1‑139 (MQ=255)
                                                   cATCCACCATTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCgg                                                                    <  2:357260/138‑1 (MQ=255)
                                                      ccaccaTTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGaca                                                                                                     <  1:200516/102‑1 (MQ=255)
                                                      ccaccaTTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGaca                                                                                                     >  2:200516/1‑102 (MQ=255)
                                                        accaTTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATcc                                                                                    <  2:824663/117‑1 (MQ=255)
                                                        accaTTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATcc                                                                                    >  1:824663/1‑117 (MQ=255)
                                                        accaTTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTgcg                                                              <  1:229597/139‑1 (MQ=255)
                                                                 gggcgggTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGCGCCGCCTgcg                                                     <  1:1039619/139‑1 (MQ=255)
                                                                 gggcgggTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGCGCCGCCTgcg                                                     <  1:720810/139‑1 (MQ=255)
                                                                 gggcgggTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGAGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGCGCCGCCTgcg                                                     <  2:471380/139‑1 (MQ=255)
                                                                                 gtgtTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGCGCCGCCTGCGCTTTCGGATACTGGcc                                     >  1:732040/1‑139 (MQ=255)
                                                                                 gtgtTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGCGCCGCCTGCGCTTTCGGATACTGGcc                                     <  1:749117/139‑1 (MQ=255)
                                                                                 gtgtTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGCGCCGCCTGCGCTTTCGGATACTGGcc                                     >  2:1056997/1‑139 (MQ=255)
                                                                                 gtgtTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGAGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGCGCCGCCTGCGCTTTCGGATACTGGcc                                     >  1:559278/1‑139 (MQ=255)
                                                                                      tGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGa                                                                                               <  2:472826/76‑1 (MQ=255)
                                                                                      tGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGa                                                                                               >  1:472826/1‑76 (MQ=255)
                                                                                            gcgTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCa                                                                                                                 >  2:70515/1‑52 (MQ=255)
                                                                                            gcgTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCa                                                                                                                 <  1:70515/52‑1 (MQ=255)
                                                                                            gcgTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAaataat                                                                                      <  1:1201461/79‑1 (MQ=255)
                                                                                            gcgTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAaataat                                                                                      >  2:1201461/1‑79 (MQ=255)
                                                                                            gcgTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGCGCCGCCTGCGCTTTCGGATACTGGCCGTATCGGGAAc                          >  1:581480/1‑139 (MQ=255)
                                                                                            gcgTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGCGCCGCCTGCGCTTTCGGATACTGGCCGTATCGGGAAc                          >  2:324690/1‑139 (MQ=255)
                                                                                            gcgTTATTAACCAGCGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGCGCCGCCTGCGCTTTCGGATACTGGCCGTATCGGGAAc                          >  1:177746/1‑139 (MQ=255)
                                                                                                          cGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGCGCCGCCTGCGCTTTCGGATACTGGCCGTATCGGGAACTCACTTCGCGACGc            >  1:1187281/1‑139 (MQ=255)
                                                                                                          cGTGTCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGCGCCGCCTGCGCTTTCGGATACTGGCCGTATCGGGAACTCACTTCGCGACGc            >  1:821009/1‑139 (MQ=255)
                                                                                                              tCATCGAGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGCGCCGCCTGCGCTTTCGGATACTGGCCGTATCGGGAACTCACTTCGCGACGc            >  1:324426/1‑135 (MQ=255)
                                                                                                                    aGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAaataat                                                                                      >  1:645720/1‑55 (MQ=255)
                                                                                                                    aGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAaataat                                                                                      <  2:645720/55‑1 (MQ=255)
                                                                                                                    aGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGCGCCGCCTGCGCTTTCGGATACTGGCCGTATCGGGAACTCACTTCGCGACGCCAGATTTgcg  >  1:1227346/1‑139 (MQ=255)
                                                                                                                    aGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGCGCCGCCTGCGCTTTCGGATACTGGCCGTATCGGGAACTCACTTCGCGACGCCAGATTTgcg  >  2:753583/1‑139 (MQ=255)
                                                                                                                    aGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGCGCCGCCTGCGCTTTCGGATACTGGCCGTATCGGGAACTCACTTCGCGACGCCAGATTTgcg  >  1:261686/1‑139 (MQ=255)
                                                                                                                    aGGGGGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGAGCCGCCTGCGCTTTCGGATACTGGCCGTATCGGGAACTCACTTCGCGACGCCAGATTTgcg  >  1:1103744/1‑139 (MQ=255)
                                                                                                                        |                                                                                                                                      
CTATTCCTTGTGGATTGATATAGAGCGGAACGATGCCATCAGTGGGTATTACATCCACCATTGCCGGGCGGGTTTTATTCAGTGTTTGCAGCGCGTTATTAACCAGCGTGTCATCGAGGGAGAAAAGCAGCGTTTTGTTTTGCATCGCCAGCGACACACGGAAAAAATAATCCGACATTAATTGATCGGGTTGCGCCGCCTGCGCTTTCGGATACTGGCCGTATCGGGAACTCACTTCGCGACGCCAGATTTGCG  >  NZ_CP009273/2328629‑2328883

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑