Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 867,718 G→C 8.0% E67Q (GAA→CAA)  gsiA → glutathione ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913867,7180GC 8.0% 154.3 / 5.8 50E67Q (GAA→CAA) gsiAglutathione ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (20/26);  new base C (2/2);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

TGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAG  >  NC_000913/867572‑867855
                                                                                                                                                   |                                                                                                                                         
tGAATATTGCCTTTATGCAGGACCAGCAGAAGATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGaa                                                                                                                                         <  1:303858/149‑1 (MQ=255)
tGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGaa                                                                                                                                         <  2:801961/149‑1 (MQ=255)
tGAATATTGCCTTATTGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGTCGCTGGGAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGCCTGCGTTGGCATTGATGCGCCTGTTGGaa                                                                                                                                         <  1:303714/149‑1 (MQ=255)
tGAATATCGCCTTTATGCAGGACCAGCAGAAAACAGCTGCGGTCCGCAATCTCGCTTGTAGTCTG‑CAACGCGGTGCGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGaa                                                                                                                                         <  2:431022/149‑1 (MQ=255)
aggttatGGCCTTTTTGCAGACCCAGCCGAAAAAAGCTGCGGTCCGCATATGCGCTCTTAG‑CTGCCAACACGGTGATACGCAGGCACTCGTTGGCGAATCCGGCTCCGGTAAGTCAGTGCCTGCGTTGGCATTGATGCGCCTGTTGGaa                                                                                                                                         <  2:431048/145‑1 (MQ=255)
  aagatTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACa                                                                                                                                       <  2:473241/146‑1 (MQ=255)
  aaTATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTAAACCATAAATAGATTAGTCGG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACa                                                                                                                                       <  2:564791/149‑1 (MQ=255)
    tatTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAgg                                                                                                                                     >  1:628805/1‑149 (MQ=255)
        gCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATCGATGCGCCTGTTCCAACAggcgca                                                                                                                                 <  1:403320/149‑3 (MQ=255)
         ccTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAggcgggc                                                                                                                                >  2:183307/1‑149 (MQ=255)
         ccTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAggcgggc                                                                                                                                >  2:467612/1‑149 (MQ=255)
           tttATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTCCAACAggcgcatca                                                                                                                              <  1:301090/149‑6 (MQ=255)
                 cAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTCCGTCAGTGACTGCGGGTACATCGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGt                                                                                                                        >  2:599071/1‑149 (MQ=255)
                  aGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCCCCTGTTGGAACAGGCGGGGGGTTTAGTa                                                                                                                       >  1:84018/1‑149 (MQ=255)
                    gACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCCGTTGGAACAGGCGGGCGGTTTAGTGCa                                                                                                                     >  2:38023/1‑149 (MQ=255)
                       cagcagAAAATAGCTGCGGTCCGCAATCTCTCTTTTTTCCTG‑CAACCCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTg                                                                                                                  <  1:346561/149‑1 (MQ=255)
                       cagcagAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTg                                                                                                                  >  1:42936/1‑149 (MQ=255)
                             aaaaTAGCTCCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATaa                                                                                                            <  1:87897/149‑1 (MQ=255)
                                   gCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCt                                                                                                      >  2:738306/1‑149 (MQ=255)
                                          tCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAg                                                                                               >  2:699932/1‑149 (MQ=255)
                                            cGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAgcg                                                                                             >  1:104226/1‑149 (MQ=255)
                                            cGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAAGGCTGTTGCAgcg                                                                                             >  2:306584/1‑149 (MQ=255)
                                            cGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGGGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATCAAATGCTGTTGCAgcg                                                                                             >  2:702781/1‑149 (MQ=255)
                                               aaTCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAgcggcg                                                                                          <  2:628805/149‑1 (MQ=255)
                                                   tctctTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGt                                                                                      >  1:415009/1‑149 (MQ=255)
                                                        tttAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGa                                                                                 <  1:470207/149‑1 (MQ=255)
                                                              tcgtcAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGa                                                                            <  2:455589/145‑1 (MQ=255)
                                                              cTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGAt                                                                           <  2:205992/149‑1 (MQ=255)
                                                                   aaCGCGGTGGGTCGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAATTGCTTTTGCAGCGGCGCAGTCGCGAAGTGATTGaa                                                                       <  1:194774/149‑1 (MQ=255)
                                                                   aaCGCGGTGAGACGCTGGCAATTGTTGGCGCATCCCGCTCCTGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGaa                                                                       <  1:149557/149‑1 (MQ=255)
                                                                     cgcgGTGAGACCCTGGCAATTGCTGGCGAATCCGGCTCCGGTAAGTCAGTGACCGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACt                                                                     <  2:405028/149‑1 (MQ=255)
                                                                              aCGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCATTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCa                                                            <  1:170628/149‑1 (MQ=255)
                                                                                       aTTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTAGTTGAACTTAGCGAGCCGGACGCGGc                                                   >  2:475211/1‑149 (MQ=255)
                                                                                        ttgttgGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATGGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGca                                                  <  2:272406/149‑1 (MQ=255)
                                                                                            tgGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGGGCAGACAGCTGCACaaa                                              >  1:71422/1‑149 (MQ=255)
                                                                                             ggCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAAt                                             <  1:180322/149‑1 (MQ=255)
                                                                                                    ccGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCAt                                      >  1:474203/1‑149 (MQ=255)
                                                                                                     cGGCTCCGGTAAGTAAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATg                                     >  2:279869/1‑149 (MQ=255)
                                                                                                                  tCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCggg                        >  2:624906/1‑148 (MQ=255)
                                                                                                                   cAGTGACTGCGTTTGTATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTTCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCTCCATGTTCGCGGTGCGGat                       <  1:581482/149‑1 (MQ=255)
                                                                                                                   cAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGat                       <  1:756433/149‑1 (MQ=255)
                                                                                                                     gTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGAGAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGatat                     >  1:531712/1‑149 (MQ=255)
                                                                                                                      tGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATg                    <  1:75686/149‑1 (MQ=255)
                                                                                                                           gCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCgat               <  1:602086/149‑1 (MQ=255)
                                                                                                                               tGGCATTGATGCGCCTGTTGGAACAGGCGGGAGGTTTAGTAAATGGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGata           <  1:543351/149‑1 (MQ=255)
                                                                                                                                    ttGATGCGCCTGTGGGAACAGGCGGGCGGTTTCGTACAGTGGGATACAATGCTGGTGCAGCGGCGCAGTCGCGAAGTGATTGACCTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCa      <  2:365095/149‑1 (MQ=255)
                                                                                                                                    accgcccgcCTGTTCCAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGCCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAAAGGGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCa      >  2:85598/7‑149 (MQ=255)
                                                                                                                                       aTGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGa   <  1:279869/149‑1 (MQ=255)
                                                                                                                                       aTGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGa   >  1:681391/1‑149 (MQ=255)
                                                                                                                                        cccgcCTGTTCCAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAAAGCGCCATGTTCGCGGTGCGGATATGGCGATGAAACTTCAGGac  >  2:276297/3‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
TGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAG  >  NC_000913/867572‑867855

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