Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_tpiA 871,104 G→T 8.0% Q46H (CAG→CAT)  gsiC → glutathione ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA871,1040GT 8.0% 136.7 / 2.4 51Q46H (CAG→CAT) gsiCglutathione ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (19/26);  new base T (2/2);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

GGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCAGGGGGATTTTGGCCTGTCGATGGTGTCGCGTCGTCCGGTTGCCGATGAGATTGCCAGCCGCTTTA  >  NC_000913_3_pae_tpiA/870958‑871252
                                                                                                                                                  |                                                                                                                                                    
gggAGTGGAATGCTTAATTACGTTATCAATCGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGt                                                                                                                                                    <  2:773258/149‑1 (MQ=255)
 agagtcaaaTGCTTAATTACGTTATCACACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGtt                                                                                                                                                   <  2:893929/142‑1 (MQ=255)
     tGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAg                                                                                                                                               >  2:37354/1‑149 (MQ=255)
      ggAATGCTTAGTTACGTGATCGAACGCTTACTGGGGTTGATTCCGACGAGGTTGAGCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGa                                                                                                                                              <  1:824906/149‑1 (MQ=255)
       gAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCCGGTTATAGaa                                                                                                                                             >  2:553455/1‑149 (MQ=255)
         aTGCTTAATTACGTGATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACt                                                                                                                                           <  1:786543/149‑1 (MQ=255)
             ttAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGAGATCAGGCGCGATTGATTGCCGGGCCGGAAGCTGGTGCGCAGGTTATAGAACTtgtg                                                                                                                                       >  2:392760/1‑149 (MQ=255)
                     gTTATCAAGCGCTTCCTGGCTTGGGTTCCGGCGCTTTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTcagca                                                                                                                               <  2:388786/149‑1 (MQ=255)
                     gTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTcagca                                                                                                                               >  1:855471/1‑149 (MQ=255)
                        aTCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCt                                                                                                                            <  1:173512/149‑1 (MQ=255)
                           aaaCGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCATGTTATAGAACTGGTGCGTCAGCAGCTggg                                                                                                                         >  1:355798/1‑149 (MQ=255)
                           aaaCGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCATGTTATAGAACTGGTGCGTCAGCAGCTggg                                                                                                                         >  1:355901/1‑149 (MQ=255)
                            gcgtcTTAGTGGGGTGTAGTCCGACGCTCTTTATCGTCGCGGTGCTGGTGTTTTTATTTGTCCTTATGGTGCCCGGCGATCCGGCGCGATTGATTGGCGGGCCCGAAGTTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTgggg                                                                                                                        <  2:549232/145‑1 (MQ=255)
                            aaCGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGATTATAGAACTGGTGCGTCGGCGGCTgggg                                                                                                                        >  1:108237/1‑149 (MQ=255)
                               gCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTg                                                                                                                     >  1:593874/1‑149 (MQ=255)
                                 ttACTGGGGTTGATTCCGACGCTGTTGATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGa                                                                                                                   <  1:343865/149‑1 (MQ=255)
                                   aCTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTTTCTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATc                                                                                                                 <  1:501791/149‑1 (MQ=255)
                                     tGGGGTTCATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAg                                                                                                               >  1:670524/1‑149 (MQ=255)
                                              ttCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTAt                                                                                                      <  1:463620/149‑1 (MQ=255)
                                                          ttATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGGGTAGTTGAAACCTCTTTTGTATTAGGAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTgg                                                                                          <  1:686290/149‑1 (MQ=255)
                                                           tATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGc                                                                                         <  2:643015/149‑1 (MQ=255)
                                                            aTCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCa                                                                                        <  2:332302/149‑1 (MQ=255)
                                                            aTCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCa                                                                                        <  2:332151/149‑1 (MQ=255)
                                                            aTCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCa                                                                                        >  1:132356/1‑149 (MQ=255)
                                                             tCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCAc                                                                                       >  1:454741/1‑149 (MQ=255)
                                                               gTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTGATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCAGCTGTATCACAAGTGCTGGCACta                                                                                     >  2:726177/1‑149 (MQ=255)
                                                                    ggTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCa                                                                                <  1:316202/149‑1 (MQ=255)
                                                                          ggTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATg                                                                          <  1:37354/149‑1 (MQ=255)
                                                                            tgtTTTTGTTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAAGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTAGAAGAAATTCTGTCTCGATATAAGCAAtccg                                                                        >  1:215970/1‑146 (MQ=255)
                                                                            tgtTTTTATTTGTCCATATGCGGCCCGGCGATCCGGCGCGATTGATTGGCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCt                                                                        <  2:606967/149‑1 (MQ=255)
                                                                              tttttATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCtgt                                                                      >  1:896518/1‑149 (MQ=255)
                                                                                tttATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGc                                                                    >  1:854798/1‑149 (MQ=255)
                                                                                 ttATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCa                                                                   <  1:933351/149‑1 (MQ=255)
                                                                                    tttGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCATGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCAggg                                                                <  2:355901/149‑1 (MQ=255)
                                                                                    tttGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCATGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCAggg                                                                <  2:355798/149‑1 (MQ=255)
                                                                                     ttGTCCCTATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCAgggg                                                               <  1:277992/149‑1 (MQ=255)
                                                                                          cctatGCTGCCCGGCGATCCCGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCAGGGGGAttt                                                          <  2:176436/147‑1 (MQ=255)
                                                                                          cagatGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCCCAACCTGTAGATTTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCAGGGGGAttt                                                          <  1:918780/146‑1 (MQ=255)
                                                                                             atGCTGCCCGGCCATCCGGCGCGATTGATTGCCGGCCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCAGGGGGATTTTgg                                                       <  2:275412/149‑1 (MQ=255)
                                                                                                   ccgggCGATCCGGCGCGATTGATTGCCGGGCACGAAGCTTATGCGCAGGTTATAGAACTGGTGCGTCAGTAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCAGGGGGATTTTGGCCTGTc                                                 <  2:814174/146‑1 (MQ=255)
                                                                                                   cccGGCGATCCGGCGCGATTGATTGCCGGGCCCGAGGATGTGGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCAGGAGGATCTTGGCCTGTc                                                 >  1:229029/1‑149 (MQ=255)
                                                                                                      ggCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCAGGGGGATTTTGGCCTGTCGAt                                              >  1:732593/1‑149 (MQ=255)
                                                                                                       gCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCAGGGGGATTTTGGCCTGTCGATg                                             <  1:392760/149‑1 (MQ=255)
                                                                                                              ggCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGAACTATATCAGCAATGCTGTGCAGGGGGATTTTGGCCTGTCGATGGTGTcgc                                      <  2:5526/149‑1 (MQ=255)
                                                                                                                      tgattgCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCAGGGGGATTTTGGCCTGTCGATGGTGTCGCGTCGTCCg                              <  2:130145/149‑1 (MQ=255)
                                                                                                                       gattgCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCAGGGGGATTTTGGCCTGTCGATGGTGTCGCGTCGTCCgg                             >  1:23984/1‑149 (MQ=255)
                                                                                                                            ccGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCAGGGGGATTTTGGCCTGTCGATGGTGTCGCGTCGTCCGGTTGcc                        >  1:841578/1‑149 (MQ=255)
                                                                                                                               ggCCCGAAGCGGAGGCGCAGGGTATAGGACGGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTTTGGAACTATATCAGCAATGCTGTGCAGGGGGATTTTGGCCTGTCGATGGTGTAGCGTCGTCCGGTTGCCGAt                     <  2:766155/149‑1 (MQ=255)
                                                                                                                                 cccGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCAGGGGGATTTTGGCCTGTCGATGGTGTCGCGTCGTCCGGTTGCCGATga                   >  2:481901/1‑149 (MQ=255)
                                                                                                                                          gATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCAGGGGGATTTTGGCCTGTCGATGGTGTCGCGTCGTCCGGTTGCCGATGAGATTGCCAg          >  2:675056/1‑149 (MQ=255)
                                                                                                                                               gcAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCAGGGGGATTTTGGCCTGTCGATGGTGTCGCGTCGTCCGGTTGCCGATGAGATTGCCAGCCGCt     >  1:542093/1‑149 (MQ=255)
                                                                                                                                                  ggTTATAGTCCTGGTGCGTCAGCAGCGTGGGTTGGATCAGCCGCTCTGGCACCATTTCTGGCACTATATCAGCAATGCTGTGCAGGGGGATTTTGGCCTGTCGATGGTGTCGCGTCGTCCGGTTGCCGATGAGATTGCCAGCCGCTTTa  <  2:241713/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACTATATCAGCAATGCTGTGCAGGGGGATTTTGGCCTGTCGATGGTGTCGCGTCGTCCGGTTGCCGATGAGATTGCCAGCCGCTTTA  >  NC_000913_3_pae_tpiA/870958‑871252

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