Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,632,029 A→T 9.1% A169A (GCT→GCA guaA ← GMP synthetase (glutamine aminotransferase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,632,0290AT 9.1% 183.7 / 8.9 55A169A (GCT→GCAguaAGMP synthetase (glutamine aminotransferase)
Reads supporting (aligned to +/- strand):  ref base A (25/25);  new base T (3/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

CGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAAC  >  NC_000913/2631903‑2632172
                                                                                                                              |                                                                                                                                               
cGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTa                                                                                                                                                 >  1:344407/1‑127 (MQ=255)
cGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTa                                                                                                                                                 <  2:344407/127‑1 (MQ=255)
        tCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAAc                                                                                                                                  >  2:63062/1‑134 (MQ=255)
        tCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAAc                                                                                                                                  <  1:63062/134‑1 (MQ=255)
            caGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATg                                                                                                                                        <  3:617425/124‑1 (MQ=255)
            caGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATg                                                                                                                                        >  4:617425/1‑124 (MQ=255)
            caGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTg                                                                                                                >  3:131253/1‑148 (MQ=255)
               ggCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATa                                                                                                                                           >  1:44362/1‑118 (MQ=255)
               ggCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATa                                                                                                                                           <  2:44362/118‑1 (MQ=255)
                  ttCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACg                                                                                                          <  4:131253/148‑1 (MQ=255)
                      cacTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTTGCCATTATgg                                                                                                                                       <  4:894952/115‑1 (MQ=255)
                      cacTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTTGCCATTATgg                                                                                                                                       >  3:894952/1‑115 (MQ=255)
                       acTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCt                                                                                                                 >  1:373443/1‑136 (MQ=255)
                       acTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCt                                                                                                                 <  2:373443/136‑1 (MQ=255)
                            cAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGt                                                                                                >  3:559957/1‑148 (MQ=255)
                                       cacaAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCg                                                                                                                      <  3:83933/115‑1 (MQ=255)
                                       cacaAAACGCTCCAGCACGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCCAACGGGCAGATTTCg                                                                                                                      >  4:83933/1‑115 (MQ=255)
                                           aaaCGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGtg                                                                                                       <  2:255588/126‑1 (MQ=255)
                                           aaaCGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGtg                                                                                                       >  1:255588/1‑126 (MQ=255)
                                             aCGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTa                                                                               <  4:559957/148‑1 (MQ=255)
                                                                cccgccagggTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCa                                                            >  4:583102/8‑148 (MQ=255)
                                                                   tGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACg                                                                                                          <  2:161497/99‑1 (MQ=255)
                                                                   tGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACg                                                                                                          >  1:161497/1‑99 (MQ=255)
                                                                       gggTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAg                                                                                                                                                <  1:819144/57‑1 (MQ=255)
                                                                       gggTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAg                                                                                                                                                >  2:819144/1‑57 (MQ=255)
                                                                         gTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACAt                                                   <  3:583102/148‑1 (MQ=255)
                                                                          tATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTaa                                                                              <  2:876085/120‑1 (MQ=255)
                                                                          tATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTaa                                                                              >  1:876085/1‑120 (MQ=255)
                                                                           aTGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTACCTTTATCGCCGTGGCTCATc                                                          <  4:370476/139‑1 (MQ=255)
                                                                           aTGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATc                                                          >  3:370476/1‑139 (MQ=255)
                                                                              aGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAAt                                                                                      <  2:113648/108‑1 (MQ=255)
                                                                              aGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAAt                                                                                      >  1:113648/1‑108 (MQ=255)
                                                                                 cACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCt                                                                                                                 <  4:1056039/78‑1 (MQ=255)
                                                                                 cACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCt                                                                                                                 >  3:1056039/1‑78 (MQ=255)
                                                                                 cACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGc                                                                                   <  2:548489/108‑1 (MQ=255)
                                                                                 cACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGc                                                                                   >  1:548489/1‑108 (MQ=255)
                                                                                 cACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCgg                                                                                 <  3:929156/110‑1 (MQ=255)
                                                                                 cACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCgg                                                                                 >  4:929156/1‑110 (MQ=255)
                                                                                               aCTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGt                                                                                                <  3:318718/81‑1 (MQ=255)
                                                                                               aCTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGt                                                                                                >  4:318718/1‑81 (MQ=255)
                                                                                                           aGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATc                 >  2:181953/1‑148 (MQ=255)
                                                                                                            gAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCt                <  1:181953/148‑1 (MQ=255)
                                                                                                                 gTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGt                          <  4:560591/133‑1 (MQ=255)
                                                                                                                 gTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGt                          >  3:560591/1‑133 (MQ=255)
                                                                                                                   ttttCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGa                                                <  4:137455/109‑1 (MQ=255)
                                                                                                                   ttttCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGa                                                >  3:137455/1‑109 (MQ=255)
                                                                                                                         tcGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTgg                                                                >  3:1054521/1‑87 (MQ=255)
                                                                                                                         tcGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTgg                                                                <  4:1054521/87‑1 (MQ=255)
                                                                                                                         tcGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTg                                                                 <  2:1082947/86‑1 (MQ=255)
                                                                                                                         tcGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTg                                                                 >  1:1082947/1‑86 (MQ=255)
                                                                                                                         ccGTTTGCCATTATGGCTAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTAc                                    >  4:1109661/2‑115 (MQ=255)
                                                                                                                         ccGTTTGCCATTATGGCTAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTAc                                    <  3:1109661/114‑1 (MQ=255)
                                                                                                                          cGTTTGCCATTATGGCTAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAAc  >  1:809472/1‑148 (MQ=255)
                                                                                                                             tAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCtt               >  2:600566/1‑132 (MQ=255)
                                                                                                                             tAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCtt               <  1:600566/132‑1 (MQ=255)
                                                                                                                             tAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACcg       <  3:1031845/140‑1 (MQ=255)
                                                                                                                             tAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACcg       >  4:1031845/1‑140 (MQ=255)
                                                                                                                              |                                                                                                                                               
CGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAAC  >  NC_000913/2631903‑2632172

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