Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,540,372 (GC)5→4 coding (442‑443/1098 nt) cysA ← sulfate/thiosulfate ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,540,3640G.100.0% 107.6 / NA 26coding (451/1098 nt)cysAsulfate/thiosulfate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (14/12);  total (14/12)
*NC_0009132,540,3650C.100.0% 107.6 / NA 26coding (450/1098 nt)cysAsulfate/thiosulfate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (14/12);  total (14/12)

GATCGTGGGTCACAAAAACGCTGGTGAATTTTAGTTCTTCATGGAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCAGCGCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATGGGCAAGCTGGACCATTTCCAGCAATTTTGTCACTTTCGCTTTGATGGCTG  >  NC_000913/2540226‑2540483
                                                                                                                                          ||                                                                                                                      
gATCGTGGGTCACAAAAACGCTGGTGAATTTTAGTTCTTCATGGAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCAcagccagc                                                                                                                        >  1:353090/1‑140 (MQ=255)
 aTCGTGGGTCACAAAAACGCTGGTGAATTTTAGTTCTTCATGGAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCAgcg                                                                                                                       <  2:708110/140‑1 (MQ=255)
  tCGTGGGTCACAAAAACGCTGGTGAATTTTAGTTCTTCATGGAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCAgcgc                                                                                                                      >  2:391507/1‑140 (MQ=255)
     tGGGTCACAAAAACGCTGGTGAATTTTAGTTCTTCATGGAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCAgcgcgcg                                                                                                                   >  1:1049847/1‑140 (MQ=255)
     tGGGTCACAAAAACGCTGGTGAATTTTAGTTCTTCATGGAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCAgcgcgcg                                                                                                                   >  2:1061867/1‑140 (MQ=255)
      gggTCACAAAAACGCTGGTGAATTTTAGTTCTTCATGGAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCAgcgcgcg                                                                                                                   <  1:198029/139‑1 (MQ=255)
      gggTCACAAAAACGCTGGTGAATTTTAGTTCTTCATGGAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCAgcgcgcg                                                                                                                   >  2:198029/1‑139 (MQ=255)
       ggTCACAAAAACGCTGGTGAATTTTAGTTCTTCATGGAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCAgcgcgcgcc                                                                                                                 >  1:845737/1‑139 (MQ=255)
        gTCACAAAAACGCTGGTGAATTTTAGTTCTTCATGGAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCAgcgcgcgcca                                                                                                                >  1:240277/1‑138 (MQ=255)
         tCACAAAAACGCTGGTGAATTTTAGTTCTTCATGGAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCAgcgcgcgccag                                                                                                               >  1:1018491/1‑137 (MQ=255)
          cacaAAAACGCTGGTGAATTTTAGTTCTTCATGGAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCAgcgcgcgccagc                                                                                                              <  2:353090/140‑5 (MQ=255)
           tcaAAAACGCTGGTGAATTTTAGTTCTTCATGGAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCAgcgcgcgccagcg                                                                                                             >  2:787107/2‑135 (MQ=255)
                  cGCTGGTGAATTTTAGTTCTTCATGGAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCAcgcg                                                                                                    <  1:1061867/140‑1 (MQ=255)
                  cGCTGGTGAATTTTAGTTCTTCATGGAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCAcgcg                                                                                                    >  1:48390/1‑140 (MQ=255)
                  cGCTGGTGAATTTTAGTTCTTCATGGAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCAcgcg                                                                                                    <  2:48390/140‑1 (MQ=255)
                      ggTGAATTTTAGTTCTTCATGGAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGt                                                                                                <  2:893774/140‑1 (MQ=255)
                                           gAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTgc                                                                           <  1:741633/140‑1 (MQ=255)
                                            aGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTgcg                                                                          <  1:787107/140‑1 (MQ=255)
                                                  cGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGAt                                                                    >  2:415846/1‑140 (MQ=255)
                                                   gCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATa                                                                   >  2:3177/1‑140 (MQ=255)
                                                                cAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGa                                                      >  2:717058/1‑140 (MQ=255)
                                                                             aCCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATGGGCAAGCTGGa                                         >  2:269278/1‑140 (MQ=255)
                                                                               cTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATggg                                                  <  1:651114/129‑1 (MQ=255)
                                                                               cTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATggg                                                  >  2:651114/1‑129 (MQ=255)
                                                                               cTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATGGGCAAGCTGGAcc                                       >  2:801234/1‑140 (MQ=255)
                                                                                tGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATGGGCAAGCTGGACCa                                      >  2:414951/1‑140 (MQ=255)
                                                                                tGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATGGGCAAGCTGGACCa                                      >  2:647180/1‑140 (MQ=255)
                                                                                  cgcATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATggg                                                  <  2:347060/126‑1 (MQ=255)
                                                                                  cgcATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATggg                                                  >  1:347060/1‑126 (MQ=255)
                                                                                      tCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATGGGCAAGCTGGACCATTTCCa                                <  2:240277/140‑1 (MQ=255)
                                                                                         aGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATGGGCAAGCTGGACCATTTCCAGCa                             >  1:816910/1‑140 (MQ=255)
                                                                                           cgcgcCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCATGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATGGGCAAGCTGGACCATTTCCAGCAAt                           >  2:608475/1‑140 (MQ=255)
                                                                                           cgcgcCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATGGGCAAGCTGGACCATTTCCAGCAAt                           >  2:305188/1‑140 (MQ=255)
                                                                                               ccAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATGGGCAAGCTGGACCATTTCCAGCAATTTTg                       <  1:1121855/140‑1 (MQ=255)
                                                                                                cAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATGGGCAAGCTGGACCATTTCCAGCAATTTTGt                      <  1:717058/140‑1 (MQ=255)
                                                                                                               gcagcaGAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATGGGCAAGCTGGACCATTTCCAGCAATTTTGTCACTTTCGCTTTGAt       <  2:816910/140‑1 (MQ=255)
                                                                                                                 agcagAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATGGGCAAGCTGGACCATTTCCAGCAATTTTGTCACTTTCGCTTTGATgg     >  1:518127/1‑140 (MQ=255)
                                                                                                                    agAATTTGCGGTTCCACAGCCA‑‑GCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATGGGCAAGCTGGACCATTTCCAGCAATTTTGTCACTTTCGCTTTGATGGCTg  <  1:125319/140‑1 (MQ=255)
                                                                                                                                          ||                                                                                                                      
GATCGTGGGTCACAAAAACGCTGGTGAATTTTAGTTCTTCATGGAGTTGACGCAGCCAGCGACGCAGCTCTTTACGCACCTGCGCATCCAGCGCGCCAAACGGTTCATCAAGCAGCAGAATTTGCGGTTCCACAGCCAGCGCGCGCGCCAGCGCCACGCGCTGTTTCTGGCCGCCGGAAAGCTGCGCCGGATAACGATCCGCCAGATGGGCAAGCTGGACCATTTCCAGCAATTTTGTCACTTTCGCTTTGATGGCTG  >  NC_000913/2540226‑2540483

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