Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 872,801 G→A *304* (TAG→TAA gsiD → glutathione ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913872,8010GA98.6% 225.3 / NA 71*304* (TAG→TAAgsiDglutathione ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (40/30);  minor base T (1/0);  total (41/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTT  >  NC_000913/872669‑872932
                                                                                                                                    |                                                                                                                                   
ggCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGt                                                                                                                               >  1:442976/1‑139 (MQ=255)
  cTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGttt                                                                                                                             >  2:1744840/1‑139 (MQ=255)
  cTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGttt                                                                                                                             >  2:1174849/1‑139 (MQ=255)
  cTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGttt                                                                                                                             >  2:1319386/1‑139 (MQ=255)
     gAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGaa                                                                                                                          >  2:112821/1‑139 (MQ=255)
         ggATATGGTTATCGAGCCGCATGTCGCTGTGTTTCCGGCCCTTGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATAtt                                                                                                                      <  1:1570841/139‑1 (MQ=255)
           atatGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGc                                                                                                                    >  2:353083/1‑139 (MQ=255)
              tGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTg                                                                                                                 <  1:353083/139‑1 (MQ=255)
                   aTCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAgg                                                                                                            >  2:158164/1‑139 (MQ=255)
                   aTCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAgg                                                                                                            >  2:1542798/1‑139 (MQ=255)
                    tCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTa                                                                                                                                  <  1:1466144/116‑1 (MQ=255)
                    tCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAggg                                                                                                           >  2:1570841/1‑139 (MQ=255)
                     cgggccggATGTCGCTGTTTTTCCGGCCCTGGCTATTTTGCTTGCCGTACTAGCGGTCCATTTGTTTGTCTGAGCTTTACGCCAGGGACCGGAATCGCAAATTAGAGGATATTTACGATTGAATATTGCTTGAAAGGGt                                                                                                          >  1:1334459/4‑139 (MQ=255)
                         ccgcATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGa                                                                                                                >  1:182415/1‑129 (MQ=255)
                         ccgcATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGa                                                                                                                <  2:182415/129‑1 (MQ=255)
                         ccgcATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATc                                                                                                      >  2:1531424/1‑139 (MQ=255)
                           gcATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCAc                                                                                                    <  2:1515446/139‑1 (MQ=255)
                           gcATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGATGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCAc                                                                                                    <  1:158164/139‑1 (MQ=255)
                               gTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTca                                                                                                <  2:1521239/139‑1 (MQ=255)
                                       ttttCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAat                                                                                        >  2:1019706/1‑139 (MQ=255)
                                       ttttCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAat                                                                                        >  1:1579113/1‑139 (MQ=255)
                                         ttCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAatta                                                                                      <  1:1174849/139‑1 (MQ=255)
                                          tCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAattat                                                                                     >  2:1684329/1‑139 (MQ=255)
                                               cccTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGccc                                                                                <  1:112821/139‑1 (MQ=255)
                                                     cTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCa                                                                          <  1:1684329/139‑1 (MQ=255)
                                                              tGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaaa                                                                                                              >  1:863815/1‑94 (MQ=255)
                                                              tGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaaa                                                                                                              <  2:863815/94‑1 (MQ=255)
                                                                  cGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGccc                                                                                >  1:763546/1‑120 (MQ=255)
                                                                  cGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGccc                                                                                <  2:763546/120‑1 (MQ=255)
                                                                    tACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAAtt                                                                                                                                   <  1:387657/67‑1 (MQ=255)
                                                                    tACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAAtt                                                                                                                                   >  2:387657/1‑67 (MQ=255)
                                                                     aCTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTct                                                          >  1:317743/1‑139 (MQ=255)
                                                                            ttCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTa                                                                                                                                  <  1:270403/60‑1 (MQ=255)
                                                                            ttCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTa                                                                                                                                  >  2:1832703/1‑60 (MQ=255)
                                                                            ttCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTa                                                                                                                                  <  1:1832703/60‑1 (MQ=255)
                                                                            ttCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTa                                                                                                                                  >  2:270403/1‑60 (MQ=255)
                                                                              cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGt                                                                                                                               <  1:1324918/61‑1 (MQ=255)
                                                                              cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGt                                                                                                                               >  2:1324918/1‑61 (MQ=255)
                                                                              cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGa                                                                                           <  2:962747/97‑1 (MQ=255)
                                                                              cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGa                                                                                           >  1:962747/1‑97 (MQ=255)
                                                                              cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGcc                                                                                 >  1:731784/1‑107 (MQ=255)
                                                                              cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGcc                                                                                 <  2:731784/107‑1 (MQ=255)
                                                                              cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTctgc                                                        >  2:257787/1‑132 (MQ=255)
                                                                              cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTctgc                                                        <  1:257787/132‑1 (MQ=255)
                                                                              cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGAtttt                                                 >  2:446178/1‑139 (MQ=255)
                                                                                   tgttgGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGa                                            >  2:837156/1‑139 (MQ=255)
                                                                                   tgttgGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGa                                            <  1:837156/139‑1 (MQ=255)
                                                                                          cGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCaa                                                                         <  1:907674/103‑1 (MQ=255)
                                                                                          cGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCaa                                                                         >  2:907674/1‑103 (MQ=255)
                                                                                             tGGTTTACGTGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGaa                                                                                                                          >  2:1484913/1‑51 (MQ=255)
                                                                                             tGGTTTACGTGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGaa                                                                                                                          <  1:1484913/51‑1 (MQ=255)
                                                                                             tGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAAt                                  <  1:1019706/139‑1 (MQ=255)
                                                                                                  tACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGt                             >  2:1706882/1‑139 (MQ=255)
                                                                                                         tGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaa                                                                                                               >  2:1744687/1‑50 (MQ=255)
                                                                                                         tGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaa                                                                                                               <  1:1744687/50‑1 (MQ=255)
                                                                                                         tGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGccc                                                                                <  1:692137/81‑1 (MQ=255)
                                                                                                         tGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGccc                                                                                >  2:692137/1‑81 (MQ=255)
                                                                                                         tGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAg                      >  2:1057376/1‑139 (MQ=255)
                                                                                                                   ccGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTaaaa                        <  1:1396006/127‑1 (MQ=255)
                                                                                                                   ccGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTaaaa                        >  2:1396006/1‑127 (MQ=255)
                                                                                                                     gAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTct                                                          >  2:196242/1‑91 (MQ=255)
                                                                                                                     gAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTct                                                          <  1:196242/91‑1 (MQ=255)
                                                                                                                      aaaaTTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGa         >  2:1552397/1‑139 (MQ=255)
                                                                                                                      aaaaTTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACACGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGa         <  1:1552397/139‑1 (MQ=255)
                                                                                                                        aaTTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACa       >  1:925175/1‑139 (MQ=255)
                                                                                                                           tAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGCTTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAAtt    >  1:1386721/1‑139 (MQ=255)
                                                                                                                             aaGGATAATTACGTTTGAGTATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTg                                                                                   >  2:204410/1‑58 (MQ=255)
                                                                                                                             aaGGATAATTACGTTTGAGTATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTg                                                                                   <  1:204410/58‑1 (MQ=255)
                                                                                                                             aaGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCa                                                                          >  2:287533/1‑67 (MQ=255)
                                                                                                                             aaGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCa                                                                          <  1:287533/67‑1 (MQ=255)
                                                                                                                             aaGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAAtttt  >  2:1088082/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                   
GGCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTT  >  NC_000913/872669‑872932

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