Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,359,877 C→A 16.1% Q344H (CAG→CAT) ‡ puuA ← glutamate‑‑putrescine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,359,8770CA16.1% 172.1 / 21.8 56Q344H (CAG→CAT) ‡puuAglutamate‑‑putrescine ligase
Reads supporting (aligned to +/- strand):  ref base C (20/27);  new base A (3/6);  total (23/33)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

TGCAAAATACCGGCAAAAATCGCTGCCATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCTTCGCCTTCCGCGTC  >  NC_000913/1359734‑1360017
                                                                                                                                               |                                                                                                                                            
tGCAAAATACCGGCAAAAATCGCTGCCATCACCAGATATGTGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGt                                                                                                                                         >  4:340086/1‑149 (MQ=255)
 gCAAAATACCGGCAAAAATCGCTGCCATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTc                                                                                                                                        >  4:119269/1‑149 (MQ=255)
   aaaaTACCGGCAAAAATCGCTGCCATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCgg                                                                                                                                      <  4:182543/149‑1 (MQ=255)
     aaTACCGGCAAAAATCGCTGCCATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCAGGCGTCgtgg                                                                                                                                    <  4:224032/149‑4 (MQ=255)
     aaTACCGGCAAAAATCGCTGCCATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCAGGCGTCgtgg                                                                                                                                    <  4:225075/149‑4 (MQ=255)
     aaTACCGGCAAAAATCGCTGCCATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCAGGCGTCgtgg                                                                                                                                    <  2:45930/149‑4 (MQ=255)
      aTACCGGCAAAAATCGCTGCCATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGcac                                                                                                                                   >  4:527015/1‑149 (MQ=255)
           ggCAAAAATCGCTGCCATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCAGGCGTCgtgg                                                                                                                                    >  3:528825/1‑140 (MQ=255)
           ggCAAAAATCGCTGCCATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCAGGCGTCgtgg                                                                                                                                    <  4:528825/143‑4 (MQ=255)
            gCAAAAATCGCTGCCATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGCAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCacatacat                                                                                                                             >  2:558911/1‑149 (MQ=255)
              aaaaaTCGCTGCCATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTc                                                                                                                           <  3:355152/149‑1 (MQ=255)
                             cACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAgcggc                                                                                                            <  3:221454/149‑1 (MQ=255)
                              aCCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCAGGCGTCgtgg                                                                                                                                    >  4:469065/1‑121 (MQ=255)
                              aCCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCAGGCGTCgtgg                                                                                                                                    <  3:469065/124‑4 (MQ=255)
                               ccAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGa                                                                                                          >  4:188962/1‑149 (MQ=255)
                                 aGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCAGGCGTCgtg                                                                                                                                     <  4:286861/120‑3 (MQ=255)
                                 aGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCAGGCGTCgtg                                                                                                                                     >  3:286861/1‑118 (MQ=255)
                                  gATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATAc                                                                                                       <  2:199529/149‑1 (MQ=255)
                                   atatGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCTGCTGGAAGCGGCGATACg                                                                                                      >  2:9013/1‑149 (MQ=255)
                                     atGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAg                                                                                                    <  2:323918/149‑1 (MQ=255)
                                               aTCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTg                                                                                          >  2:211007/1‑149 (MQ=255)
                                               aTCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTg                                                                                          <  4:365959/149‑1 (MQ=255)
                                                tCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGa                                                                                                          >  4:425556/1‑132 (MQ=255)
                                                tCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGa                                                                                                          <  3:425556/132‑1 (MQ=255)
                                                tCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTgg                                                                                         <  1:128655/149‑1 (MQ=255)
                                                     aCCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCa                                                                                    <  1:229141/149‑1 (MQ=255)
                                                     aCCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCa                                                                                    <  3:188962/149‑1 (MQ=255)
                                                       ccGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCt                                                                                                                     >  1:160853/1‑114 (MQ=255)
                                                       ccGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCt                                                                                                                     <  2:160853/114‑1 (MQ=255)
                                                        cGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGttt                                                                                           <  1:577084/139‑1 (MQ=255)
                                                        cGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGttt                                                                                           >  2:577084/1‑139 (MQ=255)
                                                         gCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCt                                                                                                                                             <  1:519223/88‑1 (MQ=255)
                                                         gCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCt                                                                                                                                             >  2:519223/1‑88 (MQ=255)
                                                               cgATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCAt                                                                          >  2:486392/1‑149 (MQ=255)
                                                                    ttCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATcgacg                                                                     <  3:527015/149‑1 (MQ=255)
                                                                             gTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCa                                                            <  2:456891/149‑1 (MQ=255)
                                                                             gTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTCTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCa                                                            <  1:459246/149‑1 (MQ=255)
                                                                                   aTGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGATTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCaa                                                      >  3:512719/1‑149 (MQ=255)
                                                                                    tGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAat                                                     <  2:65494/149‑1 (MQ=255)
                                                                                        ggTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAatcatc                                                 <  3:119269/149‑1 (MQ=255)
                                                                                         gTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATcc                                                <  3:30532/149‑1 (MQ=255)
                                                                                               gcACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGCGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGa                                          <  1:486392/149‑1 (MQ=255)
                                                                                                    gAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATc                                     >  4:205355/1‑149 (MQ=255)
                                                                                                          gCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTc                               >  2:505228/1‑149 (MQ=255)
                                                                                                          gCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTc                               <  4:258328/149‑1 (MQ=255)
                                                                                                          gCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTc                               >  4:405391/1‑149 (MQ=255)
                                                                                                          gCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTc                               >  3:310472/1‑149 (MQ=255)
                                                                                                           cAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTca                              <  2:161734/149‑1 (MQ=255)
                                                                                                            aGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTcag                             >  3:422455/1‑149 (MQ=255)
                                                                                                                           ttgttATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCTTc              <  3:163649/149‑1 (MQ=255)
                                                                                                                           ttgttATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCTTc              <  3:223312/149‑1 (MQ=255)
                                                                                                                           ttgttATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCTTc              <  4:416668/149‑1 (MQ=255)
                                                                                                                             gttATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCTTCGc            <  4:209821/149‑1 (MQ=255)
                                                                                                                                aTGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCTTCGCCtt         >  4:402040/1‑149 (MQ=255)
                                                                                                                                      ccACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCTTCGCCTTCCGCGt   >  4:124899/1‑149 (MQ=255)
                                                                                                                                       cACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGTCAGCTACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCTTCGCCTTCCGCGTc  <  3:424848/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                            
TGCAAAATACCGGCAAAAATCGCTGCCATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCTTCGCCTTCCGCGTC  >  NC_000913/1359734‑1360017

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