Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,950,559 (C)5→6 intergenic (‑37/‑273) aspS ← / → yecD aspartyl‑tRNA synthetase/isochorismatase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,950,5541.C100.0% 146.3 / NA 41intergenic (‑32/‑278)aspS/yecDaspartyl‑tRNA synthetase/isochorismatase family protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (20/21);  total (20/21)

ATGAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑CCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTAC  >  NC_000913/1950416‑1950683
                                                                                                                                           |                                                                                                                                 
aTGAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATg                                                                                                                                    >  2:132795/1‑139 (MQ=255)
                cAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAg                                                                                                                    >  1:609769/1‑139 (MQ=255)
                      cACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGAcc                                                                                                              <  1:70182/139‑1 (MQ=255)
                      cACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGAcc                                                                                                              <  1:369974/139‑1 (MQ=255)
                         gacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAg                                                                                                           >  1:2376/1‑139 (MQ=255)
                                        aaCCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                            <  2:158826/139‑1 (MQ=255)
                                        aaCCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                            <  1:132795/139‑1 (MQ=255)
                                                         tgctgcCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑cccccc                                                                                                                             >  1:243258/1‑87 (MQ=255)
                                                         tgcagcCACACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑cccccc                                                                                                                             <  2:243258/88‑2 (MQ=255)
                                                              ccccACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGaaaaa                                                                      >  1:563104/1‑139 (MQ=255)
                                                                  aCGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  <  2:2376/139‑1 (MQ=255)
                                                                         aCAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCtatt                                                           <  2:609769/139‑1 (MQ=255)
                                                                                     gTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                            >  2:1069971/1‑94 (MQ=255)
                                                                                     gTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                            <  1:1069971/94‑1 (MQ=255)
                                                                                                tCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                            >  1:1073466/1‑83 (MQ=255)
                                                                                                tCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                            <  2:1073466/83‑1 (MQ=255)
                                                                                                tCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCg                                    >  1:428368/1‑139 (MQ=255)
                                                                                                   gTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGAc                                                                                                               <  2:215986/61‑1 (MQ=255)
                                                                                                   gTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGAc                                                                                                               >  1:215986/1‑61 (MQ=255)
                                                                                                    tACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCtgctgc                                                                                                                        >  2:365589/1‑51 (MQ=255)
                                                                                                    tACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCtgctgc                                                                                                                        <  1:365589/51‑1 (MQ=255)
                                                                                                       gCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTgcgc                                                                                                                      <  2:111132/50‑1 (MQ=255)
                                                                                                       gCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTgcgc                                                                                                                      >  1:111132/1‑50 (MQ=255)
                                                                                                       gCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTa                                                                                          >  1:595111/1‑78 (MQ=255)
                                                                                                       gCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTa                                                                                          <  2:595111/78‑1 (MQ=255)
                                                                                                       gCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACGACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGc                             >  1:528403/1‑139 (MQ=255)
                                                                                                         aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGt                                                                                                                  <  2:1201714/52‑1 (MQ=255)
                                                                                                         aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGt                                                                                                                  >  1:1201714/1‑52 (MQ=255)
                                                                                                         aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCt                           >  2:430536/1‑139 (MQ=255)
                                                                                                                   ccTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCACAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCt                 >  2:1270427/1‑139 (MQ=255)
                                                                                                                      tAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGAcc                                                                                                              <  1:721188/43‑1 (MQ=38)
                                                                                                                      tAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGAcc                                                                                                              >  2:721188/1‑43 (MQ=38)
                                                                                                                      tAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGc                                                                                                    >  2:999264/1‑53 (MQ=38)
                                                                                                                      tAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGc                                                                                                    <  1:999264/53‑1 (MQ=38)
                                                                                                                       aaCTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAg                                                                                               <  1:889501/57‑1 (MQ=39)
                                                                                                                       aaCTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAg                                                                                               >  2:889501/1‑57 (MQ=39)
                                                                                                                            aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  <  2:376107/81‑1 (MQ=255)
                                                                                                                            aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  >  1:376107/1‑81 (MQ=255)
                                                                                                                            aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCa                                       <  2:216589/108‑1 (MQ=255)
                                                                                                                            aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCa                                       >  1:216589/1‑108 (MQ=255)
                                                                                                                             gCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTc       >  2:1213371/1‑139 (MQ=255)
                                                                                                                                gCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGtt    <  1:1270427/139‑1 (MQ=255)
                                                                                                                                 ccGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa   <  2:428368/139‑1 (MQ=255)
                                                                                                                                  cggcggATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTAc  <  1:1213371/139‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                 
ATGAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑CCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTAC  >  NC_000913/1950416‑1950683

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