Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,057,557 A→G 9.3% I172T (ATA→ACA)  PP_0914 ← GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,057,5570AG 9.3% 113.7 / 3.1 43I172T (ATA→ACA) PP_0914GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (19/20);  new base G (2/2);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

GAGCAGCGGTATCTGGTGCTGGCTGGGGCGGTCAGGGTTGATCTGGGTGAGGTGCTCGATGATCTTCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGCCGTAGTATTCGCTGTAAGGTCCACGGCCGACCAGTTGCGTGGTGTAGGTGCGCTCCTGGCCGAGGATGAG  >  NC_002947/1057412‑1057699
                                                                                                                                                  |                                                                                                                                              
gAGCAGCGGTATCTGGTGCTGGCTGGGGCGGTCAGGGTTGATCTGGGTGAGGTGCTCGATGATCTTCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATg                                                                                                                                              >  6:340026/1‑148 (MQ=255)
    agcGGTATCTGGTGCTGGCTGGGGCGGTCAGGGTTGATCTGGGTGAGGTGCTCGATGATCTTCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGATCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAgat                                                                                                                                         <  1:276975/149‑1 (MQ=255)
     gcGGTATCTGGTGCTGGCTGGGGCGGTCAGGGTTGATCTGGGTGAGGTGCTCGATGATCTTCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAgatg                                                                                                                                        <  5:340026/149‑1 (MQ=255)
                         gggCGGTCAGGGTTGATCTGGGTGAGGTGCTCGATGATCTTCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCAACGCGGTGATCatctt                                                                                                                    <  3:93252/149‑5 (MQ=255)
                                  gggTTGATCTGGGTGAGGTGCTCGATGATCTTCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGcc                                                                                                           >  5:225791/1‑149 (MQ=255)
                                        aTCTGGGTGAGGTGCTCGATGATCTTCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATGTGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTag                                                                                                       >  4:104717/1‑147 (MQ=255)
                                         tcagggTGAGGTGCTCGATGATCTTCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGTCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGATGATCAGGAAATCTTTGCCGAaggag                                                                                                    >  7:230532/4‑149 (MQ=255)
                                         tCTGGGTGAGGTGCTCGATGATCTTCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGTCGTaggag                                                                                                    >  2:303998/1‑149 (MQ=255)
                                            gggTGAGGTGCTCGATGATCTTCGACAGCGGTTTGGTCAGCAGGCCAGTGGTAGCCCAGTTAAACCACCAGGCCCATGGCCATGGCCCGCAGCACGCAGGATATGAAGATGAACCCCCCGTTGAACAGGAAATCTTTGCCGTaggag                                                                                                    >  1:69110/1‑147 (MQ=255)
                                            gggTGAGGTGCTCGATGATCTTCGACAGCGGTTTGGTCAGCA‑GCGAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCttt                                                                                                              <  4:248790/136‑1 (MQ=255)
                                            gggTGAGGTGCTCGATGATCTTCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCttt                                                                                                              >  3:248790/1‑136 (MQ=255)
                                                     gCTCGATGATCTTCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAg                                                                                                                       <  4:132794/118‑1 (MQ=255)
                                                         gatgatCTTCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTaggaggag                                                                                                 >  1:275688/1‑136 (MQ=255)
                                                         gatgatCTTCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTaggaggag                                                                                                 <  2:275688/136‑1 (MQ=255)
                                                          atgatCTTCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTaggaggag                                                                                                 <  5:287140/135‑1 (MQ=255)
                                                          atgatCTTCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTaggaggag                                                                                                 >  6:287140/1‑135 (MQ=255)
                                                           tgatCTTCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATGTGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATg                                                                                  <  3:104717/149‑1 (MQ=255)
                                                                ttCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTg                                                                                                             >  6:213333/1‑117 (MQ=255)
                                                                ttCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTg                                                                                                             <  5:213333/117‑1 (MQ=255)
                                                                ttCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATCAAGATGATAACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGTGGTATCGCGGGTAATGCTGAg                                                                             >  8:101205/1‑149 (MQ=255)
                                                                  cGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGTCGTAGGAGGAGGTATCGAGGGTAATGCTGAgg                                                                            <  1:303998/148‑1 (MQ=255)
                                                                  cGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGa                                                                              <  8:2646/146‑1 (MQ=255)
                                                                   gACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGt                                                                           <  6:225791/148‑1 (MQ=255)
                                                                                tCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGCCGTAGTATTCg                                                             >  5:315019/1‑149 (MQ=255)
                                                                                 cagca‑gCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGCCGTAGTATTCGc                                                            >  1:269315/1‑149 (MQ=255)
                                                                                        ccAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGCCGTAGTATTCGCTGTAAg                                                      <  8:343534/149‑1 (MQ=255)
                                                                                        ccAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGCCGTAGTATTCGCTGTAAg                                                      >  7:343534/1‑149 (MQ=255)
                                                                                               tAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGAGATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGCCGTAGTATTCGCTGTAAGGTCCACg                                               >  4:335145/1‑149 (MQ=255)
                                                                                                 gACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGcc                                                                       >  4:110485/1‑123 (MQ=255)
                                                                                                 gACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGcc                                                                       <  3:110485/123‑1 (MQ=255)
                                                                                                     aGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGCCGTAGTATTCGCTGTAAGGTCCACGGCCGAc                                         >  4:81651/1‑149 (MQ=255)
                                                                                                            gCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGTAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGTGGTCGCCGTAGTATTCGCTGTAAGGTCCACGGCCGACCAGTTGc                                  <  2:69110/149‑1 (MQ=255)
                                                                                                               ccAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGCCGTAGTATTCGCTGTAAGGTCCACGGCCGACCAGTTGCgtg                               >  7:194266/1‑149 (MQ=255)
                                                                                                                     ccatggccatggCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGCCGTAGTATTCGCTGTAAGGTCCACGGCCGACCAGTTGCGTGGTGTAg                         <  6:315019/149‑1 (MQ=255)
                                                                                                                      catggccatggcCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGCCGTAGTATTCGCTGTAAGGTCCACGGCCGACCAGTTGCGTGGTGTAgg                        <  2:6373/149‑1 (MQ=255)
                                                                                                                                gccCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGCCGTAGTATTCGCTGTAAGGTCCACGGCCGACCAGTTGCGTGgtgt                           >  3:148168/1‑136 (MQ=255)
                                                                                                                                gccCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGCCGTAGTATTCGCTGTAAGGTCCACGGCCGACCAGTTGCGTGgtgt                           <  4:148168/136‑1 (MQ=255)
                                                                                                                                 cccGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGCCGTAGTATTCGCTGTAAGGTCCACGGCCGACCAGTTGCGTGGTGTAGGTGCGCTCCTgg             <  8:194266/149‑1 (MQ=255)
                                                                                                                                       gcaCGCCGGATGTGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGCCGTAGTATTCGCTGTAAGGTCCACGGCCGACCAGtt                                    <  6:301024/120‑1 (MQ=255)
                                                                                                                                       gcaCGCCGGATGTGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGCCGTAGTATTCGCTGTAAGGTCCACGGCCGACCAGtt                                    >  5:301024/1‑120 (MQ=255)
                                                                                                                                           gCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGCCGTAGTATTCGCTGTAAGGTCCACGGCCGACCAGTTGCGTGGTGTAGGTGCGCTCCt               >  7:70412/1‑137 (MQ=255)
                                                                                                                                           gCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTAGCGAGGGTAATGCTGAGGTCGCCGTAGTATTCGCTGTAAGGTCCACGGCCGACCAGTTGCGTGGTGTAGGTGCGCTCCt               <  8:70412/137‑1 (MQ=255)
                                                                                                                                             cGGATATGAAGATGATCACCGCGTTGATCAGTAAATCTTTTCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGCCGTAGTATTCGCTGTAAGGTCCACGGCCGACCAGTTGCGTGGTTTAGGTGCGCTCCTGGCCGAGGATGAg  <  7:101205/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GAGCAGCGGTATCTGGTGCTGGCTGGGGCGGTCAGGGTTGATCTGGGTGAGGTGCTCGATGATCTTCGACAGCGGTTTGGTCAGCA‑GCCAGTGGTAGACCAGGTAAAGCACCAGGCCCATGGCCATGGCCCGCAGCACGCCGGATATGAAGATGATCACCGCGTTGATCAGGAAATCTTTGCCGTAGGAGGAGGTATCGAGGGTAATGCTGAGGTCGCCGTAGTATTCGCTGTAAGGTCCACGGCCGACCAGTTGCGTGGTGTAGGTGCGCTCCTGGCCGAGGATGAG  >  NC_002947/1057412‑1057699

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑