Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,359,878 T→G 16.5% Q344P (CAG→CCG) ‡ puuA ← glutamate‑‑putrescine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,359,8780TG16.5% 138.3 / 23.1 54Q344P (CAG→CCG) ‡puuAglutamate‑‑putrescine ligase
Reads supporting (aligned to +/- strand):  ref base T (19/26);  new base G (3/6);  total (22/32)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

TGCAAAATACCGGCAAAAATCGCTGCCATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCTTCGCCTTCCGCGTCAGAAAGCAC  >  NC_000913/1359734‑1360026
                                                                                                                                                |                                                                                                                                                    
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                                        >  3:310472/1‑149 (MQ=255)
                                                                                                          gCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTc                                        <  4:258328/149‑1 (MQ=255)
                                                                                                          gCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTc                                        >  4:405391/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                       <  4:416668/149‑1 (MQ=255)
                                                                                                                           ttgttATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCTTc                       <  3:223312/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)
                                                                                                                                                tGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCTTCGCCTTCCGCGTCAGAAAGCAc  <  3:340086/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TGCAAAATACCGGCAAAAATCGCTGCCATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCTTCGCCTTCCGCGTCAGAAAGCAC  >  NC_000913/1359734‑1360026

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