Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 49,791 G→A intergenic (+160/‑32) kefC → / → folA glutathione‑regulated potassium‑efflux system protein KefC/type 3 dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,7910GA97.9% 178.3 / ‑4.3 47intergenic (+160/‑32)kefC/folAglutathione‑regulated potassium‑efflux system protein KefC/type 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (27/19);  total (27/20)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAA  >  NZ_CP009273/49707‑49873
                                                                                     |                                                                                  
agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCg                                                                               >  1:856413/1‑90 (MQ=255)
agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCg                                                                               >  1:795375/1‑90 (MQ=255)
   gCGTGGACTCGCCAGCAGAATATAAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGAc                                                                             >  1:379124/1‑90 (MQ=255)
           tCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAAttttttt                                                                    >  2:263644/1‑90 (MQ=255)
            cGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAAtttttttt                                                                   >  1:846716/1‑90 (MQ=255)
             gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTa                                                                  >  2:155767/1‑90 (MQ=255)
                   agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGa                                                            <  1:807832/90‑1 (MQ=255)
                   agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGa                                                            <  2:348973/90‑1 (MQ=255)
                   agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGa                                                            <  1:562796/90‑1 (MQ=255)
                    gAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGaa                                                           <  1:470383/90‑1 (MQ=255)
                    gAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGaa                                                           <  2:585691/90‑1 (MQ=255)
                               ttttCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGa                                                 <  1:329136/90‑1 (MQ=255)
                               ttttCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGa                                                 >  2:21922/1‑90 (MQ=255)
                                 ttCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATc                                               >  2:525805/1‑90 (MQ=255)
                                 ttCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATc                                               >  1:178239/1‑90 (MQ=255)
                                 ttCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATc                                               <  1:21922/90‑1 (MQ=255)
                                 ttCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATc                                               >  1:2728/1‑90 (MQ=255)
                                   ccTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGaa                                                           >  1:646129/1‑76 (MQ=255)
                                   ccTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGaa                                                           <  2:646129/76‑1 (MQ=255)
                                        aCATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGa                                        >  2:424510/1‑90 (MQ=255)
                                            catcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgc                                    <  1:932382/90‑1 (MQ=255)
                                             atcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcg                                   >  1:654966/1‑90 (MQ=255)
                                             atcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcg                                   >  2:928182/1‑90 (MQ=255)
                                              tcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcgg                                  <  2:856413/90‑1 (MQ=255)
                                                cTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcggcg                                <  1:155767/90‑1 (MQ=255)
                                                  cGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGtt                              <  1:95258/90‑1 (MQ=255)
                                                  cGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGtt                              >  2:949896/1‑90 (MQ=255)
                                                  cGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGtt                              >  1:644444/1‑90 (MQ=255)
                                                  cGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGtt                              >  1:591004/1‑90 (MQ=255)
                                                      ccAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCg                          <  1:794405/90‑1 (MQ=255)
                                                        aGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGt                        >  1:844697/1‑90 (MQ=255)
                                                        aGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGt                        >  1:248043/1‑90 (MQ=255)
                                                        aGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGt                        >  2:208906/1‑90 (MQ=255)
                                                        aGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGt                        >  1:456841/1‑90 (MQ=255)
                                                        aGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGt                        >  2:448566/1‑90 (MQ=255)
                                                        aGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGt                        >  1:429655/1‑90 (MQ=255)
                                                            gacgacGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGGGtt                              <  1:27781/80‑1 (MQ=255)
                                                            gacgacGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGGGtt                              >  2:27781/1‑80 (MQ=255)
                                                                   gTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTa             >  2:268934/1‑90 (MQ=255)
                                                                     ttACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATc           <  1:418271/90‑1 (MQ=255)
                                                                     ttACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATc           <  2:429292/90‑1 (MQ=255)
                                                                         gCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCa       <  1:239482/90‑1 (MQ=255)
                                                                            ttACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATgg    >  2:273630/1‑90 (MQ=255)
                                                                            ttACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATgg    >  1:701403/1‑90 (MQ=255)
                                                                              aCGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGaa  <  1:263644/90‑1 (MQ=255)
                                                                              aCGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGaa  <  1:831784/90‑1 (MQ=255)
                                                                              aCGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGaa  <  1:945130/90‑1 (MQ=255)
                                                                                     |                                                                                  
AGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAA  >  NZ_CP009273/49707‑49873

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