Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,258,258 C→T 7.8% intergenic (+313/‑79) PP_1990 → / → asd transposase/aspartate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,258,2580CT 7.8% 143.9 / 2.2 51intergenic (+313/‑79)PP_1990/asdtransposase/aspartate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (22/25);  new base T (2/2);  total (24/27)
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.26e-01

CGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTGGATATCGCCGTCGTCGGCGCCACCGGCAGCGTCGGTGAAGCCTTGGTGC  >  NC_002947/2258110‑2258400
                                                                                                                                                    |                                                                                                                                              
cgcTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAAc                                                                                                                                                >  3:238739/1‑149 (MQ=255)
 gcTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAAc                                                                                                                                                >  5:174138/1‑148 (MQ=255)
               cATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAACTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGcaca                                                                                                                                 <  2:152783/149‑1 (MQ=255)
                aTGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGCGGcactt                                                                                                                                >  2:180929/1‑147 (MQ=255)
                   cccGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGAc                                                                                                                              <  6:255772/148‑1 (MQ=255)
                    ccGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACaa                                                                                                                            <  3:48908/149‑1 (MQ=255)
                          ccTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCACCAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAAccc                                                                                                                       >  5:144305/1‑148 (MQ=255)
                           cTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTa                                                                                                                     <  3:274642/149‑1 (MQ=255)
                           cTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTa                                                                                                                     <  5:237182/149‑1 (MQ=255)
                           cTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGAACATCCAAAACCCTa                                                                                                                     >  6:237182/1‑149 (MQ=255)
                                  cGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGATAGAGCAATAGGCAACGCAGGGCGAGGCACAGACAACACACTACAATcag                                                                                                              >  8:69446/1‑147 (MQ=21)
                                    cAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAACATCAAGATAAACAGCAACAGCAACAGCAATAACAACACCAATAGGCAACTCTGGGGGAGGCACAGACAAAACCCTa                                                                                                                     >  7:107303/1‑140 (MQ=21)
                                           tGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGa                                                                                                             >  3:176496/1‑141 (MQ=255)
                                           tGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGa                                                                                                             <  4:176496/141‑1 (MQ=255)
                                                    cctcctGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACATCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCccg                                                                                            <  6:174138/149‑1 (MQ=255)
                                                        ctGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGt                                                                                        >  3:144817/1‑149 (MQ=255)
                                                            aGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTAc                                                                                    <  4:238739/149‑1 (MQ=255)
                                                              tCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGCTCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATACGAGTAAATTGCCGCCCCCGCCGTTTACAg                                                                                  >  7:85659/1‑149 (MQ=255)
                                                                gcgAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGc                                                                                 <  1:136018/148‑1 (MQ=255)
                                                                     gatcaagatcaagatcaagatcaagatcaagatcaagatcaaCAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAgg                                                                           <  5:285965/149‑1 (MQ=255)
                                                                        caagatcaagatcaagatcaagatcaagatcaagatcaaCAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATc                                                                        <  3:174190/149‑1 (MQ=255)
                                                                           gatcaagatcaagatcaagatcaacatcaagatcaaCAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGATGGGGGAGGCACAGACAAAACCCTACATACCCAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCt                                                                     >  4:131719/1‑149 (MQ=38)
                                                                               aagatcaagatcaagatcaagatcaagatcaaCAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCa                                                                 >  8:247523/1‑149 (MQ=255)
                                                                                     atgatcaagatcaagatcaagatcaaCAGCAACAGCAACAGCAATAGCTACATCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTTCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGacac                                                           <  7:69446/147‑1 (MQ=255)
                                                                                        atcaagatcaagatcaagatcaaCAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAATGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATc                                                                        <  4:1603/133‑1 (MQ=255)
                                                                                        atcaagatcaagatcaagatcaaCAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAATGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATc                                                                        <  8:254315/133‑1 (MQ=255)
                                                                                        atcaagatcaagatcaagatcaaCAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAATGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATc                                                                        >  7:254315/1‑133 (MQ=255)
                                                                                        atcaagatcaagatcaagatcaaCAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAATGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATc                                                                        >  3:1603/1‑133 (MQ=255)
                                                                                            agatcaagatcaagatcaaCAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCttt                                                    <  4:144817/149‑1 (MQ=255)
                                                                                            agatcaagatcaagatcaaCAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCttt                                                    <  7:247523/149‑1 (MQ=255)
                                                                                             gatcaagatcaagatcaaCAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACa                                                                                   <  3:241331/117‑1 (MQ=255)
                                                                                                     tcaagatcaaCAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTGGATATCGc                                           >  2:286662/1‑149 (MQ=255)
                                                                                                        agatcaaCAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTGGATATCGCcgt                                        >  2:231707/1‑149 (MQ=255)
                                                                                                               cagcaacagcaacagcaaTAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTGGATATCGCCGTCGTCGgc                                 <  1:286662/149‑1 (MQ=255)
                                                                                                                agcaacagcaacagcaaTAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTGGATATCGCCGTCGTCGgcg                                >  5:134576/1‑149 (MQ=255)
                                                                                                                  caacagcaacagcaaTAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTGGATATCGCCGTCGTCGgcgc                               >  6:18154/1‑148 (MQ=255)
                                                                                                                  caacagcaacagcaaTAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTGGATATCGCCGTCGTCGgcgc                               >  7:40587/1‑148 (MQ=255)
                                                                                                                     cagcaacagcaaTAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGaca                                                            >  7:104090/1‑116 (MQ=255)
                                                                                                                     cagcaacagcaaTAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGaca                                                            <  8:104090/116‑1 (MQ=255)
                                                                                                                       gcaacagcaaTAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTGGATATCGCCGTCGTCGGCGCCACCgg                         >  5:164148/1‑149 (MQ=255)
                                                                                                                        caacagcaaTAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACCCCCTTTGGATATCGCCGTCGTCGGCGCCACCgg                         >  3:188687/1‑148 (MQ=255)
                                                                                                                         aacagcaaTAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTGGATATCGCCGTCGTCGGCGCCACCGGCa                       <  6:144305/149‑1 (MQ=255)
                                                                                                                                 tAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTGGATATCGCCGTCGTCGGCGCCACCGGCAGCGTCGGt               <  5:18154/149‑1 (MQ=255)
                                                                                                                                     tACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTGGATATCGCCGTCGTCGGCGCCACCGGCAGCGTCGGTGAAg           >  2:2372/1‑149 (MQ=255)
                                                                                                                                      aCAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTgg                                                  <  8:236679/109‑1 (MQ=255)
                                                                                                                                      aCAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTgg                                                  >  7:236679/1‑109 (MQ=255)
                                                                                                                                         gCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTGGATATCGCCGTCGTCGGCGCCACCGGCAGCGTCGGTGAAGCCtt       <  4:188687/149‑1 (MQ=255)
                                                                                                                                          cAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCttt                                                    <  4:232561/103‑1 (MQ=255)
                                                                                                                                           aaTAGGCAACGCTGGGTGAGGCAAAGACAAAAACCTAAAATCCTAGTAAAGTTCATACCCCTCCGTTTACAGCAAAAGGATCCCTTCCATGACACACACTTTTTATATCGCCTTATTCTGCGCAACCTGCATATTCTGTGAAGCCTTg      <  1:146364/148‑1 (MQ=255)
                                                                                                                                           aaTAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTGGATATCg                                            >  4:261193/1‑110 (MQ=255)
                                                                                                                                           aaTAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTGGATATCg                                            <  3:261193/110‑1 (MQ=255)
                                                                                                                                               ggCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTGGATATCGCCGTCGTCGGCGCCACCGGCAGCGTCGGTGAAGCCTTGGTGc  <  8:40587/148‑1 (MQ=255)
                                                                                                                                                gCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTGGATATCGCCGTCGTCGGCGCCACCGGCAGCGTCGGTGAAGCCTTGGTg   >  3:49938/1‑146 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGAGTAAAGTGCCGCCCCCGCCGTTTACAGCAAAAGGATCCCTTCCATGACACACCCTTTGGATATCGCCGTCGTCGGCGCCACCGGCAGCGTCGGTGAAGCCTTGGTGC  >  NC_002947/2258110‑2258400

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