Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 868,871 G→C G451A (GGC→GCC)  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,8710GC96.4% 86.6 / ‑5.7 28G451A (GGC→GCC) gsiAglutathione ABC transporter ATP binding subunit GsiA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (15/12);  minor base A (0/1);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTC  >  NC_000913/868794‑868945
                                                                             |                                                                          
cGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGCCCt                                                                         >  1:182208/1‑81 (MQ=255)
           gtgtACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGaa                                                               >  1:182213/1‑80 (MQ=255)
              tACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATg                                                           <  2:182165/81‑1 (MQ=255)
                cacGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGcc                                                         <  1:182166/81‑1 (MQ=255)
                 acGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCt                                                        >  1:182216/1‑81 (MQ=255)
                   ggTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCTgg                                                      >  2:182197/1‑81 (MQ=255)
                    gTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCTGGc                                                     >  1:182209/1‑81 (MQ=255)
                         ttGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCTGGCGTTAc                                                >  2:182206/1‑81 (MQ=255)
                           gCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCTGGCGTTAcc                                               <  2:182167/80‑1 (MQ=255)
                           gCCAGGTAAAGACGCGGCTGAACGCGTTGCGTGGTTGCTGGAGCGCGTGGACCTGTTACCTGAACATGCCTGGCGTTAcc                                               <  2:182168/80‑1 (MQ=255)
                             cAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCTGGCGTTACCCGc                                            <  1:182169/81‑1 (MQ=255)
                                gTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCTGGCGTTACCCGCATg                                         <  2:182170/81‑1 (MQ=255)
                                  aaaGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAg                                       >  2:182211/1‑81 (MQ=255)
                                   aaGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGt                                      >  2:182210/1‑81 (MQ=255)
                                        gcgGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTcc                                 >  2:182217/1‑81 (MQ=255)
                                         cgGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCcg                                <  2:182171/81‑1 (MQ=255)
                                              gCACGCGTTGCGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGt                           >  2:182215/1‑81 (MQ=255)
                                                aCGCGTTGCGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGt                           >  1:182225/1‑79 (MQ=255)
                                                    gTTGCGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTcagcgc                     >  2:182207/1‑81 (MQ=255)
                                                     ttGCGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTcagcgcc                    <  1:182172/81‑1 (MQ=255)
                                                        cGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTcagcgccag                  <  2:182174/80‑1 (MQ=255)
                                                        cGTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTcagcgccag                  <  1:182173/80‑1 (MQ=255)
                                                         gTGGTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTcagcgccagcg                <  2:182175/81‑1 (MQ=255)
                                                            gTTGCTGGAGCGCGTGGCCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCa              <  1:182176/80‑1 (MQ=255)
                                                                  ggAGCGCGTGGCCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCAt       >  2:182205/1‑81 (MQ=255)
                                                                    aGCGCGTGGCCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTg     >  2:182214/1‑81 (MQ=255)
                                                                      cgcgTGGCCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCt   >  2:182222/1‑81 (MQ=255)
                                                                       gcgTGGCCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTc  <  1:182177/81‑1 (MQ=255)
                                                                             |                                                                          
CGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTC  >  NC_000913/868794‑868945

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