Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 264,689 (A)5→4 intergenic (‑18/+311) ykfA ← / ← perR GTPase family protein/LysR substrate‑binding domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273264,6850A.98.8% 371.9 / ‑2.6 85intergenic (‑14/+315)ykfA/perRGTPase family protein/LysR substrate‑binding domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (37/47);  total (37/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGATTTTTT  >  NZ_CP009273/264603‑264769
                                                                                  |                                                                                    
ctctGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGaaaataat                                                                               >  1:2469392/1‑86 (MQ=255)
ctctGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGaaaataat                                                                               >  2:2206294/1‑86 (MQ=255)
  ctGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAtaataa                                                                            <  1:2842864/90‑1 (MQ=255)
     aTACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATaaaaa                                                                         <  1:449527/90‑1 (MQ=255)
     aTACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATaaaaa                                                                         <  1:522146/90‑1 (MQ=255)
     aTACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATaaaaa                                                                         <  1:760939/90‑1 (MQ=255)
     aTACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATaaaaa                                                                         <  1:2540798/90‑1 (MQ=255)
     aTACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATaaaaa                                                                         <  2:2469392/90‑1 (MQ=255)
     aTACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATaaaaa                                                                         <  2:76592/90‑1 (MQ=255)
       aCCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAAc                                                                       <  2:2459307/90‑1 (MQ=255)
          cATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGt                                                                    >  1:2777637/1‑90 (MQ=255)
          cATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGt                                                                    >  2:822127/1‑90 (MQ=255)
          cATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGt                                                                    >  2:779086/1‑90 (MQ=255)
          cATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGt                                                                    >  2:186896/1‑90 (MQ=255)
           aTTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGGATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTa                                                                   <  1:573666/90‑1 (MQ=255)
            ttGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAg                                                                  >  2:1773170/1‑90 (MQ=255)
            ttGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAg                                                                  >  2:1350584/1‑90 (MQ=255)
            ttGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAg                                                                  >  1:2237714/1‑90 (MQ=255)
            ttGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAg                                                                  >  2:424190/1‑90 (MQ=255)
              gCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGaa                                                                >  1:1986008/1‑90 (MQ=255)
              gCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGaa                                                                >  1:1493343/1‑90 (MQ=255)
              gCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGaa                                                                >  2:1397504/1‑90 (MQ=255)
                 gCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAGAATAAAAAACGGTAGAATCg                                                             >  2:2849187/1‑90 (MQ=255)
                         gCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTc                                                     <  2:411180/90‑1 (MQ=255)
                         gCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTc                                                     <  2:640034/90‑1 (MQ=255)
                           aaGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTccg                                                   >  2:1365586/1‑90 (MQ=255)
                             gCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTccgcc                                                 <  1:1678582/90‑1 (MQ=255)
                              cGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTccgccg                                                >  2:522172/1‑90 (MQ=255)
                                aTTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGtt                                              >  2:724081/1‑90 (MQ=255)
                                aTTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGtt                                              <  2:1115342/90‑1 (MQ=255)
                                aTTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGtt                                              >  1:1167229/1‑90 (MQ=255)
                                 ttGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTa                                             <  2:762852/90‑1 (MQ=255)
                                 ttGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTa                                             <  2:2180432/90‑1 (MQ=255)
                                    cTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAAtt                                          <  1:779086/90‑1 (MQ=255)
                                    cTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAAtt                                          <  2:931845/90‑1 (MQ=255)
                                    cTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAAtt                                          <  2:1493343/90‑1 (MQ=255)
                                        aaCGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGt                                      >  2:935412/1‑90 (MQ=255)
                                        aaCGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGt                                      <  1:2071863/90‑1 (MQ=255)
                                        aaCGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGt                                      <  2:2759934/90‑1 (MQ=255)
                                        aaCGACCTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGt                                      <  1:1523226/90‑1 (MQ=255)
                                         aCGACTTCAAACCTTCAGAATTATTCATCACTACTCGTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTa                                     <  1:561165/90‑1 (MQ=255)
                                         aCGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTa                                     <  1:1350584/90‑1 (MQ=255)
                                         aCGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTa                                     <  2:1986008/90‑1 (MQ=255)
                                         aCGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTa                                     <  1:2451411/90‑1 (MQ=255)
                                              ttCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTc                                <  2:1453920/90‑1 (MQ=255)
                                              ttCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTc                                <  2:2706296/90‑1 (MQ=255)
                                              ttCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTc                                <  1:2219549/90‑1 (MQ=255)
                                                cAAACCTTCAGAATTATTCATCACTCCTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCag                              <  1:2565322/90‑1 (MQ=255)
                                                 aaaCCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCaga                             <  2:2199363/90‑1 (MQ=255)
                                                  aaCCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCagag                            <  1:1773170/90‑1 (MQ=255)
                                                      ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCagag                            >  2:2154963/1‑86 (MQ=255)
                                                      ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCagag                            <  1:2154963/86‑1 (MQ=255)
                                                      ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  1:522172/90‑1 (MQ=255)
                                                      ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  2:2281798/90‑1 (MQ=255)
                                                      ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  1:186896/90‑1 (MQ=255)
                                                      ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        >  2:2904483/1‑90 (MQ=255)
                                                      ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  2:375406/90‑1 (MQ=255)
                                                      ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  1:1675676/90‑1 (MQ=255)
                                                      ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  2:443509/90‑1 (MQ=255)
                                                      ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  1:1434493/90‑1 (MQ=255)
                                                      ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  1:935412/90‑1 (MQ=255)
                                                      ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  1:424190/90‑1 (MQ=255)
                                                      ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  1:750491/90‑1 (MQ=255)
                                                      ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  1:988859/90‑1 (MQ=255)
                                                      ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  1:2226904/90‑1 (MQ=255)
                                                             ttattCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGatatat                 <  1:2571722/90‑1 (MQ=255)
                                                               attCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGatatatat               >  2:755828/1‑90 (MQ=255)
                                                                    tcaCTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGaa          >  1:2068637/1‑90 (MQ=255)
                                                                    tcaCTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGaa          >  1:1286282/1‑90 (MQ=255)
                                                                    tcaCTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGaa          >  1:593810/1‑90 (MQ=255)
                                                                    tcaCTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGaa          >  2:1847394/1‑90 (MQ=255)
                                                                    tcaCTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGaa          >  1:1533065/1‑90 (MQ=255)
                                                                    tcaCTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGaa          >  1:1623293/1‑90 (MQ=255)
                                                                    tcaCTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGaa          >  1:1034275/1‑90 (MQ=255)
                                                                    tcaCTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGaa          >  1:2683573/1‑90 (MQ=255)
                                                                      actactCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGa        >  2:464243/1‑90 (MQ=255)
                                                                       ctactCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGAt       >  2:399480/1‑90 (MQ=255)
                                                                        tactCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGAtt      <  1:1375944/90‑1 (MQ=255)
                                                                         actCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGAttt     >  2:290737/1‑90 (MQ=255)
                                                                         actCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGAttt     >  1:1363672/1‑90 (MQ=255)
                                                                         actCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGAttt     >  1:1892332/1‑90 (MQ=255)
                                                                         actCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGAttt     >  2:1546035/1‑90 (MQ=255)
                                                                            cctcTGAAAATAATAAAAAACGGTAGAATCGTGATATTCCGCCGTTa                                             >  2:1921933/1‑47 (MQ=25)
                                                                            cctcTGAAAATAATAAAAAACGGTAGAATCGTGATATTCCGCCGTTa                                             <  1:1921933/47‑1 (MQ=25)
                                                                            cctcTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  2:257286/68‑1 (MQ=255)
                                                                            cctcTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        >  1:257286/1‑68 (MQ=255)
                                                                            cctcTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGAtttttt  >  1:300799/1‑90 (MQ=255)
                                                                                  |                                                                                    
CTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGATTTTTT  >  NZ_CP009273/264603‑264769

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