Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,322,394 C→G 10.4% D185H (GAC→CAC)  trpE ← anthranilate synthase subunit TrpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,322,3940CG10.4% 108.0 / 3.5 44D185H (GAC→CAC) trpEanthranilate synthase subunit TrpE
Reads supporting (aligned to +/- strand):  ref base C (21/18);  new base G (2/3);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

GGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTC  >  NC_000913/1322252‑1322536
                                                                                                                                              |                                                                                                                                               
ggCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTc                                                                                                                                                 >  1:1184381/1‑143 (MQ=255)
 gCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGGGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCa                                                                                                                                                >  2:590773/1‑143 (MQ=255)
  gaCGGAAACCACGGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTAATGGTCaa                                                                                                                                               <  2:1054877/142‑1 (MQ=255)
  cACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCaa                                                                                                                                               <  2:218265/143‑1 (MQ=255)
          ccACTCACAGCCGCGGCGCCGCTTCTGTCCGTTGCTCACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTG‑TTTTCTTCGTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCa                                                                                                                                     <  1:959263/143‑1 (MQ=255)
             cTGGCAGCGGCGGCGCGGCGTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATT‑GACCATCAga                                                                                                                                   <  2:509608/143‑2 (MQ=255)
                               ctccGGTCAGTTGCTGACGTGGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGCGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATaaaaa                                                                                                                 <  2:1016473/140‑1 (MQ=255)
                               cTTCGGTGAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATaaaaa                                                                                                                 <  1:845203/143‑1 (MQ=255)
                               cTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGGTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGCTaaaaa                                                                                                                 >  1:1218110/1‑143 (MQ=255)
                                    gTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTTTGTAGGTCAAT‑CACCATCTGCGTTTCAGGGACATAAAAAAAGaa                                                                                                            >  2:413975/1‑143 (MQ=255)
                                           gCTGACGTGGTGCGTTCAGGCGAGCAGTGCGACGTTGTGTTTCTTCTTCATTCGGGTCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCAGCAGCGTTTCAGCGAGATAAAAACAGAAATCAggg                                                                                                     <  2:917136/143‑1 (MQ=255)
                                               aCGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGGGATAAAAACGAACATCATGGCTGt                                                                                                 >  2:138806/1‑143 (MQ=255)
                                                     ttcgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGGTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGTGAGCTAAACACAGGAATCAGGGCAGTTAtttt                                                                                           >  2:234595/1‑143 (MQ=255)
                                                     ttcgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTAtttt                                                                                           >  2:459424/1‑143 (MQ=255)
                                                     ttcgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTACTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAGCAGAAATCAGGGCAGTTCtttt                                                                                           >  2:260127/1‑143 (MQ=255)
                                                      tcgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTc                                                                                          <  1:1233334/143‑1 (MQ=255)
                                                        gttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCg                                                                                        <  1:851986/143‑1 (MQ=255)
                                                        gttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCg                                                                                        <  1:653630/143‑1 (MQ=255)
                                                          tcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCt                                                                                      >  1:612465/1‑143 (MQ=255)
                                                            aGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGa                                                                                    <  1:525638/143‑1 (MQ=255)
                                                                              tgtgTTTCTTCTTCATTCGGAGCATACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTc                                                                  <  2:939071/143‑1 (MQ=255)
                                                                                  tttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCCTCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAAt                                                              >  1:947519/1‑143 (MQ=255)
                                                                                  tttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTGTTATTC‑AATCGTGTTTGTGTCTTCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAAt                                                              <  1:590773/143‑1 (MQ=255)
                                                                                    tcttcctcATTCGGAGCAAACAGGCTGCCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAACACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATcc                                                            <  2:982827/143‑1 (MQ=255)
                                                                                      ttcttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGCCAGTTGCGGTAAATCTTCAACTCCCg                                                          >  1:491975/1‑143 (MQ=255)
                                                                                       tcttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCACATCCCGc                                                         >  2:772919/1‑143 (MQ=255)
                                                                                          tcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTGTCAGCGAGATGATAACAGAACCCAGCGCTTTTGTTTGCCGCTGACAGATGCGGTAAACATTCAAATACCGCcac                                                      >  1:47641/1‑143 (MQ=255)
                                                                                             ttCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAg                                                   <  2:1085746/143‑1 (MQ=255)
                                                                                                  aGCAAACAGGCTGGCCTGAATCTGGGTGCTTCTGTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCAt                                              <  2:1144693/143‑1 (MQ=255)
                                                                                                  aGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCAt                                              >  1:566157/1‑143 (MQ=255)
                                                                                                         aGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGaa                                       <  2:1210585/143‑1 (MQ=255)
                                                                                                            cTGGCCTGAATTCGCGGGCTCTTTTTCTGAGGGGGAGT‑CACCAGCAGCGTTGCACCGAGATAAAAACATAAATCGGTGCAGTTATTTCCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAg                                    <  2:1069535/143‑1 (MQ=255)
                                                                                                            cTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTTCGGTAAATCTTCCAATCCCGCCACCAGCACCTAATGGAACAg                                    >  2:886217/1‑143 (MQ=255)
                                                                                                                ccTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGCGAACAGgccg                                >  1:1112985/1‑143 (MQ=255)
                                                                                                                        cGGGTGCTTTTTTTCTGATGGTCAGT‑CACCATCAGCGTTTCAGCGAGCTAAAACCAGGCATCAGGGTAGTTATTTTCCGCTGTCCGTTGCGGTAACTGTTGAACTGCCGACTCAACGTCATAAGAGTACAGGTCACCgaggtc                        >  2:1017212/1‑139 (MQ=255)
                                                                                                                            tGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAACTCTTCAAATCACGCCACAAGGTCATAAGCGAACAGGCCTCCGAAGACCAtc                    >  2:1018417/1‑142 (MQ=255)
                                                                                                                             gCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATgg                   <  1:1070826/143‑1 (MQ=255)
                                                                                                                              cTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGAGAAAAACAGAAATCAGGGAAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATACGAGAACAGGCCGCAGAAGAACATGGc                  >  2:156074/1‑143 (MQ=255)
                                                                                                                               tttttttCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTTACAGTTGCGGTAAATCTTCAAATTCCGCCACAACGTCATACGAGAACAGGCTGCCGAATAACCTgacc                 >  2:124176/1‑140 (MQ=255)
                                                                                                                                   tttCTGATGGTCATT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACTGGCCGCCGAAGAACATGGCTtctc             <  2:1090985/143‑1 (MQ=255)
                                                                                                                                     tCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTCCGGTACATCTTCAAATCCCGCCACAAGCCCATCAGACACCAGGCCGCCGACAAACCTGGCCTCTCGt           >  2:620699/1‑143 (MQ=255)
                                                                                                                                     gcTGATGGTGATT‑GACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGt           >  1:1067156/2‑143 (MQ=255)
                                                                                                                                      cTGATGGTGATT‑GACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGtt          >  2:876442/1‑143 (MQ=255)
                                                                                                                                          tGGTCAAT‑CACCATCAGCGGTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGttcttc      <  2:967957/143‑1 (MQ=255)
                                                                                                                                           ggTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTcttcc     >  1:691664/1‑143 (MQ=255)
                                                                                                                                              cAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACCGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGCACATGGCTTGTCGTGctgcctcc  >  1:704230/1‑141 (MQ=255)
                                                                                                                                              |                                                                                                                                               
GGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTC  >  NC_000913/1322252‑1322536

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