Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 868,422 (C)8→9 coding (903/1872 nt) gsiA → glutathione ABC transporter ATP binding subunit GsiA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913868,4141.C90.5% 131.8 / 6.9 42coding (895/1872 nt)gsiAglutathione ABC transporter ATP binding subunit GsiA
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base C (18/20);  total (20/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.40e-01

AACAGATTTTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑CCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGT  >  NC_000913/868280‑868547
                                                                                                                                       |                                                                                                                                        
aaCAGATTTTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccccc                                                                                                                                   <  2:19612/140‑1 (MQ=255)
aaCAGATTTTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccccc                                                                                                                                   >  1:670471/1‑140 (MQ=255)
         ttCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccccccccccctcg                                                                                                                         <  2:63167/140‑7 (MQ=255)
                 ccGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGaaaa                                                                                                                  <  1:510266/140‑1 (MQ=255)
                                  ccGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGaa                                                                                                 <  2:815315/140‑1 (MQ=255)
                                  ccGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGaa                                                                                                 <  1:679273/140‑1 (MQ=255)
                                   cGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGaa                                                                                                 <  1:95846/140‑1 (MQ=255)
                                    gTGCGCTGTTAGCTGCTTTTCCCCAACTTGGTGCGTTGACAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAcc                                                                                               <  2:325190/140‑1 (MQ=255)
                                    gTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGCGCAGAAAAAGGGGGGGGGTGGctgtcc                                                                                               >  1:575762/1‑135 (MQ=255)
                                    gTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACg                                                                                                                <  1:144706/123‑1 (MQ=255)
                                    gTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACg                                                                                                                >  2:144706/1‑123 (MQ=255)
                                    gTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAcc                                                                                               >  1:679874/1‑140 (MQ=255)
                                    gTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAcc                                                                                               >  2:458008/1‑140 (MQ=255)
                                    gTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAcc                                                                                               >  2:186078/1‑140 (MQ=255)
                                    gTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAcc                                                                                               >  2:8308/1‑140 (MQ=255)
                                    gTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAcc                                                                                               <  1:571652/140‑1 (MQ=255)
                                    gTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑CCCTCCCCCATTGAGTAGAAAACGGAGGTGGATGGCGAAcc                                                                                               >  1:325190/1‑140 (MQ=255)
                                    gTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑CCCACCCCCATTGAGGAGAAAACGGGGGTGGATGGCGAAcc                                                                                               >  1:226124/1‑140 (MQ=255)
                                    gTGCGCTGTTAGCTGCTGTTCCCCAACTTGGTGCGTTGACAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAcc                                                                                               <  2:271968/140‑1 (MQ=255)
                                     tGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAcc                                                                                               <  2:814985/139‑1 (MQ=255)
                                     tGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAcc                                                                                               >  1:814985/1‑139 (MQ=255)
                                          tGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGtttt                                                                                         >  2:528815/1‑140 (MQ=255)
                                            ttAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATGGAGCAGAAAACAGTCGTGGATGGCGAACCTGTTTTAc                                                                                       >  1:396773/1‑140 (MQ=255)
                                            ttAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑CCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACg                                                                                      <  1:432726/140‑1 (MQ=255)
                                             tAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACg                                                                                      <  1:128714/140‑1 (MQ=255)
                                              aGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGa                                                                                     <  1:458008/140‑1 (MQ=255)
                                                      ttCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTaa                                                                             <  1:186078/140‑1 (MQ=255)
                                                             aCTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTc                                                                      >  2:27545/1‑140 (MQ=255)
                                                                tGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCAcc                                                                   >  2:585239/1‑140 (MQ=255)
                                                                              ggTTAGATTATCCCCGACGTTTTCCGTTGATTTCGCTTGCACATCCAGAGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCATAATCTTGTCACCCGTTTCCCTTTgcg                                                     <  2:226124/140‑1 (MQ=255)
                                                                                   gATTATCCTCGACGTTACCCGTCGATATCGCATTACCATTCAGCGAAATAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCg                                                <  2:396773/140‑1 (MQ=255)
                                                                                   gATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCg                                                <  1:528815/140‑1 (MQ=255)
                                                                                        tCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTATCACCCGTTTCCCTTTGCGCAGCGGTttg                                           >  1:394362/1‑140 (MQ=255)
                                                                                           ccGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTttgttg                                        <  1:27545/140‑1 (MQ=255)
                                                                                             gACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGaa                                      <  2:679874/140‑1 (MQ=255)
                                                                                                   tCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACGCGTTTCCCTTTGTGCAGAGGTTTGTTGAATCGCGt                                >  2:330958/1‑140 (MQ=255)
                                                                                                             atCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGaa                      <  1:174516/140‑1 (MQ=255)
                                                                                                                 cTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTACGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGGGc                  >  1:283810/1‑140 (MQ=255)
                                                                                                                    gAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATg               <  1:8308/140‑1 (MQ=255)
                                                                                                                     aaCATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGc              >  1:14437/1‑140 (MQ=255)
                                                                                                                      aCATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGcc             >  2:133972/1‑140 (MQ=255)
                                                                                                                           cAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTga        >  2:84123/1‑140 (MQ=255)
                                                                                                                            aGCGAAACCGG‑‑CCCCCTCCCATTGAGCAGAAAACGGTGGTGGATGGCGCACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTgag       <  2:394362/140‑1 (MQ=255)
                                                                                                                            aGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTgag       >  2:592877/1‑140 (MQ=255)
                                                                                                                                 aaCAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGt  <  1:585239/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                        
AACAGATTTTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑CCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGT  >  NC_000913/868280‑868547

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