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.87.5% 111.4 / 22.7 41intergenic (‑14/+315)ykfA/perRGTPase family protein/LysR substrate‑binding domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base . (14/21);  total (17/24)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGATTTTTT  >  NZ_CP009273/264598‑264769
                                                                                       |                                                                                    
aGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGaaa                                                                                    <  2:318706/90‑1 (MQ=255)
          aTACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATaaaaa                                                                         <  2:511869/90‑1 (MQ=255)
               cATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGt                                                                    >  1:825452/1‑90 (MQ=255)
                 ttGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTa                                                                   >  2:796263/1‑90 (MQ=255)
                 ttGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAg                                                                  >  2:791536/1‑90 (MQ=255)
                   gCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGaa                                                                >  2:805638/1‑90 (MQ=255)
                      gCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATc                                                              >  2:334469/1‑90 (MQ=255)
                         gATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTga                                                          >  1:445188/1‑90 (MQ=255)
                                      ttGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGtt                                              >  1:194597/1‑90 (MQ=255)
                                        gCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAAt                                           >  1:381251/1‑90 (MQ=255)
                                          tGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAAtt                                          <  2:194597/90‑1 (MQ=255)
                                             aaCGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGt                                      <  2:106836/90‑1 (MQ=255)
                                              aCGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTa                                     <  1:311350/90‑1 (MQ=255)
                                                  cTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGtt                                 >  1:778425/1‑90 (MQ=255)
                                                   ttCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTc                                <  2:141486/90‑1 (MQ=255)
                                                       aaCCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCagag                            <  2:494105/90‑1 (MQ=255)
                                                       aaCCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCagag                            <  1:791536/90‑1 (MQ=255)
                                                        aCCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGt                           <  2:565953/90‑1 (MQ=255)
                                                        aCCTTCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGt                           <  2:431654/90‑1 (MQ=255)
                                                         ccTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGt                           <  1:796263/90‑1 (MQ=255)
                                                          attCAGAATTATTCCTCACTCCTCCTCTTAAAATAATAAATAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtga                         <  2:501155/89‑1 (MQ=255)
                                                           ttCAGAATTATTCATCACTCCTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  2:303573/90‑1 (MQ=255)
                                                           ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTACTTGCGTATGTTCAGAGtgat                        <  1:559230/90‑1 (MQ=37)
                                                           ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  2:381251/90‑1 (MQ=255)
                                                           ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  2:417317/90‑1 (MQ=255)
                                                           ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  2:364631/90‑1 (MQ=255)
                                                           ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  2:825452/90‑1 (MQ=255)
                                                           ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  1:837294/90‑1 (MQ=255)
                                                           ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  1:768984/90‑1 (MQ=255)
                                                           ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        <  1:591017/90‑1 (MQ=255)
                                                           ttCAGAATTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                        >  1:254374/1‑90 (MQ=255)
                                                                aaTTATTCATCACTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGatat                   <  1:566655/90‑1 (MQ=255)
                                                                        atcaCTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGa           <  1:68838/90‑1 (MQ=255)
                                                                         tcaCTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGaa          >  2:319645/1‑90 (MQ=255)
                                                                         tcaCTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGaa          >  2:512162/1‑90 (MQ=255)
                                                                         tcaCTACTCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGaa          >  1:429683/1‑90 (MQ=255)
                                                                             tactCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGAtt      <  1:751411/90‑1 (MQ=255)
                                                                             tactCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGAtt      <  1:521686/90‑1 (MQ=255)
                                                                              actCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGAttt     >  1:303506/1‑90 (MQ=255)
                                                                              actCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGAttt     >  1:25247/1‑90 (MQ=255)
                                                                              actCCTCTGAAAATAATAAAAAACGGTAGAATCGTGAGATCCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGAttt     >  2:74040/1‑90 (MQ=37)
                                                                                 cctcTGAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGAtttttt  >  2:37080/1‑90 (MQ=255)
                                                                                       |                                                                                    
AGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGATTTTTT  >  NZ_CP009273/264598‑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: