Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,579,109 A→C G545G (GGT→GGG ydeP ← acid resistance putative oxidoreductase YdeP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,579,1090AC95.8% 229.8 / 1.6 72G545G (GGT→GGGydePacid resistance putative oxidoreductase YdeP
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base C (37/32);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-01

CAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGC  >  NZ_CP009273/1579025‑1579193
                                                                                    |                                                                                    
caGATACTCCCAGTCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCcggc                                                                                 <  1:1372443/90‑1 (MQ=255)
caGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCcggc                                                                                 <  1:360032/90‑1 (MQ=255)
caGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCcggc                                                                                 <  1:103644/90‑1 (MQ=255)
caGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCcggc                                                                                 <  2:1196997/90‑1 (MQ=255)
 aGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCcggcg                                                                                <  1:669846/90‑1 (MQ=255)
  gATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCcggcgg                                                                               <  2:620199/90‑1 (MQ=255)
        cccAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTa                                                                         >  2:430215/1‑90 (MQ=255)
            ggCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAacac                                                                     >  2:233753/1‑90 (MQ=255)
              cTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAAcacgc                                                                   <  1:1493923/90‑1 (MQ=255)
              cTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAAcacgc                                                                   >  1:1498984/1‑90 (MQ=255)
                aCCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAAcacgcca                                                                 >  2:133667/1‑90 (MQ=255)
                aCCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAAcacgcca                                                                 >  2:170525/1‑90 (MQ=255)
                aCCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAAcacgcca                                                                 >  2:467316/1‑90 (MQ=255)
                    cGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCAcgcg                                                             <  1:1553280/90‑1 (MQ=255)
                     gCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTCCCCCGGCGGGTTTTAACACGCCACGCGa                                                            <  1:1348218/90‑1 (MQ=255)
                      ctctGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAg                                                           >  2:135820/1‑90 (MQ=255)
                         tGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCa                                                        <  1:1418816/90‑1 (MQ=255)
                                gTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCcggc                                                                                 <  1:732955/58‑1 (MQ=255)
                                gTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCcggc                                                                                 >  2:732955/1‑58 (MQ=255)
                                gTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAatca                                                 >  1:1042816/1‑90 (MQ=255)
                                gTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAatca                                                 >  2:473355/1‑90 (MQ=255)
                                gTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAatca                                                 >  1:1473113/1‑90 (MQ=255)
                                gTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAatca                                                 >  1:429662/1‑90 (MQ=255)
                                gTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAatca                                                 >  1:1269643/1‑90 (MQ=255)
                                gTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGGTTTAACAACCCCCCCGGGGGATGAaaca                                                 >  1:1120714/1‑90 (MQ=255)
                                  gctgcCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAatcatc                                               >  2:578910/1‑90 (MQ=255)
                                  gctgcCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAatcatc                                               >  2:1252709/1‑90 (MQ=255)
                                  gctgcCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAatcatc                                               <  2:1269643/90‑1 (MQ=255)
                                  gctgcCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAatcatc                                               >  1:1280370/1‑90 (MQ=255)
                                  gctgcCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAatcatc                                               >  2:1487295/1‑90 (MQ=255)
                                    tgcctgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGa                                             <  2:1498984/90‑1 (MQ=255)
                                    tgcctgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGa                                             <  2:269089/90‑1 (MQ=255)
                                    tgcctgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGa                                             <  2:678296/90‑1 (MQ=255)
                                     gcctgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCcggc                                                                                 <  1:652418/53‑1 (MQ=255)
                                     gcctgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCcggc                                                                                 >  2:652418/1‑53 (MQ=255)
                                     gcctgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGAc                                            >  2:100849/1‑90 (MQ=255)
                                       ctgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACAt                                          >  1:191216/1‑90 (MQ=255)
                                         gcgcGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACAtt                                         >  2:450439/1‑89 (MQ=255)
                                         gcgcGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACAtt                                         <  1:450439/89‑1 (MQ=255)
                                         gcgcGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTg                                        >  1:1382465/1‑90 (MQ=255)
                                         gcgcGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTg                                        >  2:647340/1‑90 (MQ=255)
                                         gcgcGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTg                                        <  1:198300/90‑1 (MQ=255)
                                         gcgcGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTg                                        >  2:366627/1‑90 (MQ=255)
                                         gcgcGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTg                                        <  2:1473113/90‑1 (MQ=255)
                                         gcgcGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTg                                        >  1:424022/1‑90 (MQ=255)
                                          cgcgATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTg                                        <  2:701387/89‑1 (MQ=255)
                                          cgcgATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTg                                        >  1:701387/1‑89 (MQ=255)
                                           gcgATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGaa                                      <  2:991963/90‑1 (MQ=255)
                                           gcgATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGaa                                      <  2:1382465/90‑1 (MQ=255)
                                             gATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATc                                    >  1:665549/1‑90 (MQ=255)
                                               ttCCCGCGACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCt                                  >  1:440374/1‑90 (MQ=255)
                                                   cgcgACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAAc                              <  2:1550133/90‑1 (MQ=255)
                                                   cgcgACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAAc                              >  2:1075066/1‑90 (MQ=255)
                                                   cgcgACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAAc                              <  1:1552941/90‑1 (MQ=255)
                                                   cgcgACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAAc                              <  2:429662/90‑1 (MQ=255)
                                                    gcgACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACg                             >  1:454862/1‑90 (MQ=255)
                                                    gcgACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACg                             >  2:312656/1‑90 (MQ=255)
                                                      gACCACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGt                           >  1:656894/1‑90 (MQ=255)
                                                      gACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGt                           >  2:1007761/1‑90 (MQ=255)
                                                        ccACTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTa                         >  1:299602/1‑90 (MQ=255)
                                                          aCTGCACACTCTGTTTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTAcc                       <  2:383575/90‑1 (MQ=255)
                                                          aCTGCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTAcc                       <  1:366627/90‑1 (MQ=255)
                                                             gCACACTCTGATTTCAGCATTACCCCGGCGGGTTTTAACCCGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGGTACCGcc                    >  2:1499886/1‑90 (MQ=255)
                                                                  ctctGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATTAATCATCGACATTGAATCCTCAACGGTTACCGCCTgcgc               <  2:1422445/90‑1 (MQ=255)
                                                                  ctctGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTgcgc               >  2:834410/1‑90 (MQ=255)
                                                                    ctGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAGTCCTCAACGGTTACCGCCTGCGCAc             <  1:402409/90‑1 (MQ=255)
                                                                    ctGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCAc             <  1:170525/90‑1 (MQ=255)
                                                                    ctGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCAc             <  2:771037/90‑1 (MQ=255)
                                                                    ctGATTTCAGCATTACACCGGCGGGGTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCCGCGCAc             >  1:77363/1‑90 (MQ=255)
                                                                     tGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCAcc            >  2:410737/1‑90 (MQ=255)
                                                                            aGCATTACCCCGGCGGGATTTGATCCGCCACGCGAGGCTTGAATCATGGAGATTGAATCCTCAACGGGTACCTCCTACGCACCGgttttt     >  1:857951/1‑90 (MQ=255)
                                                                               aTTACCCCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGc  <  1:410737/90‑1 (MQ=255)
                                                                                    |                                                                                    
CAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGC  >  NZ_CP009273/1579025‑1579193

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