Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 869,338 C→T P607S (CCG→TCG)  gsiA → glutathione ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913869,3380CT100.0% 103.1 / NA 33P607S (CCG→TCG) gsiAglutathione ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/15);  total (18/15)

TGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTG  >  NC_000913/869210‑869467
                                                                                                                                |                                                                                                                                 
tGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTa                                                                                                                         <  1:369552/139‑1 (MQ=255)
     gcAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCg                                                                                                                    <  1:193721/139‑1 (MQ=255)
                         gTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGca                                                                                                                  >  2:138076/1‑121 (MQ=255)
                         gTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGca                                                                                                                  <  1:138076/121‑1 (MQ=255)
                                                     ctgTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCg                                                                    >  2:554265/1‑139 (MQ=255)
                                                      tgTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCgg                                                                   >  2:638811/1‑139 (MQ=255)
                                                               aTCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACca                                                                                                            >  1:265510/1‑89 (MQ=255)
                                                               aTCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACca                                                                                                            <  2:265510/89‑1 (MQ=255)
                                                                   tCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGc                                                      <  1:554265/139‑1 (MQ=255)
                                                                     ccAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCaa                                                    <  1:70501/139‑1 (MQ=255)
                                                                        gCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAagag                                                 <  2:150306/139‑1 (MQ=255)
                                                                          aaTATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCt                                               >  1:320054/1‑139 (MQ=255)
                                                                            tatTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGt                                             >  2:764805/1‑139 (MQ=255)
                                                                                   cTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGt                                      >  2:337236/1‑139 (MQ=255)
                                                                                      cGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAgg                                                                      >  1:564183/1‑104 (MQ=255)
                                                                                      cGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAgg                                                                      <  2:564183/104‑1 (MQ=255)
                                                                                        tGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTaca                                          <  2:64869/130‑1 (MQ=255)
                                                                                        tGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTaca                                          >  1:64869/1‑130 (MQ=255)
                                                                                        tGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTgg                                 >  2:99241/1‑139 (MQ=255)
                                                                                           cGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAAc                                                                                                              >  2:354468/1‑59 (MQ=255)
                                                                                           cGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAAc                                                                                                              <  1:354468/59‑1 (MQ=255)
                                                                                                     gCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCaa                                                    >  1:715836/1‑107 (MQ=255)
                                                                                                     gCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCaa                                                    <  2:715836/107‑1 (MQ=255)
                                                                                                       aGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGAACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTggg                  >  2:687383/1‑139 (MQ=255)
                                                                                                            tctcGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAgg                                                                      >  1:666451/1‑82 (MQ=255)
                                                                                                            tctcGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAgg                                                                      <  2:666451/82‑1 (MQ=255)
                                                                                                               cGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGa          <  1:99241/139‑1 (MQ=255)
                                                                                                                  tGCAATGCGTCGGGTTGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAg       >  2:518140/1‑139 (MQ=255)
                                                                                                                      aTGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTg                                              >  1:743903/1‑96 (MQ=255)
                                                                                                                      aTGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTg                                              <  2:743903/96‑1 (MQ=255)
                                                                                                                       tGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTg  >  1:413222/1‑139 (MQ=255)
                                                                                                                            cGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAagag                                                 >  2:300094/1‑87 (MQ=255)
                                                                                                                            cGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAagag                                                 <  1:300094/87‑1 (MQ=255)
                                                                                                                                |                                                                                                                                 
TGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTG  >  NC_000913/869210‑869467

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