Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,451,070 T→C 9.8% K116K (AAA→AAG hutU ← urocanate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,451,0700TC 9.8% 110.2 / 5.6 41K116K (AAA→AAGhutUurocanate hydratase
Reads supporting (aligned to +/- strand):  ref base T (22/15);  new base C (2/2);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

AACAGTGCCTGCTAAATCTCCATTAAAGTGTTGATTACCTAGCTCGGCAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGGCGCTTCTTCATGTGTTTTAAATACAGCAACTGGCTTTCCTGATTGTACTAACAT  >  CP000730/2450935‑2451201
                                                                                                                                       |                                                                                                                                   
aaCAGTGCCTGCTAAATCTCCATTAAAGTGTTGATTACCTAGCTCGGCAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCtttttt                                                                                                                                <  1:226704/141‑1 (MQ=255)
aaCAGTGCCTGCTAAACCGACAGTAAAGTGTTGATTCCCTAGCTCGGCAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCtttttt                                                                                                                                <  1:350429/141‑1 (MQ=255)
       ccTGCTAAATCTCCATTAAAGTGTTGATTACCTAGCTCGGCAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAAt                                                                                                                         >  2:427464/1‑141 (MQ=255)
       ccTGCTAAATCTCCATTAAAGTGTTGATTACCTAGCTCGGCAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAAt                                                                                                                         >  2:136006/1‑141 (MQ=255)
       ccTGCTAAATCTCCATTAAAGTGTTGATTACCTAGCTCGGCAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCCTTTTTATCTAAt                                                                                                                         >  1:347372/1‑141 (MQ=255)
       ccTGCTAAATCTCCATTAAAGTGTTGATTACCTAGCTAGGCAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCCTTTTTATCTAAt                                                                                                                         >  1:196373/1‑141 (MQ=255)
       ccTGCTAAATCTCCAGGAAATTGTTGATTACCTAGCTCGGCAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAAt                                                                                                                         >  2:94634/1‑141 (MQ=255)
        cTGCTAAATCTCCATTAAAGTGTTGATTACCTAGCTCGGCAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAAtt                                                                                                                        <  2:13815/141‑1 (MQ=255)
            tAAATCTCCATTAAAGTGTTGATTACCTAGCTCGGCAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCAtt                                                                                                                    <  2:303746/141‑1 (MQ=255)
            tAAATCTCCATTAAAGTGTTGATTACCTAGCTCGCCAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCAtt                                                                                                                    <  1:426414/141‑1 (MQ=255)
                   ccATTAAAGTGTTGATTACCTAGCTCGGCAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGa                                                                                                             <  2:259318/141‑1 (MQ=255)
                    cATTAAAGTGTTGATTACCTAGCTCGGCAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGAt                                                                                                            >  2:194906/1‑141 (MQ=255)
                    cATTAAAGTGTTGATTACCTAGCTCGGCAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGAt                                                                                                            >  2:291869/1‑141 (MQ=255)
                      ttAAAGTGTTGATTACCTAGCTCGGCAAATGTCTCATAAGTACCTTGCACAATACCTTGATAACCAATATATATCCATCTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATcc                                                                                                          <  1:230041/141‑1 (MQ=255)
                      ttAAAGTGTTGATTACCTAGCTCGGCAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATcc                                                                                                          <  1:290244/141‑1 (MQ=255)
                                gATTACCTAGCTCGGCAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGccc                                                                                                >  2:194417/1‑141 (MQ=255)
                                    aCCTAGCTCGGCAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTc                                                                                            >  1:96245/1‑141 (MQ=255)
                                               cAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATa                                                                                 <  1:444568/141‑1 (MQ=255)
                                               cAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATa                                                                                 <  1:324862/141‑1 (MQ=255)
                                               cAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATa                                                                                 <  1:94720/141‑1 (MQ=255)
                                                    gTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGaa                                                                            >  1:328/1‑141 (MQ=255)
                                                       tcATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGGTTGCCCATTCAGGCACTAATACTGAAttt                                                                         >  1:454496/1‑141 (MQ=255)
                                                        cATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTg                                                                        <  2:25359/141‑1 (MQ=255)
                                                           aaGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGaaa                                                                     >  2:316024/1‑141 (MQ=255)
                                                           aaGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGaaa                                                                     >  2:334129/1‑141 (MQ=255)
                                                                           aCCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGgcgc                                                     >  1:393911/1‑141 (MQ=255)
                                                                           aCCTTGAGAACCAATATATATCCAACTACCGGCTGACATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGgcgc                                                     >  1:401471/1‑141 (MQ=255)
                                                                                   aaCCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGGCGCTTCTTCAt                                             <  1:201562/141‑1 (MQ=255)
                                                                                    aCCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGGCGCTTCTTCAtg                                            <  1:238829/141‑1 (MQ=255)
                                                                                                  aaCTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGGCGCTTCTTCATGTGTTTTAAATACAg                              >  1:202666/1‑141 (MQ=255)
                                                                                                  aaCTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGGCGCTTCTTCATGTGTTTTAAATACAg                              >  1:4698/1‑141 (MQ=255)
                                                                                                       ccGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGGCGCTTCTTCATGTGTTTTAAATACAGCAACt                         >  2:108592/1‑141 (MQ=255)
                                                                                                       ccGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGGCGCTTCTTCATGTGTTTTAAATACAGCAACt                         >  1:281487/1‑141 (MQ=255)
                                                                                                       ccGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGGCGCTTCTTCATGTGTTTTAAATACAGCAACt                         >  2:48551/1‑141 (MQ=255)
                                                                                                       ccGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTAAGGCACTAATACTGAATTTGAAATTAATCCACGTGGAGATTCCTCATGTGTTTTAAATACAGCACCt                         >  1:20311/1‑141 (MQ=255)
                                                                                                        cGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGGCGCTTCTTCATGTGTTTTAAATACAGCAACTg                        >  2:122955/1‑141 (MQ=255)
                                                                                                        cGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGGCGCTTCTTCATGTGTTTTAAATACAGCAACTg                        >  2:9312/1‑141 (MQ=255)
                                                                                                                         tACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGGCGCTTCTTCATGTGTTTTAAATACAGCAACTGGCTTTCCTGATTGTACt       <  1:136006/141‑1 (MQ=255)
                                                                                                                         tACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGGCGCTTCTTCATGTGTTTTAAATACAGCAACTGGCTTTCCTGATTGTACt       <  1:427464/141‑1 (MQ=255)
                                                                                                                         tACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGGCGCTTCTTCATGTGTTTTAAATACAGCAACTGGCTTTCCTGATTGTACt       <  1:94634/141‑1 (MQ=255)
                                                                                                                         tACATGATTAAACCCTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGGCGCTTCTTCATGTGTTTTAAATACAGCAACTGGCTTTCCTGATTGTACt       <  2:347372/141‑1 (MQ=255)
                                                                                                                         tACATGATTAAACCCTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGGCGCTTCTTCATGTGTTTTAAATACAGCAACTGGCTTTCCTGATTGTACt       <  2:196373/141‑1 (MQ=255)
                                                                                                                              gATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGGCGCTTCTTCATGTGTTTTAAATACAGCAACTGGCTTTCCTGATTGTACTAACAt  >  2:293681/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                   
AACAGTGCCTGCTAAATCTCCATTAAAGTGTTGATTACCTAGCTCGGCAAATGTCTCATAAGTACCTTGCACAATACCTTGAGAACCAATATATATCCAACTACCGGCTGTCATCTGTCCATACATGATTAAACCTTTTTTATCTAATTCATTAAAATGATCCCAGTTTGCCCATTCAGGCACTAATACTGAATTTGAAATTAATACACGTGGCGCTTCTTCATGTGTTTTAAATACAGCAACTGGCTTTCCTGATTGTACTAACAT  >  CP000730/2450935‑2451201

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