Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,554,121 A→C N169T (AAC→ACC)  BW25113_RS12740 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,554,1210AC100.0% 107.4 / NA 31N169T (AAC→ACC) BW25113_RS12740hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (18/13);  total (18/13)

CTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAACACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCA  >  NZ_CP009273/2554035‑2554203
                                                                                      |                                                                                  
ctacGCTGGTACTGGCAAAAACGCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAccac                                                                                 <  1:179059/87‑1 (MQ=255)
cTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAccac                                                                                 <  2:125944/90‑1 (MQ=255)
cTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAccac                                                                                 <  2:263377/90‑1 (MQ=255)
cTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAccac                                                                                 <  2:32697/90‑1 (MQ=255)
cTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAccac                                                                                 <  2:364040/90‑1 (MQ=255)
  tCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACgg                                                                               >  1:588962/1‑90 (MQ=255)
                       cacaGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTc                                                          >  2:61406/1‑90 (MQ=255)
                               ggATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAAc                                                  <  1:267781/90‑1 (MQ=255)
                                   ggTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGGTCGTCTGGTCAACCACt                                              >  1:300978/1‑90 (MQ=255)
                                   ggTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACt                                              >  1:25632/1‑90 (MQ=255)
                                         gggTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAAc                                        >  2:84672/1‑90 (MQ=255)
                                              ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAAcc                                                 >  2:60429/1‑76 (MQ=255)
                                              ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAAcc                                                 <  1:60429/76‑1 (MQ=255)
                                              ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt                                   >  2:368108/1‑90 (MQ=255)
                                              ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt                                   >  2:530938/1‑90 (MQ=255)
                                              ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt                                   >  2:688282/1‑90 (MQ=255)
                                              ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt                                   >  1:145811/1‑90 (MQ=255)
                                              ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt                                   >  1:138863/1‑90 (MQ=255)
                                              ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCCCGACCGCCTGGTCAACCCCTCCAAACttgtt                                   >  1:227166/1‑90 (MQ=255)
                                                   gggAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgttg                                  >  1:47543/1‑86 (MQ=255)
                                                   gggAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgttg                                  <  2:47543/86‑1 (MQ=255)
                                                    ggAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGt                             >  2:67497/1‑90 (MQ=255)
                                                         aaTCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGa                        >  1:294383/1‑90 (MQ=255)
                                                          aTCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGaa                       <  1:490421/90‑1 (MQ=255)
                                                          aTCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGaa                       <  1:415262/90‑1 (MQ=255)
                                                           tCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAg                      <  2:317027/90‑1 (MQ=255)
                                                           tCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAg                      <  2:414974/90‑1 (MQ=255)
                                                            cGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAgg                     >  2:257888/1‑90 (MQ=255)
                                                            cGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAgg                     <  1:521954/90‑1 (MQ=255)
                                                             gACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGc                    >  2:372916/1‑90 (MQ=255)
                                                               cAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCta                  <  1:530938/90‑1 (MQ=255)
                                                               cAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCta                  <  1:81931/90‑1 (MQ=255)
                                                               cAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCta                  <  1:612769/90‑1 (MQ=255)
                                                                    gAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATAc             <  1:560078/90‑1 (MQ=255)
                                                                      gCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCa           <  1:257888/90‑1 (MQ=255)
                                                                               tCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCa  >  1:608960/1‑90 (MQ=255)
                                                                                      |                                                                                  
CTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAACACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCA  >  NZ_CP009273/2554035‑2554203

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