Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,224,100 C→T 7.0% W149* (TGG→TGA PP_1962 ← site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,224,1000CT 7.0% 164.0 / 2.3 57W149* (TGG→TGAPP_1962site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base C (24/29);  new base T (2/2);  total (26/31)
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.78e-01

GAAAGAGGAACCCACGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGG  >  NC_002947/2223963‑2224244
                                                                                                                                         |                                                                                                                                                 
gAAAGAGGAACCCACGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCaa                                                                                                                                        <  7:147289/149‑1 (MQ=255)
 aaaGAGGAACCCACGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGTCAAAATCGCaa                                                                                                                                        >  4:165085/1‑148 (MQ=255)
 aaaGAGGAACCCACGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGTCAAAATCGCaa                                                                                                                                        >  4:159601/1‑148 (MQ=255)
    gagGAACCCACGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACAt                                                                                                                                    >  5:12903/1‑148 (MQ=255)
       gAACCCACGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGaa                                                                                                                                 >  7:21708/1‑148 (MQ=255)
          cccACGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGaaaaa                                                                                                                              >  4:273998/1‑148 (MQ=255)
            cACGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAg                                                                                                                           <  8:179519/149‑1 (MQ=255)
             aCGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAgg                                                                                                                          <  6:12903/149‑1 (MQ=255)
             aCGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAgg                                                                                                                          <  4:175501/149‑1 (MQ=255)
                 ttGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCTAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGTGCGATAAGGATAGCCAAAATCGCAACATGAAAAACAGGCCTg                                                                                                                      <  7:102042/149‑1 (MQ=255)
                    aaTCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGTCAAAATCGCAACATGAAAAACAGGCCTGaa                                                                                                                    <  3:159601/148‑1 (MQ=255)
                    aaTCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGTCAAAATCGCAACATGAAAAACAGGCCTGaa                                                                                                                    <  3:165085/148‑1 (MQ=255)
                        gAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGaataat                                                                                                                <  8:21708/148‑1 (MQ=255)
                        gAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGaataat                                                                                                                >  7:19232/1‑148 (MQ=255)
                         aGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGa                                                                                                              <  3:273998/149‑1 (MQ=255)
                           cccTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACt                                                                                                            >  4:336450/1‑149 (MQ=255)
                             cTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCg                                                                                                          >  6:99110/1‑149 (MQ=255)
                                 aaaTTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGc                                                                                                                         >  6:36386/1‑130 (MQ=255)
                                 aaaTTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGc                                                                                                                         <  5:36386/130‑1 (MQ=255)
                                               gTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTg                                                                                        >  5:104160/1‑149 (MQ=255)
                                                 ccTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGTTAAGGAAAGCCAAAATCGCAACTTTGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCtttt                                                                                         <  1:255050/147‑1 (MQ=255)
                                                   ttGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGtt                                                                                                       >  2:8198/1‑130 (MQ=255)
                                                   ttGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGtt                                                                                                       <  1:8198/130‑1 (MQ=255)
                                                     gtgGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGTAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGACAACTTGCTTTTGcgtcg                                                                                   >  7:18732/1‑148 (MQ=255)
                                                      tgGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCgtcgtc                                                                                 >  4:138375/1‑149 (MQ=255)
                                                        gcgcTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGa                                                                               <  3:138375/149‑1 (MQ=255)
                                                                    aTCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGAtt                                                                     >  8:302800/1‑147 (MQ=255)
                                                                     tCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTaa                                                                   <  5:99110/148‑1 (MQ=255)
                                                                          tgatgTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCg                                                             >  3:223645/1‑149 (MQ=255)
                                                                          tgatgTGATATGTTCGGGTAGACTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAATAATGCCTGAATAATGAATCGGTTGTACACTTGCATGTGCTTATTCGAGTGAGCTATTAATCATCg                                                             >  8:335796/1‑149 (MQ=255)
                                                                               tgATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTc                                                                                                    >  8:127287/1‑105 (MQ=255)
                                                                               tgATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTc                                                                                                    <  7:127287/105‑1 (MQ=255)
                                                                               tgATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGAtt                                                         <  7:250249/148‑1 (MQ=255)
                                                                                gATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTg                                                       >  3:232199/1‑149 (MQ=255)
                                                                                 atttGTTCGGGTAGCCTTAAACTTTCAGGAATGGCTTGTTGAGCGATAAGGATAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGt                                                      <  4:94969/145‑1 (MQ=255)
                                                                                   atGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGtt                                                     <  8:19232/148‑1 (MQ=255)
                                                                                    tGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTa                                                    <  4:223645/148‑1 (MQ=255)
                                                                                    tGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTa                                                    <  7:150511/148‑1 (MQ=255)
                                                                                      ttCGGGTAGCCTTAAACTTTCAGGAATTAGGTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATga                                                 <  7:191696/149‑1 (MQ=255)
                                                                                         gggTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGAt                                              <  7:302800/149‑1 (MQ=255)
                                                                                                cTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGt                                        <  3:336450/148‑1 (MQ=255)
                                                                                                 ttAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGt                                                      <  3:315482/133‑1 (MQ=255)
                                                                                                 ttAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGt                                                      >  4:315482/1‑133 (MQ=255)
                                                                                                 ttAAACGTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGGTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGc                                      >  4:315561/1‑149 (MQ=255)
                                                                                                     aCTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAg                                  <  4:232199/149‑1 (MQ=255)
                                                                                                       tttCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAgg                                 >  1:173754/1‑148 (MQ=255)
                                                                                                           aGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCg                            <  6:104160/149‑1 (MQ=255)
                                                                                                             gAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCt                          >  4:176873/1‑149 (MQ=255)
                                                                                                                   gCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGTGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCa                                   <  5:263862/134‑1 (MQ=255)
                                                                                                                       gttgAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCAc                                                                                                 <  1:111511/68‑1 (MQ=255)
                                                                                                                       gttgAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCAc                                                                                                 <  1:115461/68‑1 (MQ=255)
                                                                                                                       gttgAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCAc                                                                                                 >  2:111511/1‑68 (MQ=255)
                                                                                                                       gttgAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCAc                                                                                                 >  2:115461/1‑68 (MQ=255)
                                                                                                                                   ggAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCaa      >  5:158191/1‑147 (MQ=255)
                                                                                                                                    gAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAg   >  3:278812/1‑149 (MQ=255)
                                                                                                                                     aaaGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAgg  <  3:176873/149‑1 (MQ=255)
                                                                                                                                     aaaGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGCCCCAATAgg  <  3:315561/149‑1 (MQ=255)
                                                                                                                                         ccAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTggg            <  8:117124/135‑1 (MQ=255)
                                                                                                                                         ccAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTggg            >  7:117124/1‑135 (MQ=255)
                                                                                                                                         |                                                                                                                                                 
GAAAGAGGAACCCACGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGG  >  NC_002947/2223963‑2224244

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: