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)
        gAAATTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTc                                                                                     <  1:118631/100‑1 (MQ=255)
         aaaTTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCg                                                                                    <  2:12367/100‑1 (MQ=255)
            ttATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTg                                                                                 <  1:832752/100‑1 (MQ=255)
            ttATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTg                                                                                 >  1:1239351/1‑100 (MQ=255)
               tCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGAc                                                                              <  1:825866/100‑1 (MQ=255)
               tCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGAc                                                                              <  2:273035/100‑1 (MQ=255)
                   tAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCAc                                                                          >  1:1007629/1‑100 (MQ=255)
                    aaCGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACg                                                                         <  1:1107013/100‑1 (MQ=255)
                    aaCGGGTAGCGAATCGATACGTTGTCACCCGGCAAACTTCAGGCATTACGCGGGGATATTCAGTTTATTTTTCAGG‑TCCTTACGCTTCGCTGGACCCACg                                                                         <  2:461304/100‑1 (MQ=255)
                      cGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCCGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCAAGTc                                                                       <  2:224224/100‑1 (MQ=255)
                       ggTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCa                                                                      <  2:623000/100‑1 (MQ=255)
                          cAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGAc                                                                   <  1:893962/100‑1 (MQ=255)
                           aGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGAcc                                                                  >  2:61008/1‑100 (MQ=255)
                           aGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGATTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGAc                                                                   >  2:833437/1‑100 (MQ=255)
                              gAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATc                                                               <  2:1209764/100‑1 (MQ=255)
                              gAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATc                                                               >  2:1285458/1‑100 (MQ=255)
                                 tCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGt                                                            <  2:774409/100‑1 (MQ=255)
                                  cGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTg                                                           >  1:403558/1‑100 (MQ=255)
                                       ccTTGTCACCCGGCAAACTTCAGACATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCg                                                      >  1:49148/1‑100 (MQ=255)
                                           gTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTa                                                  <  1:580997/100‑1 (MQ=255)
                                             cACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATc                                                >  2:1031256/1‑100 (MQ=255)
                                                ccGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGaa                                             >  1:210300/1‑100 (MQ=255)
                                                      aaCTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTg                                       <  1:1277157/100‑1 (MQ=255)
                                                        cTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGGTTCGATAATCGaatcaccgat                                     <  1:416354/100‑9 (MQ=255)
                                                        cTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCg                                     <  1:1221163/100‑1 (MQ=255)
                                                           cAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCgtgt                                  >  2:601530/1‑100 (MQ=255)
                                                            aGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATAGAACCGCTGCGTGTa                                 >  1:1286387/1‑100 (MQ=255)
                                                             ggCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTac                                >  1:964111/1‑100 (MQ=255)
                                                             ggCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTac                                >  1:870852/1‑100 (MQ=255)
                                                                   cttccgatctatTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGtt                          <  2:203998/91‑1 (MQ=255)
                                                                   aCGCCGGGATATTCACTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGtt                          <  1:1285469/100‑1 (MQ=255)
                                                                     gCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGttt                         <  1:290553/100‑1 (MQ=255)
                                                                     gCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTa                        >  2:989663/1‑100 (MQ=255)
                                                                       cGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTAtt                      <  2:1015863/100‑1 (MQ=255)
                                                                        gggATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTg                     <  1:810072/100‑1 (MQ=255)
                                                                        gggATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTg                     >  2:529247/1‑100 (MQ=255)
                                                                            tatTCAGTTTATTTTTCAGGA‑CCTTTCGCTTCGCTGGAGACACGTCATACTATCGGTGATTCTATTATCGAACCGCTGCGTGTATACGGTTTAtt                      <  2:450416/95‑1 (MQ=255)
                                                                            tatTCAGTTTATTTTGCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCa                  <  2:1056618/99‑1 (MQ=255)
                                                                             atTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGttagatcggaag                >  1:204000/1‑90 (MQ=255)
                                                                             atTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAgg                <  1:587784/100‑1 (MQ=255)
                                                                              ttCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGt               >  2:981582/1‑100 (MQ=255)
                                                                               tCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTa              <  2:210298/100‑1 (MQ=255)
                                                                                 aGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCCGGTaaa            >  2:594232/1‑100 (MQ=255)
                                                                                  gTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAg           >  1:120431/1‑100 (MQ=255)
                                                                                  gTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAg           >  1:1142153/1‑100 (MQ=255)
                                                                                  gTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAg           <  2:928419/100‑1 (MQ=255)
                                                                                   tttATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGa          <  2:1093237/100‑1 (MQ=255)
                                                                                       tttttCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGAcgcg      <  1:130254/100‑1 (MQ=255)
                                                                                           tCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTg  >  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: ‑