Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 787,882 +T intergenic (‑662/+630) h16_A0728 ← / ← h16_A0729 response regulator containing CheY‑like receiver, AAA‑type ATPase, and DNA‑binding domains/transcriptional regulator, LysR‑family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479787,8791.T96.8% 114.6 / ‑3.2 31intergenic (‑659/+633)h16_A0728/h16_A0729response regulator containing CheY‑like receiver, AAA‑type ATPase, and DNA‑binding domains/transcriptional regulator, LysR‑family
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (16/14);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

TCGGCGAACCCGTGCCAGTGTGTGCCGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATGGATCACAGTGCGTGG  >  AM260479/787747‑788007
                                                                                                                                      |                                                                                                                                
tCGGCGAACCCGTGCCAGTGTGTTGCCGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAaggaaaaggaaa                                                                                                                                   >  2:359730/1‑134 (MQ=255)
tCGGCGAACCCGTGCCAGTGTGTTGCCGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAaggaaaaggaaa                                                                                                                                   >  2:325765/1‑134 (MQ=255)
            tGCCAGTGTGTTGCCGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGc                                                                                                                       >  1:118936/1‑134 (MQ=255)
                        gCCGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAAtt                                                                                                            <  1:272220/133‑1 (MQ=255)
                          cGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCt                                                                                                                          <  1:235391/117‑1 (MQ=255)
                          cGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCt                                                                                                                          >  2:235391/1‑117 (MQ=255)
                                                cTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCtt                                                                                   >  2:804384/1‑134 (MQ=255)
                                                         gCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACAt                                                                          >  2:290694/1‑134 (MQ=255)
                                                            tCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATc                                                                         <  2:836342/132‑1 (MQ=255)
                                                            tCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATc                                                                         >  1:836342/1‑132 (MQ=255)
                                                            tCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTg                                                                       >  2:515433/1‑134 (MQ=255)
                                                                ttCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCg                                                                    <  2:952113/133‑1 (MQ=255)
                                                                ttCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCg                                                                    <  1:804384/133‑1 (MQ=255)
                                                                ttCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCg                                                                    <  1:359730/133‑1 (MQ=255)
                                                                        gACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCtat                                                          <  1:887646/134‑1 (MQ=255)
                                                                            tatGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATAtt                                                       >  1:510226/1‑134 (MQ=255)
                                                                                         cGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCAtt                                          >  1:958097/1‑134 (MQ=255)
                                                                                           tGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCg                                        >  2:630184/1‑134 (MQ=255)
                                                                                              cAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCg                                     <  1:325765/134‑1 (MQ=255)
                                                                                                          aTCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCt                                                                                                         <  1:783544/54‑1 (MQ=255)
                                                                                                          aTCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCt                                                                                                         >  2:783544/1‑54 (MQ=255)
                                                                                                          aTCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAAtt                          <  1:515433/133‑1 (MQ=255)
                                                                                                          aTCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAAtt                          <  2:510226/133‑1 (MQ=255)
                                                                                                           tCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAg                        <  1:290694/134‑1 (MQ=255)
                                                                                                            ccTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCg                                                                    <  2:578930/89‑1 (MQ=255)
                                                                                                            ccTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCg                                                                    >  1:578930/1‑89 (MQ=255)
                                                                                                              tttGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGgcg                     >  2:269563/1‑134 (MQ=255)
                                                                                                                tGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCa                   >  1:748753/1‑134 (MQ=255)
                                                                                                                    aaCAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTg                                                                       <  2:614077/78‑1 (MQ=255)
                                                                                                                    aaCAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTg                                                                       >  1:614077/1‑78 (MQ=255)
                                                                                                                    aaCAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAg                        <  2:335285/125‑1 (MQ=255)
                                                                                                                    aaCAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAg                        >  1:335285/1‑125 (MQ=255)
                                                                                                                    aaCAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATGGa               >  2:179145/1‑134 (MQ=255)
                                                                                                                          aggaaaaggaaatTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATgg                <  2:492812/127‑1 (MQ=255)
                                                                                                                          aggaaaaggaaatTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATgg                >  1:492812/1‑127 (MQ=255)
                                                                                                                          aggaaaaggaaatTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATGGATCACAg         >  2:404748/1‑134 (MQ=255)
                                                                                                                           ggaaaaggaaatTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATc                                                                         <  2:661129/69‑1 (MQ=255)
                                                                                                                           ggaaaaggaaatTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATc                                                                         <  2:231895/69‑1 (MQ=255)
                                                                                                                           ggaaaaggaaatTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATc                                                                         >  1:661129/1‑69 (MQ=255)
                                                                                                                           ggaaaaggaaatTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATc                                                                         >  1:231895/1‑69 (MQ=255)
                                                                                                                                 ggaaatTTTGGGCTGGGCTTGGGTAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATGGATCACAGTGCGTgg  <  2:748753/134‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                
TCGGCGAACCCGTGCCAGTGTGTGCCGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATGGATCACAGTGCGTGG  >  AM260479/787747‑788007

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