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.95.2% 157.8 / 3.7 42intergenic (‑14/+315)ykfA/perRGTPase family protein/LysR substrate‑binding domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base . (16/24);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGATTTTTT  >  NZ_CP009273/264596‑264769
                                                                                         |                                                                                    
gAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGa                                                                                      >  2:1454143/1‑90 (MQ=255)
       ctctGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGaaaataat                                                                               >  2:584627/1‑86 (MQ=255)
           gATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATaaaa                                                                          >  2:676599/1‑90 (MQ=255)
            aTACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATaaaaa                                                                         <  2:934703/90‑1 (MQ=255)
             tACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATaaaaaa                                                                        >  1:1363888/1‑90 (MQ=255)
                 cATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACgg                                                                     >  2:830640/1‑90 (MQ=255)
                   ttGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTac                                                                  <  2:1573962/90‑2 (MQ=255)
                   ttGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAg                                                                  >  1:811677/1‑90 (MQ=255)
                   ttGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAg                                                                  >  2:1528499/1‑90 (MQ=255)
                     gCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGaa                                                                >  1:1247529/1‑90 (MQ=255)
                                gCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTc                                                     <  2:756396/90‑1 (MQ=255)
                                gCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTc                                                     <  2:88683/90‑1 (MQ=255)
                                gCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTc                                                     >  2:213966/1‑90 (MQ=255)
                                       aTTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGtt                                              <  1:259864/90‑1 (MQ=255)
                                        ttGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTa                                             >  2:1046389/1‑90 (MQ=255)
                                           cTGAAACGACTTCAAACCTTCAGAATTATTCATCACTTCTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAAtt                                          <  1:751861/90‑1 (MQ=255)
                                           cTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAAtt                                          <  2:581613/90‑1 (MQ=255)
                                           cTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAAtt                                          <  2:811677/90‑1 (MQ=255)
                                               aaCGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGt                                      <  2:1247529/90‑1 (MQ=255)
                                               aaCGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGt                                      <  2:126077/90‑1 (MQ=255)
                                               aaCGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGt                                      <  2:113941/90‑1 (MQ=255)
                                               aaCGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGt                                      <  1:1127617/90‑1 (MQ=255)
                                                aCGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTGCGCCGTTAATTGCGTa                                     <  2:1134097/90‑1 (MQ=255)
                                                     ttCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTc                                <  1:1399946/90‑1 (MQ=255)
                                                       cAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCag                              <  1:434567/90‑1 (MQ=255)
                                                       cAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCag                              <  1:1428152/90‑1 (MQ=255)
                                                         aaCCTTCAGAATTATTCATCACTCCTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCagag                            <  1:1163366/90‑1 (MQ=255)
                                                         aaCCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCagag                            <  2:1524034/90‑1 (MQ=255)
                                                             ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGt                                      >  2:1301322/1‑76 (MQ=255)
                                                             ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGt                                      <  1:1301322/76‑1 (MQ=255)
                                                             ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  2:1130005/90‑1 (MQ=255)
                                                             ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  1:1046389/90‑1 (MQ=255)
                                                             ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  2:1081707/90‑1 (MQ=255)
                                                             ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  1:1568893/90‑1 (MQ=255)
                                                              tCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgatg                       >  1:1066482/1‑90 (MQ=255)
                                                                 gAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGata                    >  2:77739/1‑90 (MQ=255)
                                                                     tattCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGatatata                >  1:1018140/1‑90 (MQ=255)
                                                                     tattCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGatatata                >  1:277302/1‑90 (MQ=255)
                                                                          atcaCTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGa           <  2:1399865/90‑1 (MQ=255)
                                                                           tcaCTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGaa          >  2:774471/1‑90 (MQ=255)
                                                                           tcaCTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGaa          >  2:1367131/1‑90 (MQ=255)
                                                                                actCCTCTGAAAAATAATAAAAAGCGGTAGGATCGTGAGATTCCGCCGTAAATTGCGGATGTTCAGAGGGATGAGATATATCTGAAGAtt      >  2:204600/1‑90 (MQ=25)
                                                                                actCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGAttt     >  1:846822/1‑90 (MQ=255)
                                                                                   cctcTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGAtttttt  >  2:1399261/1‑90 (MQ=255)
                                                                                         |                                                                                    
GAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGATTTTTT  >  NZ_CP009273/264596‑264769

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