Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,322,396 A→T 11.9% I184N (ATT→AAT)  trpE ← anthranilate synthase subunit TrpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,322,3960AT11.9% 96.3 / 3.6 42I184N (ATT→AAT) trpEanthranilate synthase subunit TrpE
Reads supporting (aligned to +/- strand):  ref base A (20/15);  new base T (2/3);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-01

CACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGG  >  NC_000913/1322254‑1322538
                                                                                                                                              |                                                                                                                                               
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:653630/143‑1 (MQ=255)
                                                      gttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCg                                                                                          <  1:851986/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)
                                                                                                                                             aaT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGACGAACATGGCTTCTCGTTCTTCCTTCg   >  2:190324/1‑143 (MQ=255)
                                                                                                                                              aT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGACGAAGAACATGGCTTCTCGTTCTTCCTTCgg  >  1:545617/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                               
CACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAT‑CACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGG  >  NC_000913/1322254‑1322538

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: ‑