Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,624,386 T→C E522G (GAA→GGA)  yhiJ ← YhiJ/YhiL family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,624,3860TC100.0% 128.5 / NA 39E522G (GAA→GGA) yhiJYhiJ/YhiL family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/12);  total (27/12)

AACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTTCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTC  >  NZ_CP009273/3624301‑3624470
                                                                                     |                                                                                    
aaCATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATg                                                                                  <  2:423251/90‑1 (MQ=255)
aaCATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATg                                                                                  <  1:335119/90‑1 (MQ=255)
 aCATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGc                                                                                 <  1:649/90‑1 (MQ=255)
    ttAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTgg                                                                              >  1:309731/1‑90 (MQ=255)
      aaTATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTa                                                                            <  1:530021/90‑1 (MQ=255)
        tatTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAg                                                                          >  1:68866/1‑90 (MQ=255)
            caagttgtCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACtt                                                                      >  2:227816/4‑90 (MQ=255)
            cATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACtt                                                                      >  2:452440/1‑90 (MQ=255)
            cATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACtt                                                                      >  1:376358/1‑90 (MQ=255)
            cATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACtt                                                                      >  2:549762/1‑90 (MQ=255)
              tgttgtCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTcc                                                                    >  1:94930/1‑90 (MQ=255)
                 tgtCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAgg                                                                 >  1:567715/1‑90 (MQ=255)
                    cTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAg                                                              <  1:549762/90‑1 (MQ=255)
                         gTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGtt                                                         >  1:63790/1‑90 (MQ=255)
                          ttCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTa                                                        <  2:376358/90‑1 (MQ=255)
                           tCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTAt                                                       >  1:124039/1‑90 (MQ=255)
                           tCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTAt                                                       >  2:135697/1‑90 (MQ=255)
                            cAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTa                                                                            >  1:512090/1‑68 (MQ=255)
                            cAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTa                                                                            <  2:512090/68‑1 (MQ=255)
                            cAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGt                                                          <  1:287824/86‑1 (MQ=255)
                            cAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGt                                                          >  2:287824/1‑86 (MQ=255)
                            cAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTAtc                                                      >  1:49401/1‑90 (MQ=255)
                            cAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTAtc                                                      >  1:330451/1‑90 (MQ=255)
                              gagCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAg                                                                          <  1:363378/68‑1 (MQ=255)
                              gagCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAg                                                                          >  2:363378/1‑68 (MQ=255)
                               agCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCa                                                   <  1:315606/90‑1 (MQ=255)
                                   gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGcc                                               >  2:290351/1‑90 (MQ=255)
                                   gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGcc                                               >  1:150223/1‑90 (MQ=255)
                                   gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGcc                                               >  2:416395/1‑90 (MQ=255)
                                          tgtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAggg                                                                >  2:327099/1‑66 (MQ=255)
                                          tgtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAggg                                                                <  1:327099/66‑1 (MQ=255)
                                          tgtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGTGAGTGGTTATCTCACGCCCTGCAcc                                        >  1:545871/1‑90 (MQ=255)
                                          tgtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAcc                                        >  2:552029/1‑90 (MQ=255)
                                           gtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCa                                       >  2:527306/1‑90 (MQ=255)
                                             gcgcCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAgg                                     >  2:115895/1‑90 (MQ=255)
                                             gcgcCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAgg                                     >  1:579654/1‑90 (MQ=255)
                                                                 tAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTatgca                 <  1:461500/90‑1 (MQ=255)
                                                                       cATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGa           >  1:398143/1‑90 (MQ=255)
                                                                        aTTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGaa          >  1:193280/1‑90 (MQ=255)
                                                                                tttttcCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACAttcttc  >  1:586553/1‑90 (MQ=255)
                                                                                     |                                                                                    
AACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTTCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTC  >  NZ_CP009273/3624301‑3624470

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