Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 868,748 Δ1 bp coding (1229/1872 nt) gsiA → glutathione ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913868,7460C.95.8% 202.7 / 5.2 49coding (1227/1872 nt)gsiAglutathione ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (20/26);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

AGGGCGGCGAAATTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTG  >  NC_000913/868649‑868840
                                                                                                 |                                                                                              
aGGGCGGCGAAATTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGAcct                                                                                              <  1:727316/100‑2 (MQ=255)
 gggCGGCGAAATTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGAcctt                                                                                             >  1:1113189/1‑98 (MQ=255)
   gcggcgAAATTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGAccttac                                                                                           >  2:1141224/1‑96 (MQ=255)
        gAAATTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTc                                                                                     <  1:118631/100‑1 (MQ=255)
         aaaTTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCg                                                                                    <  2:12367/100‑1 (MQ=255)
            ttATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTg                                                                                 <  1:832752/100‑1 (MQ=255)
            ttATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTg                                                                                 >  1:1239351/1‑100 (MQ=255)
               tCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGAc                                                                              <  1:825866/100‑1 (MQ=255)
               tCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGAc                                                                              <  2:273035/100‑1 (MQ=255)
                   tAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCAc                                                                          >  1:1007629/1‑100 (MQ=255)
                    aaCGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACg                                                                         <  1:1107013/100‑1 (MQ=255)
                    aaCGGGTAGCGAATCGATACGTTGTCACCCGGCAAACTTCAGGCATTACGCGGGGATATTCAGTTTATTTTTCAGGTCCTTACGCTTCGCTGGACCCACg                                                                         <  2:461304/100‑1 (MQ=255)
                      cGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCCGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCAAGTc                                                                       <  2:224224/100‑1 (MQ=255)
                       ggTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCa                                                                      <  2:623000/100‑1 (MQ=255)
                          cAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGAc                                                                   <  1:893962/100‑1 (MQ=255)
                           aGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGAcc                                                                  >  2:61008/1‑100 (MQ=255)
                           aGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGATTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGAc                                                                   >  2:833437/1‑100 (MQ=255)
                              gAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATc                                                               <  2:1209764/100‑1 (MQ=255)
                              gAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATc                                                               >  2:1285458/1‑100 (MQ=255)
                                 tCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGt                                                            <  2:774409/100‑1 (MQ=255)
                                  cGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTg                                                           >  1:403558/1‑100 (MQ=255)
                                       ccTTGTCACCCGGCAAACTTCAGACATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCg                                                      >  1:49148/1‑100 (MQ=255)
                                           gTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTa                                                  <  1:580997/100‑1 (MQ=255)
                                             cACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATc                                                >  2:1031256/1‑100 (MQ=255)
                                                ccGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGaa                                             >  1:210300/1‑100 (MQ=255)
                                                      aaCTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTg                                       <  1:1277157/100‑1 (MQ=255)
                                                        cTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGGTTCGATAATCGaatcaccgat                                     <  1:416354/100‑9 (MQ=255)
                                                        cTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCg                                     <  1:1221163/100‑1 (MQ=255)
                                                           cAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCgtgt                                  >  2:601530/1‑100 (MQ=255)
                                                            aGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATAGAACCGCTGCGTGTa                                 >  1:1286387/1‑100 (MQ=255)
                                                             ggCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTac                                >  1:964111/1‑100 (MQ=255)
                                                             ggCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTac                                >  1:870852/1‑100 (MQ=255)
                                                                   cttccgatctatTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGtt                          <  2:203998/91‑1 (MQ=255)
                                                                   aCGCCGGGATATTCACTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGtt                          <  1:1285469/100‑1 (MQ=255)
                                                                     gCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGttt                         <  1:290553/100‑1 (MQ=255)
                                                                     gCCGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTa                        >  2:989663/1‑100 (MQ=255)
                                                                       cGGGATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTAtt                      <  2:1015863/100‑1 (MQ=255)
                                                                        gggATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTg                     <  1:810072/100‑1 (MQ=255)
                                                                        gggATATTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTg                     >  2:529247/1‑100 (MQ=255)
                                                                            tatTCAGTTTATTTTTCAGGACCTTTCGCTTCGCTGGAGACACGTCATACTATCGGTGATTCTATTATCGAACCGCTGCGTGTATACGGTTTAtt                      <  2:450416/95‑1 (MQ=255)
                                                                            tatTCAGTTTATTTTGCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCa                  <  2:1056618/99‑1 (MQ=255)
                                                                             atTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGttagatcggaag                >  1:204000/1‑90 (MQ=255)
                                                                             atTCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAgg                <  1:587784/100‑1 (MQ=255)
                                                                              ttCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGt               >  2:981582/1‑100 (MQ=255)
                                                                               tCAGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTa              <  2:210298/100‑1 (MQ=255)
                                                                                 aGTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCCGGTaaa            >  2:594232/1‑100 (MQ=255)
                                                                                  gTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAg           >  1:120431/1‑100 (MQ=255)
                                                                                  gTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAg           >  1:1142153/1‑100 (MQ=255)
                                                                                  gTTTATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAg           <  2:928419/100‑1 (MQ=255)
                                                                                   tttATTTTTCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGa          <  2:1093237/100‑1 (MQ=255)
                                                                                       tttttCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGAcgcg      <  1:130254/100‑1 (MQ=255)
                                                                                           tCAGGACCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTg  >  2:427171/1‑100 (MQ=255)
                                                                                                 |                                                                                              
AGGGCGGCGAAATTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTG  >  NC_000913/868649‑868840

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: