Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 868,422 (C)8→7 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,4220C.96.0% 99.4 / ‑3.7 25coding (903/1872 nt)gsiAglutathione ABC transporter ATP binding subunit GsiA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (10/14);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GGCGTGGTGGCAGAGATTGCCGATCGGGTACTGGTGATGTATCAGGGCGAGGCGGTGGAAACGGGTACCGTCGAACAGATTTTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGTCAGTTTTGATCTCTGGCCTGGCGAAACGCTATCGCTGGTGGGCGAGTCTGGCAGCGGTAAATCCACTACCGGGCGGGCGTTGCTGCGCCTGGTCGAATCGCAGGGCGGCGAAATTATCT  >  NC_000913/868207‑868666
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
ggCGGGGTGGCAGAAATTGCCGATCGGGTATGGTTGATTAATAAGGGCGAGGCGTTGAAAACGGGTACCGTCGAACAGATTTTTCATGCACCGAAACATCCTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGTGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGttt                                                                                                                                                                                                                 <  1:252247/251‑1 (MQ=255)
 gCGTGGTGGCAGAGATTGCCGATCGGGTACTGGTGATGTATCAGGGCGAGGCGGTGGAAACGGGTACCGTCGAACAGATTTTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGTGTTAGATTATCCCCGACGTTTCCCGTTGATAGCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGttt                                                                                                                                                                                                                 <  1:151043/251‑1 (MQ=255)
         gCAGAGATTGCCGATCGGGTACTGGTGATGTATCAGGGCGAGGCGGTGGAAACGGGTACCGTCGAACAGATTTTTCATGCACCGCAACATCCTTACCCCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTg                                                                                                                                                                                                         <  1:165968/251‑1 (MQ=255)
                           gTACTGGTGATGTATCAGGGCGAGGCGGTGGAAACGGGTAACGTCGAACAGATTTTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGGAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGGGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGt                                                                                                                                                                                       >  1:95886/1‑251 (MQ=255)
                                ggTGATGTATCAGGGCGAGGCGGTGGAAACGGGTACCGTCGAACAGATTTTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCCAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAATGCGGAACCTTGTCCCCCGGTTccc                                                                                                                                                                                  >  2:258183/1‑251 (MQ=255)
                                                   ggCGGTGGAAACGGGTACCGTCGAACAGATTTTTCATGCACCGCAAACTCCTGACGCCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCTCTTCCCTTTGCGCAGCGGTgtgtt                                                                                                                                                                <  2:65/251‑1 (MQ=255)
                                                           aaaCGGGTACCGTCGAACAGATTTTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATcgcg                                                                                                                                                        >  2:275034/1‑251 (MQ=255)
                                                                          aaCAGATTTTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGc                                                                                                                                         >  1:175859/1‑251 (MQ=255)
                                                                                    tCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTg                                                                                                                                <  1:54914/250‑1 (MQ=255)
                                                                                                                            gctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGTCAGTTTTGATCTCTGGCCTGGCGAAACGCTATc                                                                                        <  1:171136/250‑1 (MQ=255)
                                                                                                                                          tGGTGCGATTAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCTCTTGAACATCCATCGAAACAGGTCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGGTTTACGAGTGCGTAATCTTGTCACTCGGTTCCCTTTGCTCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTTATAAAGTCAGTTGTGATCTCTGGCCTGGCGAAACGCTATCGCTGGTGGGCGAg                                                                           <  1:314449/249‑1 (MQ=255)
                                                                                                                                               cGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGTCAGTTTTGATCTCTGGCCTGGCGAAACGCTATCGCTGGTGGGCGAGTCTGGc                                                                     >  2:96994/1‑250 (MQ=255)
                                                                                                                                                           tAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCAGCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAAcg                                                                                                                                                   <  1:297760/160‑1 (MQ=255)
                                                                                                                                                           tAGATTATCCCCGACGTTTACCGTTGATATCGCTTGAACATCCAGCGAACCAGGCCCCACCCTCGAGCAGAAAACGATGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTAACCAGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTCAcg                                                                                                                                                   >  2:297760/1‑160 (MQ=255)
                                                                                                                                                            aGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGTCAGTTTTGATCTCTGGCCTGGCGAAACGCTATCGCTGGTGGGCGAGTCTGGCAGCGGTAAATCCAc                                                       <  1:258183/251‑1 (MQ=255)
                                                                                                                                                                    cccGACGTTTCCCGTTGATAGCGCTAGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAACCGGTGGTGGATGGCGAATCTGTTTTACGAGTGCGTAATCCTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGTCAGTTTTGATCTCTGGCCTGGCGAAACGCTATCGCTGTTGGGCGAGTCTTGCAGCGGTAAATccccttcctggc                                               <  2:31508/251‑11 (MQ=255)
                                                                                                                                                                          gTTTCCAGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCAGTTGGGAACGTCAGTTTTGATCTCTGGCCTGGCGAAACGCTATCGCTGGTGGGCGGGTCTGGCAGCGGTAAATCCACTCCCGGGCGGGCGt                                         >  2:199224/1‑251 (MQ=255)
                                                                                                                                                                              cccGTTGATATCGCTTGAACATCCAGCGAAACACGCCCCCCCATCGAGCAGATACCGGTGGGTTAGGGCGGACCGGGGTTGCAGGTGGGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAGTCGCGGAACGAGGGAAGTGCCTGCCGTAGAGAAAGTCAGTTTTAAGCTCTGGCCTGTCGAATCGCTATCGCTGGTGGGCAAGACTGGCAGGGGTAAACCCTCTCCCGTGCGGGCGTtgct                                     <  2:229706/251‑1 (MQ=255)
                                                                                                                                                                                 gTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGTCAGTTTTGATCTCTGGACTGGCGAAACGCTATCGCTGGTGGGCGAGTCTGGCAGCGGTAAATCCACTACcgggcgg                                             >  1:148081/1‑240 (MQ=255)
                                                                                                                                                                                 gTTGATATCCCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGTCAGTTTTGATCTCTGGCCTGGCTAAACGCGAGCGCTGGTGGGCGAGTCTGGCAGCTGTAAATCCACTACcgggcgg                                             <  2:148081/240‑1 (MQ=255)
                                                                                                                                                                                      tatCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTTGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGGAAGTCAGTTTTGATCTCTTGCCTGGCGAAACGCTATCGCTGGTGGGCGAGTCTGGCAGCTGTAAATCCACTACCGGGCGGGCGTTGCTGCGCCTgg                             >  2:147123/1‑251 (MQ=255)
                                                                                                                                                                                              gAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGTCAGTTTTGATCTCTGGCCTGGCGAAACGCTATCGCTGGTGGGCGAGTCTGGCAGCGGTAAATCCACTACCGGGCGGGCGTTGCTGCGCCTGGTCGAATc                      <  1:96994/250‑1 (MQ=255)
                                                                                                                                                                                                    ccatcGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTGGTCCCCCGTTTCCCTTTGCGCATCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGGTGAGAAAGTCAGTTTTGATCTCTTGCCTGGCGAAACGCTATCGCTGGTGGGCGAGTCTGGCAGCGGTAAATCCACTACCGGGCGGGCGTTGCTGCGCCTGGTCGAATCGCAgcgc               <  2:95886/247‑4 (MQ=255)
                                                                                                                                                                                                     cAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGCTGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGTCAGTTTTGATCTCTGGCCTGGCGAAACGCTATCGCTGGTGGGCGAGTCTGGCAGCGGTAAATCCACTACCGGGCGGGCGTTGCTGCGCCTGGTCGAATCGCAGggcg              >  1:293799/1‑251 (MQ=255)
                                                                                                                                                                                                     cAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGTCAGTTTTGATCTCTGGCCTGGCGAAACGCTATCGCTGGTGGGCGAGTCTGGCAGCGGTAAATCCACTACCGGGCGGGCGTTGCTGCGCCTGGTCGAATCGCAGggcg              >  1:14748/1‑251 (MQ=255)
                                                                                                                                                                                                                  cccccccATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGTCAGTTTTGATCTCTGGCCTGGCGAAACGCTATCGCTGGTGGGCGAGTCTGGCAGCGGTAAATCCACTACCGGGCGGGCGTTGCTGCGCCTGGTCGAATCGCAGGGCGGCGAAATTATCt  <  1:275034/251‑1 (MQ=255)
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
GGCGTGGTGGCAGAGATTGCCGATCGGGTACTGGTGATGTATCAGGGCGAGGCGGTGGAAACGGGTACCGTCGAACAGATTTTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGTCAGTTTTGATCTCTGGCCTGGCGAAACGCTATCGCTGGTGGGCGAGTCTGGCAGCGGTAAATCCACTACCGGGCGGGCGTTGCTGCGCCTGGTCGAATCGCAGGGCGGCGAAATTATCT  >  NC_000913/868207‑868666

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