Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 50,988 +CG coding (235/843 nt) apaH ← diadenosine tetraphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091350,9851.G92.6% 187.6 / 6.7 54coding (238/843 nt)apaHdiadenosine tetraphosphatase
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base G (25/25);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_00091350,9852.C92.6% 186.7 / 6.8 54coding (238/843 nt)apaHdiadenosine tetraphosphatase
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base C (25/25);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGCGTGATCCCTGCGTGGGCCATCACCAGCTTTTTCTCTTCGTCGATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATT‑‑GCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGATTTCACATAGCGCAGAACATCCAGCGAACCCGGCCCGCGCGCGACCAG  >  NC_000913/50847‑51111
                                                                                                                                           ||                                                                                                                                
gcggcgTGATCCCTGCGTGGGCCATCACCAGCTTTTTCTCTTCGTCGATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATT‑‑g                                                                                                                                 <  2:249864/140‑1 (MQ=255)
 cggcgTGATCCCTGCGTGGGCCATCACCAGCTTTTTCTCTTCGTCGATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATT‑‑Gc                                                                                                                                >  2:709766/1‑140 (MQ=255)
    cgTGATCCCTGCGTGGGCCATCACCAGCTTTTTCTCTTCGTCGATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATT‑‑GCgcg                                                                                                                             >  2:878395/1‑138 (MQ=255)
                gTGGGCCATCACCAGCTTTTTCTCTTCGTCGATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGc                                                                                                                 >  1:664798/1‑140 (MQ=255)
                   ggCCATCACCAGCTTTTTCTCTTCGTCGATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCaaa                                                                                                              >  1:434609/1‑140 (MQ=255)
                   ggCCATCACCAGCTTTTTCTCTTCGTCGATTTGCAGCAGAGGCTGGCACCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCaaa                                                                                                              >  2:401778/1‑140 (MQ=255)
                       aTCACCAGCTTTTTCTCTTCGTCGATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATAcc                                                                                                          >  2:570448/1‑140 (MQ=255)
                       aTCACCAGCTTTTTCTCTTCGTCGATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATAcc                                                                                                          >  2:866448/1‑140 (MQ=255)
                        tCACCAGCTTTTTCTCTTCGTCGATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCg                                                                                                         <  1:868261/140‑1 (MQ=255)
                          aCCAGCTTTTTCTCTTCGTCGATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGcc                                                                                                       <  2:754767/140‑1 (MQ=255)
                            cAGCTTTTTCTCTTCGTCGATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCcag                                                                                                     <  1:709766/140‑1 (MQ=255)
                            cAGCTTTTTCTCTTCGTCGATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCcag                                                                                                     <  1:878395/140‑1 (MQ=255)
                            cAGCTTTTTCTCTTCGTCGATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCcag                                                                                                     <  1:563066/140‑1 (MQ=255)
                               cTTTTTCTCTTCGTCGATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCcagcag                                                                                                  >  2:232327/1‑140 (MQ=255)
                                 ttttCTCTTCGTCGATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCcag                                                                                                     <  1:432274/135‑1 (MQ=255)
                                 ttttCTCTTCGTCGATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCcag                                                                                                     >  2:432274/1‑135 (MQ=255)
                                       ctTCGTCGATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAgcagatgcagat                                                                                          <  2:434609/140‑1 (MQ=255)
                                         tcgtcgATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCcagca                                                                                                   >  1:1000883/1‑129 (MQ=255)
                                         tcgtcgATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCcagca                                                                                                   <  2:1000883/129‑1 (MQ=255)
                                                tttGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATT‑‑GCGCTGCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGc                                                                                 >  1:384036/1‑140 (MQ=255)
                                                 ttGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGcc                                                                                <  1:866448/140‑1 (MQ=255)
                                                  tGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGccc                                                                               <  2:761603/140‑1 (MQ=255)
                                                  tGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATT‑‑GCGCTGCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGccc                                                                               <  2:384036/140‑1 (MQ=255)
                                                          agGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAg                                                                       >  2:924817/1‑140 (MQ=255)
                                                          agGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAg                                                                       >  1:1063390/1‑140 (MQ=255)
                                                            gCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGAc                                                                     >  1:183055/1‑140 (MQ=255)
                                                            gCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGAc                                                                     <  2:53824/140‑1 (MQ=255)
                                                              tGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGt                                                                   <  2:8738/140‑1 (MQ=255)
                                                                 cgccgcAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACg                                                                >  1:770056/1‑140 (MQ=255)
                                                                   gggcagcaAGTTAAGCAGCCCGGCGGCGTCGGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCt                                                              <  2:997352/138‑1 (MQ=255)
                                                                   ccgcAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCt                                                              <  2:448081/140‑1 (MQ=255)
                                                                     gcAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGt                                                            <  1:232327/140‑1 (MQ=255)
                                                                      cAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTc                                                           >  1:26122/1‑140 (MQ=255)
                                                                          cAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCt                                                       <  2:981859/140‑1 (MQ=255)
                                                                            gTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGTAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTaa                                                     <  1:570448/140‑1 (MQ=255)
                                                                            gTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTaa                                                     <  2:1063390/140‑1 (MQ=255)
                                                                                 gcagcTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGAtt                                                <  1:581604/140‑1 (MQ=255)
                                                                                  cagcTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGACCGTGATTGCCCCGCACCAGACGTACGCGGTTGTCTAATGaact                                               >  1:1059375/1‑137 (MQ=255)
                                                                                   agcTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGATTTc                                              >  1:908183/1‑140 (MQ=255)
                                                                                      tcgtcgGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGATTTcaca                                           <  1:97379/140‑1 (MQ=255)
                                                                                       cgtcgGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGATTTCACAt                                          >  2:867535/1‑140 (MQ=255)
                                                                                         tcgGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGATTTCACATAg                                        <  1:637231/140‑1 (MQ=255)
                                                                                          cgGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGATTTCACATAgc                                       <  1:924220/140‑1 (MQ=255)
                                                                                          cgGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCACAGCACCAGACGTACGCTGTCGCCTAAGGATTTCACATAgc                                       >  2:686254/1‑140 (MQ=255)
                                                                                               tCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGATTTCACATAGCGCAGa                                  >  2:904678/1‑140 (MQ=255)
                                                                                                  ggCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGATTTCACATAGCGCAGAACa                               >  2:744865/1‑140 (MQ=255)
                                                                                                         ccAGCAGCGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGATTTCACATAGCGCAGAACATCCAGCg                        >  2:24994/1‑140 (MQ=255)
                                                                                                             cagcGGTGTCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGATTTCACATAGCGCAGAACATCCAGCGAAcc                    <  1:867535/140‑1 (MQ=255)
                                                                                                                  gtgtCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGATTTCACATAGCGCAGAACATCCAGCGAACCCGGcc               >  2:402662/1‑140 (MQ=255)
                                                                                                                    gtCAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGATTTCACATAGCGCAGAACATCCAGCGAACCCGGCCcg             <  2:67734/140‑1 (MQ=255)
                                                                                                                      cAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATCCCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGATTTCACATAGCGCAGAACATCCAGCGAACCCGGCCcgcg           <  2:5733/140‑1 (MQ=255)
                                                                                                                      cAGGCGATCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGATTTCACATAGCGCAGAACATCCAGCGAACCCGGCCcgcg           >  2:201500/1‑140 (MQ=255)
                                                                                                                            aTCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGATTTCACATAGCGCAGAACATCCAGCGAACCCGGCCCGCGCGCGAc     >  1:220308/1‑140 (MQ=255)
                                                                                                                            aTCTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGGGATTGCCCAGCACCATACGTACGCTGTCGCCTAAGGATTTCACATAGCGCAGAACATCCAGCTAACCCTGCCCGCGCGCGAc     >  2:986193/1‑140 (MQ=255)
                                                                                                                              cTTTCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGATTTCACATAGCGCAGAACATCCAGCGAACCCGGCCCGCGCGCGACCa   <  2:26122/140‑1 (MQ=255)
                                                                                                                               tttCGGTTTATTGCGCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGATTTCACATAGCGCAGAACATCCAGCGAACCCGGCCCGCGCGCGACCAg  >  1:938347/1‑140 (MQ=255)
                                                                                                                                           ||                                                                                                                                
GCGGCGTGATCCCTGCGTGGGCCATCACCAGCTTTTTCTCTTCGTCGATTTGCAGCAGAGGCTGGCGCCGCAGCCAGTTAAGCAGCTCGTCGGCATCCGGCGCTTCCAGCAGCGGTGTCAGGCGATCTTTCGGTTTATT‑‑GCG‑‑GCTGATCCCGGCAAATACCGCCAGCAGATGCAGATCGTGATTGCCCAGCACCAGACGTACGCTGTCGCCTAAGGATTTCACATAGCGCAGAACATCCAGCGAACCCGGCCCGCGCGCGACCAG  >  NC_000913/50847‑51111

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