Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 992,865 A→T 9.0% intergenic (+23/‑129) der → / → PP_0858 GTPase Der/methionine/glutamine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947992,8650AT 9.0% 140.6 / 3.6 52intergenic (+23/‑129)der/PP_0858GTPase Der/methionine/glutamine aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (24/23);  new base T (2/3);  total (26/26)
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.73e-01

CAAGGGCGGCGAGAACCCGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAAGAGGGCTCCATGATCCG  >  NC_002947/992719‑993001
                                                                                                                                                   |                                                                                                                                        
cAAGGGCGGCGAGAACCCGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAg                                                                                                                                        <  2:80907/149‑1 (MQ=255)
   gggCGGCGAGAACCCGTACGAGGGTAAGAAGAACACCCTTACCGATCCCCAGGTCAACAAGAAACGC‑CGCCTGATGGCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGAcc                                                                                                                                     <  3:54026/149‑1 (MQ=255)
                ccGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACg                                                                                                                        >  5:238437/1‑149 (MQ=255)
                      gAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCAAAGCAGCGCCATCTGGGGGCGGGCAAGACCGCCCCCACAGACCCCCttt                                                                                                                  >  7:79639/1‑149 (MQ=255)
                        gggTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGa                                                                                                                <  6:257513/149‑1 (MQ=255)
                         ggTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCGCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTg                                                                                                                 >  1:256298/1‑148 (MQ=255)
                           tAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATcc                                                                                                             <  6:238437/149‑1 (MQ=255)
                             agaagaACACCCTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGTCTTGCCCGCCCCca                                                                                                                              >  6:302257/1‑130 (MQ=255)
                             agaagaACACCCTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGTCTTGCCCGCCCCca                                                                                                                              <  5:302257/130‑1 (MQ=255)
                              gaagaaCACCCTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCaga                                                                                                          <  5:232652/149‑1 (MQ=255)
                               aagaaCACCCTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCagag                                                                                                         >  4:178589/1‑149 (MQ=255)
                                  aaCACCCTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCaa                                                                                                      >  8:27884/1‑149 (MQ=255)
                                  aaCACCCTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCaa                                                                                                      <  7:27884/149‑1 (MQ=255)
                                    cacCCTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCtt                                                                                                                   >  2:36336/1‑134 (MQ=255)
                                    cacCCTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCtt                                                                                                                   <  1:36336/134‑1 (MQ=255)
                                    cacCCTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCt                                                                                                                    <  5:166087/133‑1 (MQ=255)
                                    cacCCTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCt                                                                                                                    >  6:166087/1‑133 (MQ=255)
                                       ccTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATcccc                                                                                                 <  2:256298/149‑1 (MQ=255)
                                       ccTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATcccc                                                                                                 >  5:185385/1‑149 (MQ=255)
                                         ttACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATccccc                                                                                                >  3:97257/1‑148 (MQ=255)
                                                          aacaaGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGTCTTGCCCGCCCCca                                                                                                                              <  5:289479/101‑1 (MQ=39)
                                                          aacaaGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGTCTTGCCCGCCCCca                                                                                                                              >  6:289479/1‑101 (MQ=39)
                                                            caaGAACCGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACACACGCTACGCAGCGCCATGTGGGGGCGGGCAAGCCCGCCCCCCCAGCCGCCCTGTGCTCCCGAGCAATCCCCCGATGGTTTCACCCGGTTTcc                                                                            >  5:289485/1‑149 (MQ=255)
                                                            caaGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTcc                                                                            >  5:277252/1‑149 (MQ=255)
                                                             aaGAAACGC‑CGCCTGATGTCGCACCACACGAAGGCCGAGAAGAAGCGGCGCGACAAACGCTACGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCt                                                                           <  6:289485/149‑1 (MQ=255)
                                                                        gccTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCCCAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGa                                                                 >  5:268754/1‑149 (MQ=255)
                                                                         ccTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGAt                                                                >  1:47383/1‑149 (MQ=255)
                                                                          cTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCTGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTg                                                                          <  2:91158/138‑1 (MQ=255)
                                                                             aTGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGc                                                            <  6:185385/149‑1 (MQ=255)
                                                                              tGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGc                                                            >  6:245264/1‑148 (MQ=255)
                                                                                   caccacAAGAAGGCCGACAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGCGGCGGGCAACACCCCCCCCACAGACGCCCTTTGATCCAGACCAATCCCCCTATTGTTTCAACCATTATCCTCACCGCTCGATCCGCTATGCt                                                      >  2:164547/1‑149 (MQ=255)
                                                                                      cacaAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTc                                                                             <  1:191532/123‑1 (MQ=255)
                                                                                      cacaAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTc                                                                             >  2:191532/1‑123 (MQ=255)
                                                                                      cacaAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCg                                                    >  1:351199/1‑148 (MQ=255)
                                                                                       acaAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGa                                                  <  6:277252/149‑1 (MQ=255)
                                                                                       acaAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCATATTGTTTCAAc                                                                                    >  4:316562/1‑115 (MQ=255)
                                                                                        caAGAAGGCCGAGAAGAAGCGTAGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGTCTAGAACTCCCACACAGACGCCCTTTGATCCAGAGCAATCGCGCTATTTTTTCAACGTTTTTCCTGACCGTTGGATACTCTTTGCtgtgta                                                 <  8:79639/149‑6 (MQ=255)
                                                                                            aagGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCa                                             <  1:164547/149‑1 (MQ=255)
                                                                                              ggCCGAGAAGTAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGCCGCCCTGTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCAcc                                           >  3:102087/1‑149 (MQ=255)
                                                                                               gCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCAccc                                          >  8:217744/1‑149 (MQ=255)
                                                                                                ccGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCAcccc                                         <  3:178589/149‑1 (MQ=255)
                                                                                                ccGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCAcccc                                         <  8:282433/149‑1 (MQ=255)
                                                                                                 cGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCt                                                           >  6:347497/1‑130 (MQ=255)
                                                                                                 cGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCt                                                           <  5:347497/130‑1 (MQ=255)
                                                                                                 cGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCg                                        >  4:187724/1‑149 (MQ=255)
                                                                                                          gCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCGCAGACGCCCTTTGCTGCAGAGCCATCCCCGTATTGTGTAATCGTTTTTACTGACCGCTCGATCCGATATGATCGGACTCCCCCCCGTGCGGGcta                               >  6:226422/1‑149 (MQ=255)
                                                                                                            gTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGAtata                             <  2:47383/149‑1 (MQ=255)
                                                                                                             tCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATAc                            <  2:351199/149‑1 (MQ=255)
                                                                                                                      aCGATAAGCAGCGCCATGTGGGGGCGGGCAAGCCCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTTTTGTTTCAACCTTTTTCCTGACGCCTCGATCCGCTATGCTCGGACTCCACCACGTGCCGTATATACCCCAGGa                     >  8:245763/1‑147 (MQ=255)
                                                                                                                                   ccATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAAGAGGGCTCCATGa      <  7:217744/149‑1 (MQ=255)
                                                                                                                                     aTGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAAGAGGGCTCCATGATc    <  5:245264/149‑1 (MQ=255)
                                                                                                                                        ggggggCGGGCAAGACCGCACCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCTCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAAGAGGGCTCCATGATCCg  <  4:102087/147‑1 (MQ=255)
                                                                                                                                                   aaGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAagag                <  2:250259/123‑1 (MQ=255)
                                                                                                                                                   aaGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACACCAGGAAagag                >  1:250259/1‑123 (MQ=255)
                                                                                                                                                   |                                                                                                                                        
CAAGGGCGGCGAGAACCCGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGC‑CGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAAGAGGGCTCCATGATCCG  >  NC_002947/992719‑993001

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