Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,903,061 A→T intergenic (‑38/‑121) mgrB ← / → yobH PhoP/PhoQ regulator MgrB/YobH family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,903,0610AT98.0% 188.6 / ‑3.4 50intergenic (‑38/‑121)mgrB/yobHPhoP/PhoQ regulator MgrB/YobH family protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (25/24);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01

AGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTAGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCTA  >  NZ_CP009273/1902979‑1903147
                                                                                  |                                                                                      
cgcCAACACCACCCCAACCCGAACCCCCCATCGAAACTTTTTCACTCCACTCTCCCTTTCACCTCCCTTATGTCATTAAACTTGCATGCg                                                                                 <  2:1154111/89‑1 (MQ=255)
  ccAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCgtt                                                                               <  2:357568/90‑1 (MQ=255)
  ccAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCgtt                                                                               <  2:45265/90‑1 (MQ=255)
   cAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCgtta                                                                              >  2:759150/1‑90 (MQ=255)
     acacCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCgttagt                                                                            >  1:277220/1‑90 (MQ=255)
          acGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc                                                                       <  1:1015823/90‑1 (MQ=255)
          acGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc                                                                       <  1:981514/90‑1 (MQ=255)
          acGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc                                                                       >  2:1439833/1‑90 (MQ=255)
          acGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc                                                                       <  1:320509/90‑1 (MQ=255)
          acGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc                                                                       <  2:594272/90‑1 (MQ=255)
          acGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc                                                                       <  2:655233/90‑1 (MQ=255)
          acGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc                                                                       <  2:852763/90‑1 (MQ=255)
          acGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc                                                                       <  1:1086532/90‑1 (MQ=255)
           cGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCgtt                                                                               >  2:873146/1‑81 (MQ=255)
           cGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCgtt                                                                               <  1:873146/81‑1 (MQ=255)
           cGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACa                                                                      >  2:1205691/1‑90 (MQ=255)
           cGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACa                                                                      <  2:933959/90‑1 (MQ=255)
           cGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACa                                                                      <  1:532351/90‑1 (MQ=255)
               aaCCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAgcgc                                                                  >  1:782011/1‑90 (MQ=255)
                aCCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGcc                                                                 >  1:1175302/1‑90 (MQ=255)
                 ccAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCt                                                                <  2:349651/90‑1 (MQ=255)
                    gAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAAc                                                             >  2:1335184/1‑90 (MQ=255)
                         aCCCATCGAGACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGt                                                        >  2:1173775/1‑90 (MQ=255)
                         aCCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGt                                                        >  2:413254/1‑90 (MQ=255)
                          cccATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTc                                                       >  1:210081/1‑90 (MQ=255)
                           ccATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCg                                                      >  1:1028131/1‑90 (MQ=255)
                                    cTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc                                                                       <  2:1555220/64‑1 (MQ=255)
                                    cTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc                                                                       >  1:1555220/1‑64 (MQ=255)
                                       tttCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTAGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCtt                                          >  1:541125/1‑90 (MQ=255)
                                               cAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGaa                                  <  1:628314/90‑1 (MQ=255)
                                                      cGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCg                           >  2:916649/1‑90 (MQ=255)
                                                      cGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCg                           >  1:922106/1‑90 (MQ=255)
                                                          tCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCgcttg                       >  2:766154/1‑90 (MQ=255)
                                                              ctacctTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTc                   <  2:1051720/90‑1 (MQ=255)
                                                               tacctTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCa                  >  2:496663/1‑90 (MQ=255)
                                                                acctTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAg                 <  2:1475153/90‑1 (MQ=255)
                                                                acctTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAg                 <  1:693133/90‑1 (MQ=255)
                                                                    tATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCa             >  2:1262256/1‑90 (MQ=255)
                                                                     aTGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCat            >  1:1099403/1‑90 (MQ=255)
                                                                      tGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCata           >  2:508471/1‑90 (MQ=255)
                                                                      tGTCATCAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCgcttgctt                    >  2:1298843/1‑81 (MQ=255)
                                                                      tGTCATCAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCgcttgctt                    <  1:1298843/81‑1 (MQ=255)
                                                                         cATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTg        <  1:871953/90‑1 (MQ=255)
                                                                          aTTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGc       >  2:479757/1‑90 (MQ=255)
                                                                          aTTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGc       >  1:1508520/1‑90 (MQ=255)
                                                                           ttAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGcc      >  1:370239/1‑90 (MQ=255)
                                                                           ttAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGcc      >  1:816347/1‑90 (MQ=255)
                                                                             aaaCTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGc    <  2:210081/90‑1 (MQ=255)
                                                                               aCTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCTa  <  1:959429/90‑1 (MQ=255)
                                                                               aCTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCTa  <  1:285571/90‑1 (MQ=255)
                                                                                  |                                                                                      
AGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTAGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCTA  >  NZ_CP009273/1902979‑1903147

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