Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 28,269 C→T intergenic (+62/‑105) rihC → / → dapB ribonucleoside hydrolase 3/dihydrodipicolinate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091328,2690CT100.0% 150.9 / NA 48intergenic (+62/‑105)rihC/dapBribonucleoside hydrolase 3/dihydrodipicolinate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/26);  total (22/26)

ATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCC  >  NC_000913/28144‑28392
                                                                                                                             |                                                                                                                           
aTTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTa                                                                                                                <  1:1195889/139‑1 (MQ=255)
    gATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAAttt                                                                                                            <  2:946806/139‑1 (MQ=255)
      tCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCt                                                                                                          <  1:186390/139‑1 (MQ=255)
              tgAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCAt                                                                                                                  <  1:1087511/123‑1 (MQ=255)
              tgAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCAt                                                                                                                  >  2:1087511/1‑123 (MQ=255)
                    gCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGttt                                                                                            <  1:1473786/139‑1 (MQ=255)
                    gCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGttt                                                                                            >  2:40266/1‑139 (MQ=255)
                          taagagacagggCTGAGGTGCTGACTCTGGCGTCGTAACCTGTCACATGTTATTCGCATGCAGTCATTCATCGACTCATGCCTTGCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtggct                                                                                      <  2:1090297/129‑1 (MQ=255)
                            cagTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtggctgg                                                                                    <  2:384296/139‑1 (MQ=255)
                                           gTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCa                                                                                                                   <  2:202577/93‑1 (MQ=255)
                                           gTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCa                                                                                                                   >  1:202577/1‑93 (MQ=255)
                                           gTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATgcg                                                                     >  1:124808/1‑139 (MQ=255)
                                               tGGCTCTGGCGTCGTAACGTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGcctgtcgctta                                                                 >  1:1090297/1‑129 (MQ=255)
                                                   tctGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTAc                                                             >  2:437301/1‑139 (MQ=255)
                                                     tGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTActc                                                           <  1:1486368/139‑1 (MQ=255)
                                                      ggCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTActct                                                          >  2:749510/1‑139 (MQ=255)
                                                           cgtAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCGGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaaa                                                     <  1:940086/139‑1 (MQ=255)
                                                             tAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACg                                                   <  1:437301/139‑1 (MQ=255)
                                                                            tATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTaa                                    <  2:124808/139‑1 (MQ=255)
                                                                                ggCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATgc                                                                      <  1:1233371/101‑1 (MQ=255)
                                                                                ggCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATgc                                                                      >  2:1233371/1‑101 (MQ=255)
                                                                                  cATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACg                                                   <  2:844060/118‑1 (MQ=255)
                                                                                  cATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACg                                                   >  1:844060/1‑118 (MQ=255)
                                                                                   aTGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAaga                             <  2:1153250/139‑1 (MQ=255)
                                                                                     gCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAgaga                           >  1:429753/1‑139 (MQ=255)
                                                                                      cAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAAt                                                                                                              <  2:1433699/55‑1 (MQ=255)
                                                                                      cAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAAt                                                                                                              >  1:1433699/1‑55 (MQ=255)
                                                                                      cAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                          <  1:1057893/139‑1 (MQ=255)
                                                                                      cAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                          <  2:429753/139‑1 (MQ=255)
                                                                                        gTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATa                        >  1:1470204/1‑139 (MQ=255)
                                                                                           attcatCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCt                                                                                                          <  2:884969/54‑1 (MQ=255)
                                                                                           attcatCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCt                                                                                                          >  1:884969/1‑54 (MQ=255)
                                                                                           attcatCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGtt                                                               <  1:206873/97‑1 (MQ=255)
                                                                                           attcatCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGtt                                                               >  2:206873/1‑97 (MQ=255)
                                                                                           attcatCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGt                                                                <  1:268529/96‑1 (MQ=255)
                                                                                           attcatCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGt                                                                >  1:1426060/1‑96 (MQ=255)
                                                                                           attcatCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGt                                                                <  2:1426060/96‑1 (MQ=255)
                                                                                           attcatCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGt                                                                >  2:268529/1‑96 (MQ=255)
                                                                                           attcatCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCt                     >  2:451043/1‑139 (MQ=255)
                                                                                                tCGACTCAGGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCa                <  2:1625182/139‑1 (MQ=255)
                                                                                                   aCTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGt                                                                <  1:1526686/88‑1 (MQ=255)
                                                                                                   aCTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGt                                                                >  2:1526686/1‑88 (MQ=255)
                                                                                                     tCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCatgatg           >  2:735043/1‑139 (MQ=255)
                                                                                                      cATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaa                                                      >  1:157194/1‑95 (MQ=255)
                                                                                                      cATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaa                                                      <  2:157194/95‑1 (MQ=255)
                                                                                                              tCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATcc  >  2:161168/1‑139 (MQ=255)
                                                                                                                     tatCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGgcggcg                                                                               <  1:1598887/55‑1 (MQ=255)
                                                                                                                     tatCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGgcggcg                                                                               >  2:1598887/1‑55 (MQ=255)
                                                                                                                             |                                                                                                                           
ATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCC  >  NC_000913/28144‑28392

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: