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,8010GA100.0% 91.7 / NA 29*304* (TAG→TAAgsiDglutathione ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/14);  total (15/14)

TGAGGCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAG  >  NC_000913/872666‑872912
                                                                                                                                       |                                                                                                               
tGAGGCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTa                                                                                                              <  1:32473/139‑1 (MQ=255)
tGAGGCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTa                                                                                                              <  1:956527/139‑1 (MQ=255)
tGAGGCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTa                                                                                                              <  1:484182/139‑1 (MQ=255)
   ggCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGt                                                                                                           <  1:509142/139‑1 (MQ=255)
    gCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGtt                                                                                                          <  2:279484/139‑1 (MQ=255)
     cTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCGGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGttt                                                                                                         >  2:1026302/1‑139 (MQ=255)
           cGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAatat                                                                                                   >  1:571573/1‑139 (MQ=255)
                      aTCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAgg                                                                                        <  2:571573/139‑1 (MQ=255)
                       tCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAggg                                                                                       >  2:349728/1‑139 (MQ=255)
                        cgcgCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGt                                                                                      <  1:349728/139‑1 (MQ=255)
                          cgccgcATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTaa                                                                                    >  1:809187/1‑139 (MQ=255)
                            ccgcATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAata                                                                                                    >  2:356267/1‑121 (MQ=255)
                            ccgcATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAata                                                                                                    <  1:356267/121‑1 (MQ=255)
                                           tttCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAatt                                                                   >  2:783016/1‑139 (MQ=255)
                                              ccGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaaa                                                                                          <  2:268239/113‑1 (MQ=255)
                                              ccGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaaa                                                                                          >  1:268239/1‑113 (MQ=255)
                                                                               ttCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGAtt                               >  2:481470/1‑139 (MQ=255)
                                                                                 cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGa                                                                       >  1:380034/1‑97 (MQ=255)
                                                                                 cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGa                                                                       <  2:380034/97‑1 (MQ=255)
                                                                                      tgttgGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGa                        <  1:783016/139‑1 (MQ=255)
                                                                                            gCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGc                                                              <  1:911447/95‑1 (MQ=255)
                                                                                            gCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGc                                                              >  2:911447/1‑95 (MQ=255)
                                                                                            gCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCggg                  >  2:918701/1‑139 (MQ=255)
                                                                                                  gTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaa                                                                                           <  1:767684/60‑1 (MQ=255)
                                                                                                  gTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaa                                                                                           >  2:767684/1‑60 (MQ=255)
                                                                                                            tGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAg  >  2:106528/1‑139 (MQ=255)
                                                                                                            tGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAg  >  2:253002/1‑139 (MQ=255)
                                                                                                                 tGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAg                                                    >  2:586347/1‑84 (MQ=255)
                                                                                                                 tGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAg                                                    <  1:586347/84‑1 (MQ=255)
                                                                                                                                       |                                                                                                               
TGAGGCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAG  >  NC_000913/872666‑872912

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: