Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,542,393 G→A A40A (GCC→GCT cysU ← sulfate/thiosulfate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,542,3930GA100.0% 96.0 / NA 28A40A (GCC→GCTcysUsulfate/thiosulfate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/12);  total (16/12)

GCCATCAGCAGACCGAAAACGCCGTTAAAAATCGATGCCACAAACGCCGACAGCAGCGTTACTTTGTAGGCCGCGACCACCTGCGGGTTGGTGATCACCTCCCAGTACTGCGCCCAGCTCATCTGGGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCAAAATCAGGCACACAAACAGCAGACTGGTGCCGAGGCTTAAGGTAAAGCCCGGCAGCACGCGTCTGGAGGAGACAGCAAACATCAGTT  >  NC_000913/2542268‑2542517
                                                                                                                             |                                                                                                                            
gCCATCAGCAGACCGAAAACGCCGTTAAAAATCGATGCCACAAACGCCGACAGCAGCGTTACTTTGTAGGCCGCGACCACCTGCGGGTTGGTGATCACCTCCCAGTACTGCGCCCAGCTCATCTGAGCCAGt                                                                                                                        >  2:343079/1‑132 (MQ=255)
gCCATCAGCAGACCGAAAACGCCGTTAAAAATCGATGCCACAAACGCCGACAGCAGCGTTACTTTGTAGGCCGCGACCACCTGCGGGTTGGTGATCACCTCCCAGTACTGCGCCCAGCTCATCTGAGCCAGt                                                                                                                        <  1:343079/132‑1 (MQ=255)
    tCAGCAGACCGAAAACGCCGTTAAAAATCGATGCCACAAACGCCGACAGCAGCGTTACTTTGTAGGCCGCGACCACCTGCGGGTTGGTGATCACCTCCCAGTACTGCGCCCAGCTCATCTGAGCCAGTTGCATCACCAg                                                                                                             >  1:56473/1‑139 (MQ=255)
            ccGAAAACGCCGTTAAAAATCGATGCCACAAACGCCGACAGCAGCGTTACTTTGTAGGCCGCGACCACCTGCGGGTTGGTGATCACCTCCCAGTACTGCGCCCAGCTCATCTGAGCCAGTTGCATCa                                                                                                                 <  2:14703/127‑1 (MQ=255)
            ccGAAAACGCCGTTAAAAATCGATGCCACAAACGCCGACAGCAGCGTTACTTTGTAGGCCGCGACCACCTGCGGGTTGGTGATCACCTCCCAGTACTGCGCCCAGCTCATCTGAGCCAGTTGCATCa                                                                                                                 >  1:14703/1‑127 (MQ=255)
                              aTCGATGCCACAAACGCCGACAGCAGCGTTACTTTGTAGGCCGCGACCACCTGCGGGTTGGTGATCACCTCCCAGTACTGCGCCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCAAAATCa                                                                                   >  2:185568/1‑139 (MQ=255)
                                        caAACGCCGACAGCAGCGTTACTTTGTAGGCCGCGACCACCTGCGGGTTGGTGATCACCTCCCAGTACTGCGCCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCAAAATCAGGcac                                                                              <  1:74464/134‑1 (MQ=255)
                                        caAACGCCGACAGCAGCGTTACTTTGTAGGCCGCGACCACCTGCGGGTTGGTGATCACCTCCCAGTACTGCGCCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCAAAATCAGGcac                                                                              >  2:74464/1‑134 (MQ=255)
                                             gCCGACAGCAGCGTTACTTTGTAGGCCGCGACCACCTGCGGGTTGGTGATCACCTCCCAGTACTGCGCCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCAAAATCAGGCACACAAAcagca                                                                    <  1:185568/139‑1 (MQ=255)
                                                                  tAGGCCGCGACCACCTGCGGGTTGGTGATCACCTCCCAGTACTGCGCCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCAAAATCAGGCACACAAACAGCAGACTGGTGCCGAGGCTTAAgg                                               >  1:682528/1‑139 (MQ=255)
                                                                       cgcgACCACCTGCGGGTTGGTGATCACCTCCCAGTACTGCGCCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGagag                                                                                                    <  2:510042/81‑1 (MQ=255)
                                                                       cgcgACCACCTGCGGGTTGGTGATCACCTCCCAGTACTGCGCCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGagag                                                                                                    >  1:510042/1‑81 (MQ=255)
                                                                       cgcgACCACCTGCGGGTTGGTGATCACCTCCCAGTACTGCGCCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCgg                                                                                                 <  2:291393/84‑1 (MQ=255)
                                                                       cgcgACCACCTGCGGGTTGGTGATCACCTCCCAGTACTGCGCCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCgg                                                                                                 >  1:291393/1‑84 (MQ=255)
                                                                                           tGATCACCTCCCAGTACTGCGCCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCAAAATCAGGCACACAAACAGCAGACTGGTGCCGAGGCTTAAGGTAAAGCCCGGCAGCACGCGTCTgga                      <  1:640931/139‑1 (MQ=255)
                                                                                                            tGCGCCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGgc                                                                                                <  1:357440/48‑1 (MQ=255)
                                                                                                            tGCGCCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGgc                                                                                                >  2:357440/1‑48 (MQ=255)
                                                                                                            tGCGCCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCAAAATCAGGCACACAAACAGCAGACTGGTGCCg                                                         >  1:681897/1‑87 (MQ=255)
                                                                                                            tGCGCCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCAAAATCAGGCACACAAACAGCAGACTGGTGCCg                                                         <  2:681897/87‑1 (MQ=255)
                                                                                                            tGCGCCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCAAAATCAGGCACACAAACAGCAGACTGGTGCCGAGGCTTAAGGTAAAGCCCGGCAGCACGCGTCTGGAGGAGACAGCAAACATCa     >  1:325124/1‑139 (MQ=255)
                                                                                                             gcgcCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCAAAATCAGGCACACAAACAGCAGACTGGTGCCGAGGCTTAAGGTAAAGCCCGGCAGCACGCGTCTGgaggag                  <  1:628903/125‑1 (MQ=255)
                                                                                                             gcgcCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCAAAATCAGGCACACAAACAGCAGACTGGTGCCGAGGCTTAAGGTAAAGCCCGGCAGCACGCGTCTGgaggag                  >  2:628903/1‑125 (MQ=255)
                                                                                                              cgcCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCTAAATCAGGCACACAAACAGCAGACTGGTGCCGAGGCTTAAGGTAAAGCCCGGCAGCTCGCGTCTGGAGGAGACAGCAAACATCAGt   >  2:472612/1‑139 (MQ=255)
                                                                                                              cgcCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCAAAATCAGGCACACAAACAGCAGACTGGTGCCGAGGCTTAAGGTAAAGCCCGGCAGCACGCGTCTGGAGGAGACAGCAAACATc      >  2:581419/1‑136 (MQ=255)
                                                                                                              cgcCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCAAAATCAGGCACACAAACAGCAGACTGGTGCCGAGGCTTAAGGTAAAGCCCGGCAGCACGCGTCTGGAGGAGACAGCAAACATc      <  1:581419/136‑1 (MQ=255)
                                                                                                               gcCCAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCAAAATCAGGCACACAAACAGCAGACTGGTGCCGAGGCTTAAGGTAAAGCCCGGCAGCACGCGTCTGGAGGAGACAGCAAACATCAGtt  >  2:162012/1‑139 (MQ=255)
                                                                                                                  cAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCAAAATCAGGCACACAAACAGCAGACTGGTGCCGAGGCTTAAGGTAAAGCCCGGCAGCAcg                              >  1:56465/1‑108 (MQ=255)
                                                                                                                  cAGCTCATCTGAGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCAAAATCAGGCACACAAACAGCAGACTGGTGCCGAGGCTTAAGGTAAAGCCCGGCAGCAcg                              <  2:56465/108‑1 (MQ=255)
                                                                                                                             |                                                                                                                            
GCCATCAGCAGACCGAAAACGCCGTTAAAAATCGATGCCACAAACGCCGACAGCAGCGTTACTTTGTAGGCCGCGACCACCTGCGGGTTGGTGATCACCTCCCAGTACTGCGCCCAGCTCATCTGGGCCAGTTGCATCACCAGCGCGGAGAGCGGCAGCAGCAAAATCAGGCACACAAACAGCAGACTGGTGCCGAGGCTTAAGGTAAAGCCCGGCAGCACGCGTCTGGAGGAGACAGCAAACATCAGTT  >  NC_000913/2542268‑2542517

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