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,4760AG98.4% 231.0 / ‑3.2 62intergenic (‑78/+128)rspA/ynfAstarvation‑sensing protein RspA/YnfA family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (30/31);  minor base . (1/0);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.85e-01

TCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGC  >  NZ_CP009273/1649392‑1649561
                                                                                    |                                                                                     
tCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGgaaaaa                                                                                  >  1:677253/1‑90 (MQ=255)
   tCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGcc                                                                               >  2:2767877/1‑90 (MQ=255)
    cATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCg                                                                              <  1:2599910/90‑1 (MQ=255)
         tCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTg                                                                         >  2:556509/1‑90 (MQ=255)
           gATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGAt                                                                       >  1:1050556/1‑90 (MQ=255)
             ttCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATca                                                                     <  1:1492298/90‑1 (MQ=255)
                cTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACaa                                                                  >  2:774039/1‑90 (MQ=255)
                   gCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACaaaaa                                                               <  1:122563/90‑1 (MQ=255)
                    cATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACaaaaaa                                                              <  2:1471204/90‑1 (MQ=255)
                    cATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACaaaaaa                                                              <  1:1667922/90‑1 (MQ=255)
                    cATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACaaaaaa                                                              >  2:1340543/1‑90 (MQ=255)
                     aTCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACaaaaaaa                                                             <  2:1050556/90‑1 (MQ=255)
                        gCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGc                                                          >  1:2968542/1‑90 (MQ=255)
                        gCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGc                                                          >  1:1520151/1‑90 (MQ=255)
                        gCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGc                                                          >  2:122825/1‑90 (MQ=255)
                         cTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCg                                                         <  1:926711/90‑1 (MQ=255)
                          ttGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGt                                                        <  2:123580/90‑1 (MQ=255)
                                    gCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTgcgc                                              <  2:1514516/90‑1 (MQ=255)
                                    gCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTgcgc                                              <  2:1980472/90‑1 (MQ=255)
                                    gCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTgcgc                                              <  1:2414533/90‑1 (MQ=255)
                                      aTGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCtg                                            >  2:2071929/1‑90 (MQ=255)
                                      aTGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCtg                                            >  2:2935272/1‑90 (MQ=255)
                                           aTCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAAt                                       <  2:1377179/90‑1 (MQ=255)
                                           aTCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAAt                                       <  1:597330/90‑1 (MQ=255)
                                           aTCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAAt                                       <  1:1641137/90‑1 (MQ=255)
                                              tACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGt                                    >  1:2428234/1‑90 (MQ=255)
                                              tACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGt                                    >  1:884897/1‑90 (MQ=255)
                                                    aCTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACaa                              <  2:1800884/90‑1 (MQ=255)
                                                    aCTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACaa                              >  2:476779/1‑90 (MQ=255)
                                                    aCTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACaa                              <  2:2968542/90‑1 (MQ=255)
                                                    aCTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACaa                              <  2:845246/90‑1 (MQ=255)
                                                      tGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACaaaa                            <  1:2071929/90‑1 (MQ=255)
                                                       gAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAg                           >  2:1100843/1‑90 (MQ=255)
                                                       gAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAg                           >  1:2278529/1‑90 (MQ=255)
                                                        aGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGa                          <  1:293813/90‑1 (MQ=255)
                                                         gCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGAt                         <  2:1551008/90‑1 (MQ=255)
                                                         gCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGAt                         <  2:833462/90‑1 (MQ=255)
                                                           tactacCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGa                       <  2:1425838/90‑1 (MQ=255)
                                                              tacCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAAt                    <  1:3023450/90‑1 (MQ=255)
                                                                 cATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGa                 >  1:2077606/1‑90 (MQ=255)
                                                                 cATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGa                 >  2:2488946/1‑90 (MQ=255)
                                                                 cATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGa                 >  2:1430415/1‑90 (MQ=255)
                                                                 cATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGa                 >  1:224118/1‑90 (MQ=255)
                                                                  aTACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGaa                >  2:1607443/1‑90 (MQ=255)
                                                                   tACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAAt               <  2:129918/90‑1 (MQ=255)
                                                                   tACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAAt               >  1:395193/1‑90 (MQ=255)
                                                                      aaGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAAt                                       <  1:800806/63‑1 (MQ=255)
                                                                      aaGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAAt                                       >  2:800806/1‑63 (MQ=255)
                                                                       aGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGa                       <  2:294169/78‑1 (MQ=255)
                                                                       aGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGa                       >  1:294169/1‑78 (MQ=255)
                                                                       aGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCt           >  1:647504/1‑90 (MQ=255)
                                                                       aGTATAAAGATCGGAAAAGCCGGAGTGATCACAAAAAAAGACGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCtt          >  2:944314/1‑90 (MQ=255)
                                                                         tataAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGa                       <  2:1914053/76‑1 (MQ=255)
                                                                         tataAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGa                       >  1:1914053/1‑76 (MQ=255)
                                                                         tataAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCttt         >  2:1488269/1‑90 (MQ=255)
                                                                          ataAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCtttt        <  1:476779/90‑1 (MQ=255)
                                                                          ataAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCtttt        >  1:2649426/1‑90 (MQ=255)
                                                                            aaaGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGt      <  2:2564545/90‑1 (MQ=255)
                                                                             aaGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTc     <  1:862949/90‑1 (MQ=255)
                                                                             aaGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTc     >  1:1461836/1‑90 (MQ=255)
                                                                               gATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAg   <  1:2935272/90‑1 (MQ=255)
                                                                               gATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAg   >  2:857028/1‑90 (MQ=255)
                                                                                aTCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGc  >  2:1060618/1‑90 (MQ=255)
                                                                                    |                                                                                     
TCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGC  >  NZ_CP009273/1649392‑1649561

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