Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,591,603 G→A *119* (TAG→TAA yffB → putative ArsC family reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,591,6030GA100.0% 105.7 / NA 33*119* (TAG→TAAyffBputative ArsC family reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/17);  total (16/17)

CCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGA  >  NC_000913/2591467‑2591728
                                                                                                                                        |                                                                                                                             
ccGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATc                                                                                                                             >  2:865268/1‑139 (MQ=255)
           tctGCGGCGACATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGcc                                                                                                                  <  1:115010/139‑1 (MQ=255)
              gcggcgGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCgg                                                                                                               <  2:515087/139‑1 (MQ=255)
                              cTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGATCTGAcaca                                                                                               >  2:689470/1‑139 (MQ=255)
                              cTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGAcaca                                                                                               <  1:148618/139‑1 (MQ=255)
                                gAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACacaac                                                                                             >  2:863313/1‑139 (MQ=255)
                                                tCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGccc                                                                             >  2:710633/1‑139 (MQ=255)
                                                              ctctGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCtg                                                               >  1:878351/1‑139 (MQ=255)
                                                                   cGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGAcaca                                                                                               <  2:597064/102‑1 (MQ=255)
                                                                   cGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGAcaca                                                                                               >  1:597064/1‑102 (MQ=255)
                                                                     tGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAATTATCAGCAATTTTTCCATGAGGTGTCATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTCAGTCCTGATGATGc                                                        <  1:689470/139‑1 (MQ=255)
                                                                        ccGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGacac                                                                                                <  2:11449/96‑1 (MQ=255)
                                                                        ccGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGacac                                                                                                >  1:11449/1‑96 (MQ=255)
                                                                           gTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATg                                                  <  2:477236/139‑1 (MQ=255)
                                                                             aaGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTcc                                                                         <  1:532284/114‑1 (MQ=255)
                                                                             aaGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTcc                                                                         >  2:532284/1‑114 (MQ=255)
                                                                               gCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAg                                              <  1:863313/139‑1 (MQ=255)
                                                                                     gctgctGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACacaac                                                                                             <  1:960651/86‑1 (MQ=255)
                                                                                     gctgctGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACacaac                                                                                             >  2:960651/1‑86 (MQ=255)
                                                                                         ctgGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGttga                                    >  2:1003772/1‑139 (MQ=255)
                                                                                          tgGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTAtt                                                                                     <  2:1097803/89‑1 (MQ=255)
                                                                                          tgGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTAtt                                                                                     >  1:1097803/1‑89 (MQ=255)
                                                                                                tCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCgccgcc                                                                              >  1:1058943/1‑90 (MQ=255)
                                                                                                tCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCgccgcc                                                                              <  2:1058943/90‑1 (MQ=255)
                                                                                                        tCCAGTTATCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAgg                     <  2:878351/139‑1 (MQ=255)
                                                                                                                tCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcg                                                                                   <  2:1084648/69‑1 (MQ=255)
                                                                                                                tCAGCAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTcg                                                                                   >  1:1084648/1‑69 (MQ=255)
                                                                                                                    cAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCgccgcc                                                                              >  2:1029370/1‑70 (MQ=255)
                                                                                                                    cAATTTTTCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCgccgcc                                                                              <  1:1029370/70‑1 (MQ=255)
                                                                                                                       tttttCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCttt                                         <  2:1098862/104‑1 (MQ=255)
                                                                                                                       tttttCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCttt                                         >  1:1098862/1‑104 (MQ=255)
                                                                                                                       tttttCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTAACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCg      >  2:802867/1‑139 (MQ=255)
                                                                                                                           tCCATGAGGTGTAATCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGa  <  1:710633/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                             
CCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGA  >  NC_000913/2591467‑2591728

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