Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,649,476 A→G intergenic (‑78/+128) rspA ← / ← ynfA starvation‑sensing protein RspA/YnfA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,649,4760AG95.8% 144.0 / ‑0.9 48intergenic (‑78/+128)rspA/ynfAstarvation‑sensing protein RspA/YnfA family protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (26/20);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

CGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAG  >  NZ_CP009273/1649389‑1649560
                                                                                       |                                                                                    
cGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGgaa                                                                                    <  1:69929/90‑1 (MQ=255)
   tCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGgaaaaa                                                                                 >  1:568499/1‑90 (MQ=255)
   tCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGgaaaaa                                                                                 >  2:1398018/1‑90 (MQ=255)
      tCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGcc                                                                              >  2:1279560/1‑90 (MQ=255)
      tCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGcc                                                                              >  1:18121/1‑90 (MQ=255)
          gTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAg                                                                          >  2:693728/1‑90 (MQ=255)
          gTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAg                                                                          <  1:619176/90‑1 (MQ=255)
            tCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTg                                                                        >  2:1397154/1‑90 (MQ=255)
            tCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTg                                                                        >  1:885530/1‑90 (MQ=255)
            tCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTg                                                                        >  1:802006/1‑90 (MQ=255)
            tCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTg                                                                        >  2:1199033/1‑90 (MQ=255)
              gATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGAt                                                                      >  2:310120/1‑90 (MQ=255)
                ttCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATca                                                                    <  1:1317305/90‑1 (MQ=255)
                       cATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACaaaaaa                                                             <  2:802006/90‑1 (MQ=255)
                           gCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGc                                                         >  1:1060000/1‑90 (MQ=255)
                             ttGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGt                                                       <  2:1437131/90‑1 (MQ=255)
                             ttGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGt                                                       <  1:844104/90‑1 (MQ=255)
                             ttGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGt                                                       <  1:1154806/90‑1 (MQ=255)
                                      tGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTgcg                                              <  2:1022075/90‑1 (MQ=255)
                                       gCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTgcgc                                             <  1:1379623/90‑1 (MQ=255)
                                       gCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTgcgc                                             <  2:298163/90‑1 (MQ=255)
                                       gCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTgcgc                                             <  1:938832/90‑1 (MQ=255)
                                       gCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTgcgc                                             <  1:693728/90‑1 (MQ=255)
                                            aaaTCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGa                                        >  2:1232699/1‑90 (MQ=255)
                                              aTCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAAt                                      <  1:1398018/90‑1 (MQ=255)
                                              aTCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAAt                                      <  1:1185578/90‑1 (MQ=255)
                                                cTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATgg                                    <  2:330335/90‑1 (MQ=255)
                                                 tACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGt                                   >  1:254688/1‑90 (MQ=255)
                                                  aCGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGtt                                  <  2:1075912/90‑1 (MQ=255)
                                                  aCGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGtt                                  >  1:138393/1‑90 (MQ=255)
                                                          gAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAg                          >  2:616423/1‑90 (MQ=255)
                                                            gCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGAt                        <  2:1326146/90‑1 (MQ=255)
                                                              tactacCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGa                      <  1:1199033/90‑1 (MQ=255)
                                                                ctacCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAAt                   <  2:652386/90‑1 (MQ=255)
                                                                    cATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGa                >  1:118002/1‑90 (MQ=255)
                                                                    cATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGa                >  1:46585/1‑90 (MQ=255)
                                                                    cATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGa                >  1:805123/1‑90 (MQ=255)
                                                                      tACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAAt              <  1:1279560/90‑1 (MQ=255)
                                                                          aGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTgcgc                                             >  1:272071/1‑55 (MQ=255)
                                                                          aGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTgcgc                                             <  2:272071/55‑1 (MQ=255)
                                                                          aGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCt          >  2:638282/1‑90 (MQ=255)
                                                                          aGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCt          >  1:1510798/1‑90 (MQ=255)
                                                                           gTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCtt         >  1:1313221/1‑90 (MQ=255)
                                                                            tataAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCttt        >  1:454943/1‑90 (MQ=255)
                                                                             ataAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCtttt       >  1:378423/1‑90 (MQ=255)
                                                                             ataAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCtttt       >  2:130565/1‑90 (MQ=255)
                                                                             ataAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCtttt       >  2:1121897/1‑90 (MQ=255)
                                                                                  gATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAg  >  1:615074/1‑90 (MQ=255)
                                                                                  gATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAg  >  2:1509843/1‑90 (MQ=255)
                                                                                  gATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAg  >  2:751642/1‑90 (MQ=255)
                                                                                    tCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTgcg                                              <  2:1057907/44‑1 (MQ=255)
                                                                                       |                                                                                    
CGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAG  >  NZ_CP009273/1649389‑1649560

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: