Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,153,219 G→A 12.5% intergenic (+39/‑81) fabF → / → pabC 3‑oxoacyl‑[acyl‑carrier‑protein] synthase II/4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,153,2190GA12.5% 129.7 / 8.7 48intergenic (+39/‑81)fabF/pabC3‑oxoacyl‑[acyl‑carrier‑protein] synthase II/4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex
Reads supporting (aligned to +/- strand):  ref base G (20/22);  new base A (3/3);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

TGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCGGGCAACGCAGTTTGGTGA  >  NC_000913/1153071‑1153367
                                                                                                                                                    |                                                                                                                                                    
tGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCggg                                                                                                                                                      >  1:51997/1‑149 (MQ=255)
    tttCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCtt                                                                                                                                                  >  1:346371/1‑149 (MQ=255)
      tCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGCCCGCAAGCGGACCtttt                                                                                                                                                <  1:382153/149‑1 (MQ=255)
         tACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCtttttt                                                                                                                                              <  4:244974/148‑1 (MQ=255)
         tACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCtttttt                                                                                                                                              >  4:103809/1‑148 (MQ=255)
          aCCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGCCCGCAAGCGGACCttttata                                                                                                                                             <  3:245776/148‑4 (MQ=255)
            cGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGCCCGCAAGCGGACCtttt                                                                                                                                                >  2:270012/1‑143 (MQ=255)
            cGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGCCCGCAAGCGGACCtttt                                                                                                                                                <  1:270012/143‑1 (MQ=255)
                cGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCtt                                                                                                                                      <  2:331722/149‑1 (MQ=255)
                   aGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTa                                                                                                                                   <  1:295515/149‑1 (MQ=255)
                         tCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGa                                                                                                                             <  2:19365/149‑1 (MQ=255)
                                              cacTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACAtg                                                                                                         <  3:103809/148‑1 (MQ=255)
                                                         aCTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTaa                                                                                             <  2:51997/149‑1 (MQ=255)
                                                            ccTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTg                                                                                                                         <  3:188248/118‑1 (MQ=255)
                                                            ccTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTg                                                                                                                         >  4:188248/1‑118 (MQ=255)
                                                                 ggcttcggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATtctc                                                                                     <  1:210399/149‑1 (MQ=255)
                                                                  gcttcggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTctct                                                                                    >  2:71437/1‑149 (MQ=255)
                                                                        gTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATaa                                                                               <  4:383850/148‑1 (MQ=255)
                                                                          ggCACTAATGGTTCTTTAATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTCGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCCCTGATAAGGa                                                                            <  2:59370/149‑1 (MQ=255)
                                                                                    gTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTctct                                                                                    <  4:3635/131‑1 (MQ=255)
                                                                                    gTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTctct                                                                                    >  3:3635/1‑131 (MQ=255)
                                                                                     ttCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGtt                                                                 >  2:329708/1‑149 (MQ=255)
                                                                                      tCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTc                                                                <  1:123385/149‑1 (MQ=255)
                                                                                      tCTTTGATCTTTAAAAAGACCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTc                                                                >  1:70022/1‑149 (MQ=255)
                                                                                        tttGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTaa                                                                                             >  1:89380/1‑118 (MQ=255)
                                                                                        tttGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTaa                                                                                             <  2:89380/118‑1 (MQ=255)
                                                                                        tttGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCtt                                                              >  1:299236/1‑149 (MQ=255)
                                                                                               tttAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGtt                                                                                                                       >  2:251401/1‑85 (MQ=255)
                                                                                               tttAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGtt                                                                                                                       <  1:251401/85‑1 (MQ=255)
                                                                                               tttAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGaa                                                                                                             <  2:41154/95‑1 (MQ=255)
                                                                                               tttAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGaa                                                                                                             >  1:41154/1‑95 (MQ=255)
                                                                                                     aaGACCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCAt                                                 <  2:70022/149‑1 (MQ=255)
                                                                                                       gATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATa                                                <  3:113124/148‑1 (MQ=255)
                                                                                                         tCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCGTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAAGATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAg                                              >  4:248756/1‑148 (MQ=255)
                                                                                                            aaGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAgg                                          >  2:373069/1‑149 (MQ=255)
                                                                                                            aaGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAg                                           >  4:133071/1‑148 (MQ=255)
                                                                                                             aGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGa                                         <  2:346371/149‑1 (MQ=255)
                                                                                                                tGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAAt                                       <  3:248756/148‑1 (MQ=255)
                                                                                                                      tttCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCTGACTTGTTCCGTAGTGAACATGCGTCCACACTATCAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTgg                                 >  4:184151/1‑148 (MQ=255)
                                                                                                                           aCCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTaa                           >  1:346317/1‑149 (MQ=255)
                                                                                                                              cTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAAc                                                       >  2:271149/1‑118 (MQ=255)
                                                                                                                              cTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAAc                                                       <  1:271149/118‑1 (MQ=255)
                                                                                                                              cTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGc                         >  4:20919/1‑148 (MQ=255)
                                                                                                                                   aaaaGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCgt                   <  2:346317/149‑2 (MQ=255)
                                                                                                                                   aaaaGGCCCGCAAGCGGACCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCgg                   >  2:77870/1‑149 (MQ=255)
                                                                                                                                    aaaGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCggg                  <  1:253118/149‑1 (MQ=255)
                                                                                                                                    aaaGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCggg                  >  2:348343/1‑149 (MQ=255)
                                                                                                                                     aaGGCCCGCAAGCGGACCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTa                                                                     >  4:324227/1‑97 (MQ=255)
                                                                                                                                         tCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCGGGCatc              >  3:242742/1‑146 (MQ=255)
                                                                                                                                                    gCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTTTCTGATAAGGAGCGGGTATGTTCTTAATTAACGGTCATAAGCAGGAATAGCTGGCAGTAAGCGATCGGGCAACGCAGTTTGgaga  >  2:389697/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCGGGCAACGCAGTTTGGTGA  >  NC_000913/1153071‑1153367

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: