Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,359,877 C→A 13.0% 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,8770CA13.0% 148.0 / 10.1 46Q344H (CAG→CAT) ‡puuAglutamate‑‑putrescine ligase
Reads supporting (aligned to +/- strand):  ref base C (20/20);  new base A (3/3);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGCAAAAATCGCTGCCATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCTTCGCCTTCCGCGTCAGAAAGC  >  NC_000913/1359743‑1360024
                                                                                                                                      |                                                                                                                                                   
ccGGCAAAAATCGCTGCCATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCacat                                                                                                                                        <  2:447054/148‑1 (MQ=255)
                ccATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGaa                                                                                                                        >  1:210976/1‑148 (MQ=255)
                 cATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAg                                                                                                                       >  4:129511/1‑148 (MQ=255)
                  aTCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCAGGCGTCgtgg                                                                                                                                           >  4:558865/1‑124 (MQ=255)
                  aTCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCAGGCGTCgtgg                                                                                                                                           <  3:558865/127‑4 (MQ=255)
                  aTCACCAGATATGGGTTAGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAgc                                                                                                                      >  2:433536/1‑148 (MQ=255)
                   tCACCAGATATGGGTTAGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAgcg                                                                                                                     <  1:433536/148‑1 (MQ=255)
                    cACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCAGGCGTCgtgg                                                                                                                                           <  3:946732/125‑4 (MQ=255)
                    cACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCAGGCGTCgtgg                                                                                                                                           >  4:946732/1‑122 (MQ=255)
                      ccAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAgcggcg                                                                                                                  <  3:129511/148‑1 (MQ=255)
                         gATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATa                                                                                                               >  4:328591/1‑148 (MQ=255)
                             tGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAg                                                                                                           <  1:167963/148‑1 (MQ=255)
                                     cATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGtt                                                                                                   <  3:328591/148‑1 (MQ=255)
                                       tCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCgg                                                                                                                              >  2:979297/1‑119 (MQ=255)
                                       tCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCgg                                                                                                                              <  1:979297/119‑1 (MQ=255)
                                        cGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCgg                                                                                                                              >  3:1136761/1‑118 (MQ=255)
                                        cGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCgg                                                                                                                              <  4:1136761/118‑1 (MQ=255)
                                            aCCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGcc                                                                                            <  2:210976/148‑1 (MQ=255)
                                                gCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGc                                                                                             >  3:1134681/1‑143 (MQ=255)
                                                gCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGc                                                                                             <  4:1134681/143‑1 (MQ=255)
                                                    cgcgATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGc                                                                                    >  4:545079/1‑148 (MQ=255)
                                                        atatTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCgg                                                                                                                              <  4:508914/102‑1 (MQ=255)
                                                        atatTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCgg                                                                                                                              >  3:508914/1‑102 (MQ=255)
                                                               aCGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCAGGCGTCgtg                                                                                                                                            <  4:1121569/81‑3 (MQ=255)
                                                               aCGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCAGGCGTCgtg                                                                                                                                            >  3:1121569/1‑79 (MQ=255)
                                                                cgcgGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTg                                                                                              >  2:665222/1‑126 (MQ=255)
                                                                cgcgGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTg                                                                                              <  1:665222/126‑1 (MQ=255)
                                                                                       cACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTgg                                                                                                <  4:511619/101‑1 (MQ=255)
                                                                                       cACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTgg                                                                                                >  3:511619/1‑101 (MQ=255)
                                                                                           gAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCAt                                             >  1:256005/1‑148 (MQ=255)
                                                                                             aTACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCa                                                                   <  3:713532/124‑1 (MQ=255)
                                                                                             aTACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCa                                                                   >  4:713532/1‑124 (MQ=255)
                                                                                                  cAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGt                                                                <  3:348915/122‑1 (MQ=255)
                                                                                                  cAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGt                                                                >  4:348915/1‑122 (MQ=255)
                                                                                                  cAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTc                                      >  3:977545/1‑148 (MQ=255)
                                                                                                  cAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTc                                      <  3:545079/148‑1 (MQ=255)
                                                                                                    gggCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTc                                                                                                        <  4:849061/80‑1 (MQ=255)
                                                                                                    gggCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTc                                                                                                        >  3:849061/1‑80 (MQ=255)
                                                                                                              gCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGaa                          >  3:1015345/1‑148 (MQ=255)
                                                                                                                 gttgttATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCt                       >  1:149398/1‑148 (MQ=255)
                                                                                                                  ttgttATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCtt                      <  2:149398/148‑1 (MQ=255)
                                                                                                                  ttgttATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCtt                      <  4:1015345/148‑1 (MQ=255)
                                                                                                                        tGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCTTCGCCtt                >  1:554841/1‑148 (MQ=255)
                                                                                                                             ccACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAgg                                                                 >  4:811383/1‑94 (MQ=255)
                                                                                                                             ccACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAgg                                                                 <  3:811383/94‑1 (MQ=255)
                                                                                                                              cACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCTTCGCCTTCCGCGt          <  2:554841/148‑1 (MQ=255)
                                                                                                                                      cTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCTTCGCCTTCCGCGTCAGAAAGc  <  4:977545/148‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                                   
CCGGCAAAAATCGCTGCCATCACCAGATATGGGTTGGCATCGGCACCCGCCACGCGATATTCCACGCGGTGATTATGACGGTCGCCGCACGGAATACGCAGGGCGACGGTGCGGTTGTTATGGCCCCACGACGCCTGCGTCGGCACATACATTCCCGGCTGGAAGCGGCGATACGAGTTCACGTTTGGTGCCAGCAACGCCATCGACGACGGCATCAGGTCAATCATCCCGGCGAGCATCTTTTTCAGCAGCGGCGAATCTTCGCCTTCCGCGTCAGAAAGC  >  NC_000913/1359743‑1360024

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