Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 861,896 C→T 6.6% N459N (AAC→AAT pgm → bisphosphoglycerate mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730861,8960CT 6.6% 178.5 / 2.9 61N459N (AAC→AATpgmbisphosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base T (2/2);  total (31/30)
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.93e-01

GAAGCGGTTGATGAATGTTTAGGAGAAGTGGTTGATAAGATTTTAGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTAGAACAACCTGAAGATATGACAGGTGAATCTTTAATT  >  CP000730/861756‑862028
                                                                                                                                            |                                                                                                                                    
gAAGCGGTTGATGAATGTTTAGGAGAAGTGGTTGATAAGATTTTAGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAAACCATGACTACGCATACAACGAAc                                                                                                                                      <  1:565895/141‑1 (MQ=255)
 aaGCGGTTGATGAATGTTTAGGAGAAGTGGTTGATAAGATTTTAGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAAcc                                                                                                                                     >  1:541269/1‑141 (MQ=255)
    cGGTTGATGAATGTTTAGGAGAAGTGGTTGATAAGATTTTAGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAATCCAg                                                                                                                                  >  1:404449/1‑141 (MQ=255)
    cGGTTGATGAATGTTTAGGAGAAGTGGTTGATAAGATTTTAGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAATCCAg                                                                                                                                  >  1:404445/1‑141 (MQ=255)
    cGGTTGATGAATGTTTAGGAGAAGTGGTTGATAAGATTTTAGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAg                                                                                                                                  >  1:166220/1‑141 (MQ=255)
      gTTGATGAATGTTTAGGAGAAGTGGTTGATAAGATTTTAGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTa                                                                                                                                <  1:511772/141‑1 (MQ=255)
          atgaATGTTTATGTGAAGCGTTTTATAAGATTTTAGATCTGGACGGTTATGCAATTCTTACTGCTGACCATGGTAACTCTCATCAACTGTTGACGGATGATGATCAACCAATGACTACGCACACAACGAACCCAGTACCAg                                                                                                                            <  2:250828/141‑1 (MQ=255)
                  ttAGGAGAAGTGGTGGATAAGATTTTAGACATGGACGGTTATCCAATTATTACTGCTCACCATCGTAATGCTGATCAAGTATTGCCGGCTCACGATCAACCAATGACTACGCATCCAACGAACCCAGTACCAGTGATTGTa                                                                                                                    <  2:143949/141‑1 (MQ=255)
                     ggAGAAGTGGTTGATAAGATTTTAGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCCGTGATTGTAACa                                                                                                                 >  1:494630/1‑141 (MQ=255)
                         aaGTGGTTGATAAGATTTTAGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAaag                                                                                                             >  1:462374/1‑141 (MQ=255)
                           tttgTTTATAAGCTTTTAGACATGGACGGTTATGCAATTATTACTGCTGATCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAaagaa                                                                                                           <  2:298335/138‑1 (MQ=255)
                            tGGTTGATAAGATTTAAGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAaagaag                                                                                                          <  1:509127/141‑1 (MQ=255)
                                 gATAAGATTTTAGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGtt                                                                                                     <  1:361083/141‑1 (MQ=255)
                                  aTAAGATTTTAGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCCGCACAAGTGATTGTAACAAAAGAAGGCGtaa                                                                                                    >  2:130499/1‑139 (MQ=255)
                                        ttttAGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGGGTTACACATa                                                                                              >  2:68787/1‑141 (MQ=255)
                                        ttttAGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTa                                                                                              <  2:541269/141‑1 (MQ=255)
                                        ttttAGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTa                                                                                              <  1:145372/141‑1 (MQ=255)
                                          ttAGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTaga                                                                                            >  1:411534/1‑141 (MQ=255)
                                           tAGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTagag                                                                                           <  2:396572/141‑1 (MQ=255)
                                            aGACATGGACGGTTGTGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAgaga                                                                                          >  2:387387/1‑141 (MQ=255)
                                            aGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAgaga                                                                                          >  1:147686/1‑141 (MQ=255)
                                             gACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGaa                                                                                         >  2:130961/1‑141 (MQ=255)
                                               cATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAAc                                                                                       >  1:406482/1‑141 (MQ=255)
                                               cATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAAc                                                                                       >  2:422860/1‑141 (MQ=255)
                                                  ggACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTgg                                                                                    >  1:175707/1‑141 (MQ=255)
                                                     cGGTTATGCAATTATTACTGCTGACTATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCg                                                                                 <  2:10675/141‑1 (MQ=255)
                                                           tGCATTTATTACTGCTGACCATGGTAACTCTGATCAAGTATGGACGGATTATGATCAACTACTGACAAAGCATACAACGAACCCAGTACCCGGCAGAGTCACAAAAGAAGGCGTTACACTTAAAGCAACTGGTCGCTTAgg                                                                           >  2:225373/1‑141 (MQ=255)
                                                            gCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAggg                                                                          >  2:329471/1‑140 (MQ=255)
                                                            gCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGt                                                                          <  1:144403/141‑1 (MQ=255)
                                                                ttattaCTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACt                                                                      <  1:245097/141‑1 (MQ=255)
                                                                       tgctgACCATGGTATCTCTGCTCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCGTAGGTGACTTAGCAcc                                                               <  1:284815/141‑1 (MQ=255)
                                                                       tgctgACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCAcc                                                               >  2:382219/1‑141 (MQ=255)
                                                                           gACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACa                                                           <  1:130961/141‑1 (MQ=255)
                                                                             ccATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACatt                                                         <  1:318067/141‑1 (MQ=255)
                                                                                tGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCATCTACattatt                                                      <  1:414155/141‑1 (MQ=255)
                                                                                tGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACattatt                                                      >  1:122817/1‑141 (MQ=255)
                                                                                tGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACattatt                                                      <  1:381383/141‑1 (MQ=255)
                                                                                     aCTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGAtt                                                 >  2:294206/1‑141 (MQ=255)
                                                                                       tctGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAATCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACGCTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATtta                                               <  2:404445/141‑1 (MQ=255)
                                                                                       tctGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAATCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATtta                                               <  2:404449/141‑1 (MQ=255)
                                                                                       tctGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAGCTGGTCGCTTAGGTGACTACGCACCTACATTATTAGAttaa                                               >  2:163160/1‑139 (MQ=255)
                                                                                       tctGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATtta                                               >  2:421939/1‑141 (MQ=255)
                                                                                                gTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGCTTTATTAAATGTa                                      >  1:53661/1‑141 (MQ=255)
                                                                                                  aTTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGATAAACCGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTAGa                                    <  1:474067/141‑1 (MQ=255)
                                                                                                  aTTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTAGa                                    >  1:561126/1‑141 (MQ=255)
                                                                                                     gACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTAGaaca                                 <  1:374492/141‑1 (MQ=255)
                                                                                                     gACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTAGaaca                                 >  1:312751/1‑141 (MQ=255)
                                                                                                     gACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGAATTAGCACCTACATTATTGGCTTTATTGAATGTGGaaac                                 >  1:490893/1‑139 (MQ=255)
                                                                                                      aCGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTAGaacaa                                >  2:383151/1‑141 (MQ=255)
                                                                                                       cGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTAGaacaac                               <  1:446358/141‑1 (MQ=255)
                                                                                                        gcatgatgatgAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGGAACAAAAGAAGGCGTTCCACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTAGAACAAcc                              <  2:315504/139‑1 (MQ=255)
                                                                                                          atgatgatCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTAGAACAACCTg                            >  2:240662/1‑141 (MQ=255)
                                                                                                           tgatgaaCAACCAATGTCTACGCATACAACGAACCCAGTACCAGTGATTGTGTCAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTAGAACAACCTGa                           <  1:130499/141‑1 (MQ=255)
                                                                                                                 tCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTAGAACAACCTGAAGatat                     <  1:363425/141‑1 (MQ=255)
                                                                                                                     ccAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTAGAACAACCTGAAGATATGACa                 >  2:162375/1‑141 (MQ=255)
                                                                                                                      cAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTGACAAAAGAAGGCGTTACGCTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTGGAACAACCTGAAGATATGACAg                >  2:44239/1‑141 (MQ=255)
                                                                                                                      cAAGGACTACGAAAACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTAGAACAACCTGAAGATATGACAg                <  1:68787/141‑1 (MQ=255)
                                                                                                                          gACTACGCATACAACGAACCCAGTACCAGTGATTGTAACTAAAGAAGGCGTTACACTTAGAGAAAATGGTCGCTTAGGTGACTTCGCACCTACATTATTAGATTTATTAAATGTAGGACAACCTGAAGATATGACAGGTGa            >  2:282719/1‑141 (MQ=255)
                                                                                                                          gACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTAGGACAACGTGACGATCTGAGAGGTGa            >  1:162324/1‑141 (MQ=255)
                                                                                                                              aCGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTAGAACAACCTGAAGATATGACAGGTGAATCt        <  1:489481/141‑1 (MQ=255)
                                                                                                                                gCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTAGAACAACCTGAAGATATGACAGGTGAATCttt      <  1:496029/141‑1 (MQ=255)
                                                                                                                                   tACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTAGAACAACCTGAAGATATGACAGGTGAATCTttaat   <  1:382219/141‑1 (MQ=255)
                                                                                                                                    aataCGAACCCAGTACCAGTGATTGTAACAAAAGATCGCGATACGCGGAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGGAGAACAACCGGAAGATATGACAGGTGAATCTttaatt  <  1:13961/138‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                    
GAAGCGGTTGATGAATGTTTAGGAGAAGTGGTTGATAAGATTTTAGACATGGACGGTTATGCAATTATTACTGCTGACCATGGTAACTCTGATCAAGTATTGACGGATGATGATCAACCAATGACTACGCATACAACGAACCCAGTACCAGTGATTGTAACAAAAGAAGGCGTTACACTTAGAGAAACTGGTCGCTTAGGTGACTTAGCACCTACATTATTAGATTTATTAAATGTAGAACAACCTGAAGATATGACAGGTGAATCTTTAATT  >  CP000730/861756‑862028

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: