Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,629,389 (T)8→7 intergenic (+42/+18) creD → / ← arcA cell envelope integrity protein CreD/two‑component system response regulator ArcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,629,3820T.96.4% 107.3 / ‑2.3 28intergenic (+35/+25)creD/arcAcell envelope integrity protein CreD/two‑component system response regulator ArcA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (15/12);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTC  >  NZ_CP009273/4629295‑4629466
                                                                                       |                                                                                    
tGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGttt                                                                                    >  1:939070/1‑90 (MQ=255)
  aaaGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGttttt                                                                                  >  1:816309/1‑90 (MQ=255)
       cAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCTTTTTTTTGAc                                                                            <  1:883892/90‑1 (MQ=255)
        aGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGtttttttgacg                                                                            <  1:320614/90‑5 (MQ=255)
        aGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGtttttttgacg                                                                            <  1:909751/90‑5 (MQ=255)
        aGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGtttttttgacg                                                                            <  1:55177/90‑5 (MQ=255)
        aGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGtttttttgacg                                                                            <  1:564638/90‑5 (MQ=255)
        aGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGtttttttgacg                                                                            <  2:116233/90‑5 (MQ=255)
               aagTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTaa                                                                    <  2:145109/90‑1 (MQ=255)
                agTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTaaa                                                                   <  1:265372/90‑1 (MQ=255)
                agTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTaaa                                                                   >  2:591034/1‑90 (MQ=255)
                 gTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAg                                                                  <  1:86661/90‑1 (MQ=255)
                          gACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAAt                                                         <  2:20714/90‑1 (MQ=255)
                          gACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAAt                                                         <  2:469725/90‑1 (MQ=255)
                            cGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAATCt                                                       <  2:751059/90‑1 (MQ=255)
                              aTGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTc                                                     >  1:720074/1‑90 (MQ=255)
                                    tACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATc                                               >  2:540183/1‑90 (MQ=255)
                                     aCGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACAGTGGTAAAGCCGATTAATCTTCCAGATCa                                              >  2:102309/1‑90 (MQ=255)
                                         aTCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGc                                          <  2:552953/90‑1 (MQ=255)
                                            tGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGa                                       >  2:102448/1‑90 (MQ=255)
                                            tGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGa                                       <  2:527827/90‑1 (MQ=255)
                                            tGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGa                                       >  1:846808/1‑90 (MQ=255)
                                            tGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGa                                       >  1:658795/1‑90 (MQ=255)
                                                  tAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGAt                                 >  1:247837/1‑90 (MQ=255)
                                                    aGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATaa                               >  1:903117/1‑90 (MQ=255)
                                                        tGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACgg                                                                          >  2:491/1‑43 (MQ=38)
                                                        tGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACgg                                                                          <  1:491/43‑1 (MQ=38)
                                                        tGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGGTAACCtt                           >  1:952457/1‑90 (MQ=255)
                                                        tGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCtt                           >  1:785893/1‑90 (MQ=255)
                                                        tGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCtt                           >  2:656521/1‑90 (MQ=255)
                                                                 aaaaCGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGaa                  <  1:495749/90‑1 (MQ=255)
                                                                  aaaCGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAAt                 <  1:70150/90‑1 (MQ=255)
                                                                            aaaaaGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCgatg       <  2:799088/90‑1 (MQ=255)
                                                                                 gcgcCGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtc  >  1:151291/1‑90 (MQ=255)
                                                                                 gcgcCGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtc  >  1:308190/1‑90 (MQ=255)
                                                                                       |                                                                                    
TGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTC  >  NZ_CP009273/4629295‑4629466

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: