Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 787,110 C→T 16.6% R146C (CGC→TGC)  pabB → anthranilate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730787,1100CT16.6% 51.2 / 5.1 27R146C (CGC→TGC) pabBanthranilate synthase
Reads supporting (aligned to +/- strand):  ref base C (11/11);  new base T (2/2);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.47e-01

ACATATGAACATCAAACTTCTTATGTATCAAAGTATCATTTTTCATTTGTTGAATCTTCTGAGGTTATGATGACTAATATTAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGCGCTTAACGGATAACATTTATTATCCTATTAGTACTTTATATGAACGATTAACTCAATTTAGTAATGGTAATTATACTGCGTTATTACAAACTGATGAAATCCAAGTAGCGTCTATCTCACCAGAATTATTTTTTCAAAAA  >  CP000730/786972‑787250
                                                                                                                                          |                                                                                                                                            
atatatGAACATCAAACTTCTTATGTATCAAATTATCATTTTTCATTTGGTGAATTTTCTGAGGTTATGATGACTAATATTAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGgcgcgc                                                                                                                                            <  2:24099/139‑1 (MQ=255)
    atGAACATCAAACTTCTTATGTATCAAAGTATCATTTTTCATTTGTTGAATCTTCTGAGGTTATGATGACTAATATTAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGCGCTTaa                                                                                                                                        >  2:80560/1‑141 (MQ=255)
      gAACATCAAACTTCGTATGTATCAAAGTATCATTTTTCATTTGTTGAATCTTCTGAGGTTATGATGACTAATATTAAACGTGTCTAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGCGCTTAACg                                                                                                                                      <  1:8924/141‑1 (MQ=255)
       aaCATCAAACTTCTTATGTATCAAAGTATCATTTTTCATTTGTTGAATCTTCTGAGGTTATGATGACTAATATTAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGCGCTTAACgg                                                                                                                                     >  1:281286/1‑141 (MQ=255)
        aCATCAAACTTCTTATGTATCAAAGTATCATTTTTCATTTGTTGAATCTTCTGAGGTTATGATGACTAATATTAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGCGCTTAACGGa                                                                                                                                    >  2:43263/1‑141 (MQ=255)
                    ttATGTATCAAAGTATCATTTTTCATTTGTTGAATCTTCTGAGGTTATGATGACTAATATTAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAGGTGAACTGTACGGCGCGCTTAACGGATAACATttatta                                                                                                                        >  1:177635/1‑141 (MQ=255)
                    ttATGTATCAAAGTATCATTTTTCATTTGTTGAATCTTCTGAGGTTATGATGACTAATATTAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATCCGGCGCGCTTAACGGATAACATttatta                                                                                                                        >  1:309510/1‑141 (MQ=255)
                           tCAAAGTATCATTTTTCATTTGTTGAATCTTCTGAGGTTATGATGACTAATATTAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATTAAGTGTACTATACGGCGCGCTTAACGGATAACATTTATTATCCTAtt                                                                                                                 <  2:226342/141‑1 (MQ=255)
                                  aTCATTTTTCATTTGTTGAATCTTCTGAGGTTATGATGACTAATATTAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGCGCTTAACGGATAACATTTATTATCCTATTAGTACtt                                                                                                          <  1:43263/141‑1 (MQ=255)
                                       ttttCATTTGTTGAATCTTCTGAGGTTATGATGACTAATATTAAACGTGTCCAACAAGCAATTATTGAGGGCGAAACGTATCAAGTGAACTATACGGCGGGCTTAACGGATAACATTTATTAGCCTATTAGTACTTTagat                                                                                                     >  1:191369/1‑141 (MQ=255)
                                              ttgttgAATCTTCTGAGGTTATGATGACTAATATTAAACGTGTCCAACAAGCAATTGTTGAGGGCGAAACGTATCAAGTGAACTATACGGCGCGCTCAACGGATAACATTTATTCTCCTATTAGTACTATTTGTGAACGAt                                                                                              >  2:249743/1‑141 (MQ=255)
                                               tgtttAATCTTCTGAGGTTATGATGACTAATATTAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGCGCTTAACGGATAACATTTATTATCCTATTAGTACTTTATATGAACGAtt                                                                                             <  1:295992/141‑1 (MQ=255)
                                                        ttctGAGGTTATGATGACTAATATTAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGCGCTTAACGGATAACATTTATTATCCTATTAGTACTTTATATGAACGATTAACTCAAtt                                                                                    <  1:150670/141‑1 (MQ=255)
                                                          ctGAGGTTATGATGACTAATATTAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGCGCTTAACGGATAACATTTATTATCCTATTAGAACTTGATATGAACGATTAACTCAATTTa                                                                                  >  2:314200/1‑141 (MQ=255)
                                                                 tATGATGACTAATATTAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGCGCTTAACGGATAACATTTATTATCCTATTAGTACTTTATATGAACGATTAACTCAATTTAGTAATgg                                                                           >  1:169477/1‑141 (MQ=255)
                                                                      tgaCTAATATTAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGCGCTTAACGGATAACATTTATTATCCTATTAGTACTTTATATGAACGATTAACTCAATTTAGTAATGGTAAtt                                                                      <  1:277700/141‑1 (MQ=255)
                                                                          tAATATTAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGCGCTTAACGGATAACATTTATTATCCTATTAGTACTTTATATGAACGATTAACTCAATTTAGTAATGGTAATTATAc                                                                  <  1:32631/141‑1 (MQ=255)
                                                                            atatTAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGTGCTTAACGGATAACATTTATTATCCTATTAGTACTTTATATGAACGATTAACTCAATTTAGTAATGGTAATTATACTg                                                                >  2:14754/1‑141 (MQ=255)
                                                                            atatTAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGTGCTTAACGGATAACATTTATTATCCTATTAGTACTTTATATGAACGATTAACTCAATTTAGTAATGGTAATTATACTg                                                                >  2:5281/1‑141 (MQ=255)
                                                                               atAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGCGCTTAACGGATAACATTTATTATCCTATTAGTACTTTATATGAACGATTAACTCAATTTAGTAATGGTAATTATACTGCGt                                                             <  1:313126/140‑1 (MQ=255)
                                                                                   aCGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGCGCTTAACGGATAACATTTCTTCTATATTTGGTACTTTATATGAACGATTAACTCAATTTAGTAATGGTAATTATACTGCGttatt                                                         >  1:79220/1‑141 (MQ=255)
                                                                                    ggtgtCCAACAAGCAATTGTTGAAGGGGAAATGGAACAAGTGAACTATCCGGCGCGCTTAACGGATAACATTTATTATCCTATTAGTACTTTATATGAACGATTAACTCAATTTAGTAATGGTAATTATACTGCGttatta                                                        <  2:121380/140‑1 (MQ=255)
                                                                                        tCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGCGCTTAACGGATAACATTTATTATCCTATTAGTACTTTATATGAACGATTAACTCAATTTAGTAATGGTAATTATACTGCGTTATTACaaa                                                    <  1:93576/141‑1 (MQ=255)
                                                                                                   aTTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGCGCTTAACGGATAACATTTATTATCCTATTAGTACTTTATATGAACGATTAACTCAATTTAGTAATGGTAATTATACTGCGTTATTACAAACTGATGAAATc                                         >  1:14345/1‑141 (MQ=255)
                                                                                                                          aGTGAACTATACGGCGCGCTTAACGGATAACATTTATTATCCTATTAGTACTTTATATGAACGATTAACTCAATTTAGTAATGGTAATTATACTGCGTTATTACAAACTGATGAAATCCAAGTAGCGTCTATCTCACCAGa                  <  2:102500/141‑1 (MQ=255)
                                                                                                                           gTGAACTATACGGCGTGCTTAACGGATAACATTTATTATCCTATTAGTACTTTATATGAACGATTAACTCAATTTAGTAATGGTAATTATACTGCGTTATTACAAACTGATGAAATCCAAGTAGCGTCTATCTCACCAGaa                 <  1:14754/141‑1 (MQ=255)
                                                                                                                           gTGAACTATACGGCGTGCTTAACGGATAACATTTATTATCCTATTAGTACTTTATATGAACGATTAACTCAATTTAGTAATGGTAATTATACTGCGTTATTACAAACTGATGAAATCCAAGTAGCGTCTATCTCACCAGaa                 <  1:5281/141‑1 (MQ=255)
                                                                                                                                   taCGGCGCGCTTAACGGATAACATTTATTATCCTATTAGTACTTTATATGAACGATTAACTCAATTTAGTAATGGTAATTATACTGCGTTATTACAAACTGATGAAATCCAAGTAGCGTCTATCTCACCACAATTAttttt         >  2:171978/1‑141 (MQ=255)
                                                                                                                                          cgcTTAACGGATAACATTTATTATCCTATTAGTACTTTATATGAACGATTAACTCAATTTAGTAATGGTAATTATACTGCGTTATTACAAACTGATGAAATCCAAGTAGCGTCTATCTCACCAGAATTATTTTTTCaaaaa  >  1:31859/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                            
ACATATGAACATCAAACTTCTTATGTATCAAAGTATCATTTTTCATTTGTTGAATCTTCTGAGGTTATGATGACTAATATTAAACGTGTCCAACAAGCAATTGTTGAAGGCGAAACGTATCAAGTGAACTATACGGCGCGCTTAACGGATAACATTTATTATCCTATTAGTACTTTATATGAACGATTAACTCAATTTAGTAATGGTAATTATACTGCGTTATTACAAACTGATGAAATCCAAGTAGCGTCTATCTCACCAGAATTATTTTTTCAAAAA  >  CP000730/786972‑787250

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