Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 118,882 G→T E102* (GAA→TAA)  pdhR → pyruvate dehydrogenase complex transcriptional repressor PdhR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273118,8820GT84.0% 62.3 / 8.5 25E102* (GAA→TAA) pdhRpyruvate dehydrogenase complex transcriptional repressor PdhR
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base T (9/12);  total (10/15)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

TTTTGTCCAGAGCAGCCTATGGCAAAGCTTCAGCGATCCGCTGGTGGAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCGAAACACGACACGCCCTGGAAGGTATCGCCGCTTATTACGCCGCGCTGCGTAGTACCGATGAAGACAAGGAACGCATCCGTGAACTCCAC  >  NZ_CP009273/118794‑118971
                                                                                        |                                                                                         
ttttGTCCAGAGCAGCCTATGGCAAAGCTTCAGCGATCCGCTGGTGGAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCGa                                                                                          <  1:756090/90‑1 (MQ=255)
       cAGAGCAGCCTATGGCAAAGCTTCAGCGATCCGCTGGTGGAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCTAAacacga                                                                                   <  1:125353/90‑1 (MQ=255)
         gagCAGCCTATGGCAAAGCTTCAGCGATCCGCTGGTGGAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCTAAacacgaca                                                                                 >  2:253951/1‑90 (MQ=255)
         gagCAGCCTATGGCAAAGCTTCAGCGATCCGCTGGTGGAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCTAAacacgaca                                                                                 <  1:762709/90‑1 (MQ=255)
          agcagcCTATGGCAAAGCTTCAGCGATCCGCTGGTGGAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCTAAacacgacac                                                                                <  2:739912/90‑1 (MQ=255)
          agcagcCTATGGCAAAGCTTCAGCGATCCGCTGGTGGAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCTAAacacgacac                                                                                <  2:521059/90‑1 (MQ=255)
          agcagcCTATGGCAAAGCTTCAGCGATCCGCTGGTGGAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCTAAacacgacac                                                                                <  1:194230/90‑1 (MQ=255)
           gcagcCTATGGCAAAGCTTCAGCGATCCGCTGGTGGAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCTAAacacgacacg                                                                               <  2:511362/90‑1 (MQ=255)
           gcagcCTATGGCAAAGCTTCAGCGATCCGCTGGTGGAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCTAAacacgacacg                                                                               <  2:337201/90‑1 (MQ=255)
              gcCTATGGCAAAGCTTCAGCGATCCGCTGGTGGAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCTAAACACGACACGccc                                                                            >  2:608134/1‑90 (MQ=255)
                 tATGGCAAAGCTTCAGCGATCCGCTGGTGGAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCTAAACACGACACGCCCTgg                                                                         <  1:253951/90‑1 (MQ=255)
                   tGGCAAAGCTTCAGCGATCCGCTGGTGGAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCTAAACACGACACGCCCTGGaa                                                                       >  2:119404/1‑90 (MQ=255)
                   tGGCAAAGCTTCAGCGATCCGCTGGTGGAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCTAAACACGACACGCCCTGGaa                                                                       >  1:110679/1‑90 (MQ=255)
                           cTTCAGCGATCCGCTGGTGGAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCTAAACACGACACGCCCTGGAAGGTATcgc                                                               <  2:101622/90‑1 (MQ=255)
                             tCAGCGATCCGCTGGTGGAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCTAAACACGACACGCCCTGGAAGGTATcgccg                                                             <  2:330263/90‑1 (MQ=255)
                                    tCCGCTGGTGGAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCTAAACACGACACGCCCTGGAAGGTATCGCCGCttatta                                                      >  1:232078/1‑90 (MQ=255)
                                            tggAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCTAAACACGACACGCCCTGGAAGGTATCGCCGCTTATTACGCCgcgc                                              <  2:232078/90‑1 (MQ=255)
                                            tggAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCGAAACACGACACGCCCTGGAAGGTATCGCCGCTTATTACGCCgcgc                                              <  2:277554/90‑1 (MQ=255)
                                            tggAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCGAAACACGACACGCCCTGGAAGGTATCGCCGCTTATTACGCCgcgc                                              <  2:277555/90‑1 (MQ=255)
                                                      ctcCGACCATCCTGAGTCACAGTATGACTTGCTCTAAACACGACACGCCCTGGAAGGTATCGCCGCTTATTACGCCGCGCTGCGTAGTAc                                    <  1:148576/90‑1 (MQ=255)
                                                              aTCCTGAGTCACAGTATGACTTGCTCTAAACACGACACGCCCTGGAAGGTATCGCCGCTTATTACGCCGCGCTGCGTAGTACCGATGAAg                            >  2:679044/1‑90 (MQ=255)
                                                                         caGTATGACTTGCTCTAAACACGACACGCCCTGGAAGGTATCGCCGCTTATTACGCCGCGCTGCGTAGTACCGATGAAGACAAGGAACGc                 >  1:420671/1‑90 (MQ=255)
                                                                         caGTATGACTTGCTCTAAACACGACACGCCCTGGAAGGTATCGCCGCTTATTACGCCGCGCTGCGTAGTACCGATGAAGACAAGGAACGc                 >  1:139908/1‑90 (MQ=255)
                                                                                   tGCTCTAAACACGACACGCCCTGGAAGGTATCGCCGCTTATTACGCCGCGCTGCGTAGTACCGATGAAGACAAGGAACGCATCCGTGAAc       >  2:193322/1‑90 (MQ=255)
                                                                                   tGCTCGAAACACGACACGCCCTGGAAGGTATCGCCGCTTATTACGCCGCGCTGCGTAGTACCGATGAAGACAAGGAACGCATCCGTGAAc       >  2:602126/1‑90 (MQ=255)
                                                                                        gAAACACGACACGCCCTGGAAGGTATCGCCGCTTATTACGCCGCGCTGCGTAGTACCGATGAAGACAAGGAACGCATCCGTGAACTccac  <  1:123649/90‑1 (MQ=255)
                                                                                        |                                                                                         
TTTTGTCCAGAGCAGCCTATGGCAAAGCTTCAGCGATCCGCTGGTGGAGCTGCTCTCCGACCATCCTGAGTCACAGTATGACTTGCTCGAAACACGACACGCCCTGGAAGGTATCGCCGCTTATTACGCCGCGCTGCGTAGTACCGATGAAGACAAGGAACGCATCCGTGAACTCCAC  >  NZ_CP009273/118794‑118971

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