Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,869,922 A→T 17.6% N68K (AAT→AAA)  U712_19940 ← Deoxyribose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,869,9220AT17.6% 100.7 / 2.7 58N68K (AAT→AAA) U712_19940Deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (20/26);  new base T (6/5);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGATTGACACATACTGAAGCAAATTTGTATGTTTTCGCTTCTTCAAT  >  CP006881/3869786‑3870053
                                                                                                                                        |                                                                                                                                    
ccACTACACCGCGGATATCAGCTTCCACCCCATCGTCTTCCTTGTCTGTTAAAGCGACAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGATATGTCGTCGTTTGTTTCGCACGCGTTTGTTTCAGTTGTagt                                                                                                                                    >  6:143702/1‑137 (MQ=255)
         cgcgGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCGTTTGAAATGGCGTCTTTTGTTTCGAAAGCGTTTGTTTGAGTTTTTTTGGCACCGAg                                                                                                                           >  4:197146/1‑139 (MQ=255)
           cgGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATTGCGTGTGTTTTTTCGCACGCGTTTGTTTGATGTGTATTTGGACCGAGCg                                                                                                                         >  3:247794/1‑139 (MQ=255)
           cgGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTGTTTTGTTTCGAACGCGTTTGTTTCAGGTGTATTGGCACCGAGCg                                                                                                                         >  8:60933/1‑139 (MQ=255)
           cgGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCg                                                                                                                         <  7:36213/139‑1 (MQ=255)
           cgGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGAGGCAATAATAATGACCATATCCACTTCAGTGGCGCCTATTTAAATGGCGTGGTTTTTTTGGTCAGCCTTTGTTTAAGTTTTTATTGCACCAAGCg                                                                                                                         >  6:1409/1‑139 (MQ=255)
             gATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTATAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCggg                                                                                                                       <  4:149096/139‑1 (MQ=255)
             gATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCggg                                                                                                                       <  8:62647/139‑1 (MQ=255)
                atCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTAAGTTTTATTTGCACCGCGCGGGAAg                                                                                                                    >  8:105165/1‑139 (MQ=255)
                   aGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTTAAATGGCGGCTTTTTGTTGGGCCGCGTTTGTTTTAGTTTTATTGGCCCCCGGCGGGGCGCCg                                                                                                                 >  8:186001/1‑139 (MQ=255)
                        ccaccaCATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGa                                                                                                             <  7:105165/138‑1 (MQ=255)
                        ccaccaCATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTAGAACGCTTTTGTTTCAGTTGTATTTGCACCCAGCGGGAATCCGAttaa                                                                                                            >  8:45539/1‑136 (MQ=255)
                        ccaccaCATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCGTTTTAAATTGCGTGGTTTTTTTTGTACGCGTTTTTTTTGTTTTTTTTTGGACCCAGCCGGTAGCCCATGaa                                                                                                            >  4:203045/1‑138 (MQ=255)
                              caTCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGACGCCGATGACCGTACa                                                                                                      <  7:356693/139‑1 (MQ=255)
                               aTCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGGGTTTTTTTTTGCACGCTTTTTGTTAAGTTGTTGTTGCCCCCACCCCGCAGCCGCTGACCCTACa                                                                                                      >  7:75427/1‑138 (MQ=255)
                                   tCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTGAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑t                                                                                                >  8:106376/1‑139 (MQ=255)
                                   tCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACAATATCCACTTCAGTGGCGACATTTGAAATTGCGTGGTTTGTTTAGAACCCTTTTGTTAGATATGTAATGGCAACCAGCGGGAAGCAGATGACCCGACAAACga                                                                                                 >  2:116750/1‑138 (MQ=255)
                                    cTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCTGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑Tc                                                                                               <  5:321919/139‑1 (MQ=255)
                                      tCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTTTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATTAACGTACAAACG‑TCGa                                                                                             >  8:356693/1‑139 (MQ=255)
                                       ccTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑Tc                                                                                               <  7:45539/136‑1 (MQ=255)
                                            tCTTTTAAAGCGGCAAATTTAATGACCAAATCCACTGCTCTGGCGCCTTTTGAAATGGCGTATTTTGTTTCTAACGCTTTGGTTTCTGGTGTAGTGGCCCCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAg                                                                                       <  5:120685/139‑1 (MQ=255)
                                            gcTTTTAAAGCGGAAATATTAATGACCATATCCACTTCAGGGGCGCCTTTTTAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCTGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAg                                                                                       <  8:223647/138‑1 (MQ=255)
                                             cTTTTAAAGCGGAAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGGATTTTTTTTCGAACGCTTTTGTTTCAGTTGTTTTGCCACCGAGCGGGAAGCCTATGACCGTACAAACG‑TCGACTCCAg                                                                                       <  5:1409/138‑1 (MQ=255)
                                                 tAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑Tc                                                                                               <  7:22358/126‑1 (MQ=255)
                                                 tAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑Tc                                                                                               >  8:22358/1‑126 (MQ=255)
                                                 tAAAGCGGCAATATTAATGACCATATCAACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGTAAGATCTTGTTTAAGTTGTATTGTCACCCCGCCGGATTCCGAATTACGGACAAAAG‑TCGACACCAGTTccc                                                                                  >  2:44466/1‑139 (MQ=255)
                                                      cGGCAATATTAATGACCATATCCACTTCCGTGGCGCCTTTTGAAATGGCGGCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGCGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAg                                                                             <  3:271476/139‑1 (MQ=255)
                                                        gCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGGTCCCGTAAGcc                                                                           >  7:223647/1‑138 (MQ=255)
                                                                 attcccATATACAATTCAGTGGCGCCTTTGGATAGGGAGTCTTTTGTAGGTACCGCTTTTGTTTCTGTTGTCTTGCCCCCGCGCGGGACGCCCTTGACCGTACAACCG‑GCGACTCCAGTTCCCTTAAGCTCTTTTgcag                                                                  <  2:177560/135‑1 (MQ=255)
                                                                      cATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGc                                                             >  7:366546/1‑139 (MQ=255)
                                                                      cATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGc                                                             >  8:205465/1‑139 (MQ=255)
                                                                         atCCACTTCAGTGGCGCCTTTTGAACTGGCGTGTTTTGTTTCGAACGCTTTTGTTTCTGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTc                                                           <  1:116750/138‑1 (MQ=255)
                                                                         atCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGt                                                                                      <  8:121152/111‑1 (MQ=255)
                                                                         atCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGt                                                                                      >  7:121152/1‑111 (MQ=255)
                                                                               ttCAGTGGCGCCTTTTGAAATGGCGTCGTTTGTTTGGCACGCTTTTGTTTGAGTTGTTTTGGCACCGCGCGGGGCTCCGATGACCGTACAAACC‑GAGCCACCAGTTCCCGTAAGCGCTTTTTCAGCCAGCACCCCCCaa                                                    >  8:176265/1‑138 (MQ=255)
                                                                                tCAGTGGCGCCTTTTGAAATGGCGTGTTTTGTTTCGAAAGCGTTTGTTTGAGGTGTTTTGGCACCGAGCGGGGATCCGCAGACCGGTCAAACG‑TCGGCTGCAGTGGCCGTCAGCGCTGTTGCTGCCAGGGCCCCCCagg                                                   >  4:174142/1‑137 (MQ=255)
                                                                                  aGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTc                                                 >  1:218076/1‑139 (MQ=255)
                                                                                  aGTGGCGCCTTTTGAAATGGCGTCATTTGTTTCGAACGCATTTGTTTCAGTTGTATTGGCACCGAGCGGCAAGCCGAAGACCGTACAAACG‑TCGCCACCAGTTCCCATAAGCTCTTTTGCAGCAACCTCCACCCATGTc                                                 >  7:44426/1‑139 (MQ=255)
                                                                                    tGGCGCTTTTTTAATTGGGGTTTTTTTTTTGGAAGGCTTTTGTTTCTGTTGTTTTGGCACGTCGGGGTACGCCTATGACCGTAAAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGt                                                  <  6:49052/136‑1 (MQ=255)
                                                                                    tGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTGAGTTGTATTTGCACCGCGCGGGTAGCCGATGACCGTACAAACG‑TCGACTCCAGTTGCCGTTAGCTCTTTTGCAGCAAGCAGCCCCCATGTCg                                                >  4:338913/1‑138 (MQ=255)
                                                                                     ggCGCCTTTTTAAATGGGGTTTTTTGTTTCCAACGCTTTTGTTTCTGTTGTTTTGGCACCGAGCGGTAAGCCGATGACCTTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGa                                              <  1:245802/139‑1 (MQ=255)
                                                                                     ggCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGa                                              <  7:205465/139‑1 (MQ=255)
                                                                                     ggCGCATTTTGAAATGGCGTGTTTTGTTTCGAACGCTTTTGTTTGAGTTGTTTTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGa                                              <  8:75427/139‑1 (MQ=255)
                                                                                      gcgcCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCCCCGAGCGGGAAGCCGATGACCGTACAAACGCTCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGa                                              <  2:33133/139‑1 (MQ=255)
                                                                                                  aTGGCGTATTTTGTTTCCACCGCTTTTGTTGCAGTTGTATTGGCACCGTGCGGTAAGCCGATAACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGATTGACACATACTg                                 <  1:44466/139‑1 (MQ=255)
                                                                                                               tttAGAACGCTTTTGTTTCAGTTTTATTTGCACCGAGCGGGAAGCCGATGACCGTACAAAAG‑TCGACTCCAGTTCCCTTAATCTCTTTTTCAGCAAGCTCCACCCAAGTAGGATTTACACATAAAGAATCAAAATTGTa                    >  3:213685/1‑139 (MQ=255)
                                                                                                                    aaCGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGTAAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCAc                                                        <  2:203291/98‑1 (MQ=255)
                                                                                                                    aaCGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCAc                                                        >  1:203291/1‑98 (MQ=255)
                                                                                                                    aaCGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGATTGACACATACTGAAGCAAATTTGTATGttt               >  1:177479/1‑139 (MQ=255)
                                                                                                                            tGTTTCAGTTGTATTGGCACCGAGCGGGTAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGATTGACACATACTGAAGCAAATTTGTATGTTTTCGcttc        <  4:247794/138‑1 (MQ=255)
                                                                                                                                 ctgTTGTATTGGCACGGAGCGGGAAGCCGATGACCTTAAAAACC‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGATTGACACATACTGAAGCAAATTTGTATGTTTTCGCTTCTTCAAt  <  2:368376/137‑1 (MQ=255)
                                                                                                                                  aGTTTTATTGGCACGGAGGGGGAAGCCGATGACCGTACATTCG‑TCGACTCCAGTTCACTTAAGCTCTTTTGCAGCAAGCTACACCCATGTCGGATTGACACATACTGAAGCAAATTTGTATGTTTTCGCTTCTTCAAt  <  1:282454/138‑1 (MQ=255)
                                                                                                                                  aGTTGTTTTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGATTGACACATACTGAAGCAAATTTGTATGTTTTCGCTTCTTCAAt  <  3:203045/138‑1 (MQ=255)
                                                                                                                                  aGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGATTGACACATACTGAAGCAAATTTGTATGTTTTCGCTTCTTCAAt  <  7:186001/138‑1 (MQ=255)
                                                                                                                                  aGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGATTGACACATACTGAAGCAAATTTGTATGTTTTCGCTTCTTCAAt  <  7:106376/138‑1 (MQ=255)
                                                                                                                                  aGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGATTGACACATACTGAAGCAAATTTGTATGTTTTCGCTTCTTCAAt  <  3:197146/138‑1 (MQ=255)
                                                                                                                                  aGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGATTGACACATACTGAAGCAAATTTGTATGTTTTCGCTTCTTCAAt  <  2:177479/138‑1 (MQ=255)
                                                                                                                                  aGTTGTATTGGCACCGAGCGGGAACCCGATGACCGTACAACGG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGATTGACACATACTGAAGCAAATTTGTATGTTTTCGCTTCTTCAAt  <  5:172562/138‑1 (MQ=255)
                                                                                                                                  aGTTGTATGGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGATTGACACATACTGAAGCAAATTTGTATGTTTTCGCTTCTTCAAt  <  4:282960/138‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                    
CCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACG‑TCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGATTGACACATACTGAAGCAAATTTGTATGTTTTCGCTTCTTCAAT  >  CP006881/3869786‑3870053

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