Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 867,717 G→C 8.0% L66F (TTG→TTC)  gsiA → glutathione ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913867,7170GC 8.0% 149.0 / 5.8 50L66F (TTG→TTC) 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 = 5.06e-01

ATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAG  >  NC_000913/867569‑867855
                                                                                                                                                     |                                                                                                                                          
aTCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTg                                                                                                                                            >  1:297374/1‑149 (MQ=255)
aTCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTg                                                                                                                                            >  2:602086/1‑149 (MQ=255)
   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)
                                                                                                                                                     |                                                                                                                                          
ATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAG  >  NC_000913/867569‑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: ‑