Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,321,474 A→T 7.8% V491V (GTT→GTA)  trpE ← component I of anthranilate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,321,4740AT 7.8% 146.0 / 3.0 51V491V (GTT→GTA) trpEcomponent I of anthranilate synthase
Reads supporting (aligned to +/- strand):  ref base A (22/25);  new base T (2/2);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

TGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTT  >  NC_000913/1321326‑1321622
                                                                                                                                                    |                                                                                                                                                    
tGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGa                                                                                                                                                      <  8:137276/149‑1 (MQ=255)
     tGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAg                                                                                                                                                  >  8:257335/1‑148 (MQ=255)
     tGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCAATCAGAAAGTCTCCTGTGCATGATGCGCGCTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCCACTGCGGAACAg                                                                                                                                                  >  6:83718/1‑148 (MQ=255)
       ccAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGTGGAACAGAAt                                                                                                                                               >  6:78845/1‑149 (MQ=255)
        cAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATATCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAAt                                                                                                                                               >  7:248293/1‑148 (MQ=255)
         aGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATc                                                                                                                                              >  6:72535/1‑148 (MQ=255)
          ggTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCATTAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCaa                                                                                                                                            <  3:121999/149‑1 (MQ=255)
           gTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAg                                                                                                                                           <  5:248584/149‑1 (MQ=255)
           gTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAg                                                                                                                                           <  6:326056/149‑1 (MQ=255)
           gTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAACGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAg                                                                                                                                           >  5:49151/1‑149 (MQ=255)
           gTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAACGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAg                                                                                                                                           <  6:49151/149‑1 (MQ=255)
                 cGTAAAAGAGGCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGGCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTa                                                                                                                                      <  6:295832/148‑1 (MQ=255)
                   tAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTacac                                                                                                                                   <  7:257335/149‑1 (MQ=255)
                     aaaGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTaccc                                                                                                                                   >  8:188019/1‑145 (MQ=255)
                        gagTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCa                                                                                                                              >  2:331737/1‑149 (MQ=255)
                                    tCGAGCAGCAGAATGTCAGCCATCAGAAAGTCGCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACg                                                                                                                  >  7:57356/1‑149 (MQ=255)
                                      gAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCaa                                                                                                                                            >  2:385595/1‑121 (MQ=255)
                                      gAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCaa                                                                                                                                            <  1:385595/121‑1 (MQ=255)
                                      gAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGGTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTaca                                                                                                                                    >  5:278538/1‑129 (MQ=255)
                                          agcagAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGTACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGc                                                                                                             >  6:77642/1‑148 (MQ=255)
                                          agcagAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGTACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGc                                                                                                             >  6:72247/1‑148 (MQ=255)
                                           gcagAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGa                                                                                                           <  5:72535/149‑1 (MQ=255)
                                             agAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATa                                                                                                         <  1:331737/149‑1 (MQ=255)
                                                   tCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGTACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGttt                                                                                                   <  5:77642/149‑1 (MQ=255)
                                                   tCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGTACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGttt                                                                                                   <  5:72247/149‑1 (MQ=255)
                                                     aGCCATCAGAAAGGCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTc                                                                                                 <  1:74620/149‑1 (MQ=255)
                                                      gCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTcc                                                                                                <  8:248293/149‑1 (MQ=255)
                                                         aTCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCaa                                                                                                                                            >  8:145459/1‑102 (MQ=255)
                                                         aTCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTccacc                                                                                             <  7:188019/149‑1 (MQ=255)
                                                         aTCAGAAAGTCTCCTGTGCATGATGCGCGGGGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCaa                                                                                                                                            <  7:145459/102‑1 (MQ=255)
                                                            aGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGTCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTac                                                                                                                                     <  7:14143/106‑1 (MQ=255)
                                                                  tctcCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATa                                                                                                         <  4:59979/128‑1 (MQ=255)
                                                                  tctcCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATa                                                                                                         >  3:59979/1‑128 (MQ=255)
                                                                    tcCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCg                                                                                  >  4:346629/1‑149 (MQ=255)
                                                                       tgtgCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGa                                                                                >  3:240650/1‑148 (MQ=255)
                                                                         tgCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGa                                                                                <  4:240650/146‑1 (MQ=255)
                                                                         tgCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGAt                                                                               <  2:249247/147‑1 (MQ=255)
                                                                         tgCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATc                                                                              <  7:285513/148‑1 (MQ=255)
                                                                           cATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATcaca                                                                           <  3:346629/149‑1 (MQ=255)
                                                                            atgatgCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGTACTACACCAGCACCCGCTTGCAGGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCaaaa                                                                          <  6:30125/149‑1 (MQ=255)
                                                                                tgCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATcaca                                                                           <  4:107806/144‑1 (MQ=255)
                                                                                tgCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATcaca                                                                           >  3:107806/1‑144 (MQ=255)
                                                                                               gcgcgcAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGc                                                        <  8:57356/148‑1 (MQ=255)
                                                                                                 gcgcAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCAt                                                     <  5:83718/149‑1 (MQ=255)
                                                                                                      gTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGcgcg                                                <  2:109430/149‑1 (MQ=255)
                                                                                                       tACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACa                                                                                                                                                   >  1:250398/1‑49 (MQ=255)
                                                                                                                     tGTTACGGGTTTCGTCGGCTTCCGACTGTGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTAc                                  <  5:78845/148‑1 (MQ=255)
                                                                                                                             gTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCCGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGCGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCgccgcc                         >  6:13200/1‑149 (MQ=255)
                                                                                                                                   tcgGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCt                   >  3:165423/1‑149 (MQ=255)
                                                                                                                                          ccGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCgcgg            >  5:369110/1‑149 (MQ=255)
                                                                                                                                             aCTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGcga         >  5:121406/1‑149 (MQ=255)
                                                                                                                                                  ggAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGAcc    >  3:152018/1‑149 (MQ=255)
                                                                                                                                                    aaCAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCtt  >  7:363671/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTT  >  NC_000913/1321326‑1321622

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑