Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 869,724 A→C 9.6% N105T (AAC→ACC)  gsiB → glutathione ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913869,7240AC 9.6% 163.1 / 8.5 51N105T (AAC→ACC) gsiBglutathione ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (26/20);  new base C (3/2);  total (29/22)
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.06e-01

ATCGTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTT  >  NC_000913/869578‑869872
                                                                                                                                                  |                                                                                                                                                    
aTCGTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAAcg                                                                                                                                                    >  2:176486/1‑149 (MQ=255)
          cAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCgcggcgtt                                                                                                                                          >  4:532009/1‑147 (MQ=255)
              ggCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCGCGGCGTTGaaa                                                                                                                                      <  3:445423/149‑1 (MQ=255)
              ggCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGaaa                                                                                                                                      <  1:413238/149‑1 (MQ=255)
                 tGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCGCGGCGTTGAAATCg                                                                                                                                   <  3:532009/149‑1 (MQ=255)
                       gTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTg                                                                                                                             <  2:504610/149‑1 (MQ=255)
                             aTAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCgg                                                                                                                       >  1:302808/1‑149 (MQ=255)
                                aaGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGcc                                                                                                                    >  2:469939/1‑149 (MQ=255)
                                   agaTGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGc                                                                                                                     >  4:589361/1‑145 (MQ=255)
                                   agaTGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGc                                                                                                                     <  3:589361/145‑1 (MQ=255)
                                    gaTGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCg                                                                                                                >  3:247469/1‑149 (MQ=255)
                                      tGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGAt                                                                                                              >  3:320721/1‑149 (MQ=255)
                                        aaaCTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATcc                                                                                                            <  2:302808/149‑1 (MQ=255)
                                         aaCTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTGAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCg                                                                                                           >  4:368320/1‑149 (MQ=255)
                                          aCTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCgg                                                                                                          <  1:605945/149‑1 (MQ=255)
                                            tGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCg                                                                                                        >  3:62693/1‑149 (MQ=255)
                                                     tGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCtt                                                                                               <  1:34442/149‑1 (MQ=255)
                                                        tGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTaaa                                                                                            <  4:247469/149‑1 (MQ=255)
                                                            ggAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCt                                                                                        >  1:585281/1‑149 (MQ=255)
                                                                  ttATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCt                                                                                                                              >  2:416144/1‑105 (MQ=255)
                                                                  ttATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCt                                                                                                                              <  1:416144/105‑1 (MQ=255)
                                                                  ttATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAAcc                                                                                  >  1:86320/1‑149 (MQ=255)
                                                                  ttATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGACATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAAcc                                                                                  >  4:372317/1‑149 (MQ=255)
                                                                     taCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGt                                                                               >  1:121152/1‑149 (MQ=255)
                                                                      aCCGTTTCCGATGACGGCATTACTTACACCGTGACATTGCGGGAAGACATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGta                                                                              <  3:372317/149‑1 (MQ=255)
                                                                       ccGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGtat                                                                             <  4:569621/149‑1 (MQ=255)
                                                                        cGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGtata                                                                            >  1:71318/1‑149 (MQ=255)
                                                                               gATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAatat                                                                     <  3:129638/149‑1 (MQ=255)
                                                                                aTGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAACTTCCAGGATGGCACCGATTTCACCGCTGCGGCGGTGAAAGCGAATATGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATAtt                                                                    >  2:365353/1‑149 (MQ=255)
                                                                                  gACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGc                                                                  <  2:124521/149‑1 (MQ=255)
                                                                                  gACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGc                                                                  <  3:368320/149‑1 (MQ=255)
                                                                                       cATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTaaaa                                                             <  4:62693/149‑1 (MQ=255)
                                                                                       cATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTaaaa                                                             <  3:432406/149‑1 (MQ=255)
                                                                                          tacttacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACgg                                                          <  1:469939/149‑1 (MQ=255)
                                                                                           acttacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGa                                                         >  4:323767/1‑149 (MQ=255)
                                                                                             ttacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAg                                                       >  3:273627/1‑149 (MQ=255)
                                                                                              tacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGc                                                      <  2:585281/149‑1 (MQ=255)
                                                                                                  ccGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGcgatc                                                  <  4:320721/149‑1 (MQ=255)
                                                                                                    gTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGcgatcga                                                >  4:447536/1‑149 (MQ=255)
                                                                                                                 aaGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTa                                   >  2:620133/1‑149 (MQ=255)
                                                                                                                 aaGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTa                                   >  3:560623/1‑149 (MQ=255)
                                                                                                                    gCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAg                                >  3:162632/1‑149 (MQ=255)
                                                                                                                      aTTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGAt                              >  1:325443/1‑149 (MQ=255)
                                                                                                                        tAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTa                            <  4:274927/149‑1 (MQ=255)
                                                                                                                          aaTTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTAcc                          >  2:29456/1‑149 (MQ=255)
                                                                                                                                     ggCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGcc               <  2:86320/149‑1 (MQ=255)
                                                                                                                                      gCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCg              >  1:606515/1‑149 (MQ=255)
                                                                                                                                        aCCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGtt            >  3:70352/1‑149 (MQ=255)
                                                                                                                                        aCCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGtt            >  1:2184/1‑149 (MQ=255)
                                                                                                                                          cGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACACTCAAACAGCCGTtct          >  3:529671/1‑149 (MQ=255)
                                                                                                                                             tttCACCGCCGCGGCGTTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAg       >  1:457212/1‑149 (MQ=255)
                                                                                                                                                 aaCGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGtt   >  2:39264/1‑149 (MQ=255)
                                                                                                                                                 aaCGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGtt   <  1:356008/149‑1 (MQ=255)
                                                                                                                                                 aaCGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGtt   <  1:176486/149‑1 (MQ=255)
                                                                                                                                                  aCGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGttt  <  1:29456/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATCGTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTT  >  NC_000913/869578‑869872

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