Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 992,868 A→C 7.9% intergenic (+26/‑126) 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,8680AC 7.9% 142.4 / 4.4 51intergenic (+26/‑126)der/PP_0858GTPase Der/methionine/glutamine aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (21/26);  new base C (2/2);  total (23/28)
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.08e-01

AAGGGCGGCGAGAACCCGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAAGAGGGCTCCATGATCCGCAGCAAACT  >  NC_002947/992720‑993010
                                                                                                                                                    |                                                                                                                                              
aaGGGCGGCGAGAACCCGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGa                                                                                                                                                <  5:361500/149‑1 (MQ=255)
aaGGGCGGCGAGAACCCGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGa                                                                                                                                                <  7:326128/149‑1 (MQ=255)
aaGGGCGGCGAGAACCCGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGa                                                                                                                                                >  2:33369/1‑149 (MQ=255)
     cggcgAGAACCCGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCTCCACAAGACGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGcc                                                                                                                                           <  2:281763/149‑1 (MQ=255)
     aggcgAGAAGCCGTTCGAGTGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCGTGATGTCGCACCACAAGAAGGCCGAGAAGTAGCGTCGCGACAAACGCTATGCTGCGCCATGTGGGGGCGGGCAAGACCGcc                                                                                                                                           <  7:268250/148‑1 (MQ=255)
          agaACCCGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGACGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCcac                                                                                                                                      >  1:47176/1‑149 (MQ=255)
            aaCCCGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCca                                                                                                                                       >  6:188656/1‑146 (MQ=255)
              cccGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGTCTTGCCCGCCCCcacatgg                                                                                                                                  <  2:230750/149‑4 (MQ=255)
              cccGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGcc                                                                                                                                           >  5:357635/1‑140 (MQ=255)
              cccGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGcc                                                                                                                                           <  6:357635/140‑1 (MQ=255)
              cccGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGAc                                                                                                                                  >  1:136975/1‑149 (MQ=255)
               ccGTAAGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACg                                                                                                                                 >  7:238386/1‑149 (MQ=255)
                      aGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTg                                                                                                                          <  5:188656/149‑1 (MQ=255)
                          tAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATc                                                                                                                       >  6:12086/1‑148 (MQ=255)
                            agaagaACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATc                                                                                                                       >  4:281779/1‑146 (MQ=255)
                            agaagaACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCag                                                                                                                    <  5:12086/149‑1 (MQ=255)
                             gaagaaCACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCaga                                                                                                                   <  2:136975/149‑1 (MQ=255)
                             gaagaaCACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCaga                                                                                                                   <  4:183451/149‑1 (MQ=255)
                                 aaCACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGAcc                                                                                                                                              >  2:38266/1‑118 (MQ=255)
                                 aaCACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGAcc                                                                                                                                              <  1:38266/118‑1 (MQ=255)
                                   cacCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAACAAGGCCGAGAAGAAGCTTCGCGACAAACGCTAAGCAGCGCCATCAGGCGGCGGGCAAGACCTCCCCCACACACGCCCTTAGATCCAGAGCAATc                                                                                                             >  8:27126/1‑149 (MQ=255)
                                    acCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATcc                                                                                                            <  8:238386/149‑1 (MQ=255)
                                          aCCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACTGACGCCCTTTGATCCAGAGCAATCCCCCTa                                                                                                       <  1:33369/148‑1 (MQ=255)
                                                         aacaaGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATGGGTTCAACCttttt                                                                                       >  2:111856/1‑149 (MQ=255)
                                                          acaaGAAACGCCGCCTGATGTCGCACCACAAGACGGCCGATAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTc                                                                                      <  7:27126/149‑1 (MQ=255)
                                                           caaGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGCCGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTcc                                                                                     >  5:24426/1‑149 (MQ=255)
                                                                  cgccgcCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGNCCGc                                                                              <  5:65971/149‑1 (MQ=255)
                                                                     cgccTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTc                                                                            >  1:180103/1‑148 (MQ=255)
                                                                            tGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCt                                                                    <  3:281779/149‑1 (MQ=255)
                                                                             gTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTa                                                                   >  2:270610/1‑149 (MQ=255)
                                                                             gTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGAGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCACCACAGACCCCCTTAGATCAAGAGCAATCCCCCTATTGTTTCAACCTTTTTCATGACGGCTCCATCCCCTa                                                                   >  7:272028/1‑149 (MQ=255)
                                                                                 caccacAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGc                                                                >  8:224852/1‑148 (MQ=255)
                                                                                  accacAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTc                                                              <  7:224852/149‑1 (MQ=255)
                                                                                    cacaAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCACGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCg                                                             >  4:214701/1‑148 (MQ=255)
                                                                                       aagaagGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACt                                                         >  4:324839/1‑149 (MQ=255)
                                                                                        agaagGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTc                                                        <  1:270610/149‑1 (MQ=255)
                                                                                        agaagGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTc                                                        <  3:217059/149‑1 (MQ=255)
                                                                                                      aagCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCGCGAGCAATCCCCCTTTTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCg                                                             <  5:259128/130‑1 (MQ=255)
                                                                                                      aagCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGCGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGCTCCGCTATGCACg                                                             >  6:259128/1‑130 (MQ=255)
                                                                                                       agCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGat                                         <  6:24426/149‑1 (MQ=255)
                                                                                                             gcgACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATAccc                                   >  3:285090/1‑149 (MQ=255)
                                                                                                             gcgACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATAccc                                   >  4:355393/1‑149 (MQ=255)
                                                                                                               gACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCa                                 <  2:47176/149‑1 (MQ=255)
                                                                                                               gACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCa                                 <  3:355393/149‑1 (MQ=255)
                                                                                                               caaaaaCGCCAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCa                                 <  7:306331/146‑1 (MQ=255)
                                                                                                                 cAAACGCTAAGCAGCGCCATGTGGTGGCGGGCAAGACCGCCCCGACAGACGCCATTTGATCCAGAGCAATCCCGCTATTGTTTCAACCTTTTTCCTGTCCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAg                                <  7:311415/148‑1 (MQ=255)
                                                                                                                 cAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCt                                                                                    <  7:267903/96‑1 (MQ=255)
                                                                                                                 cAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCt                                                                                    >  8:267903/1‑96 (MQ=255)
                                                                                                                           gcagcGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAAGAGGGCt                     <  3:324839/149‑1 (MQ=255)
                                                                                                                                gcCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCATATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAAGAGGGCTCCATg                <  1:185407/149‑1 (MQ=255)
                                                                                                                                      tgGGGGCGGTCTTGCCCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAagag                         <  8:338993/134‑1 (MQ=255)
                                                                                                                                      tgGGGGCGGTCTTGCCCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAagag                         >  7:338993/1‑134 (MQ=255)
                                                                                                                                        ggggCGGTCTTGCCCGCCCCCACAGACGCCCTTTGATCCAGAGCAATACCCCTATTGTTTCAACCTGTTTCCTGACCGCTCGATCCGCGATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAAGAGGGCTCCATGATCCgcag        >  4:344103/1‑149 (MQ=255)
                                                                                                                                              ggcAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAAGAGGGCTCCATGATCCGCAGCAAACt  <  2:180103/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
AAGGGCGGCGAGAACCCGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAAGAGGGCTCCATGATCCGCAGCAAACT  >  NC_002947/992720‑993010

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: ‑