Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 235,575 G→A 8.3% R420R (CGC→CGT)  argH ← argininosuccinate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947235,5750GA 8.3% 131.7 / 2.2 48R420R (CGC→CGT) argHargininosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base G (24/20);  new base A (2/2);  total (26/22)
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

GCCAACAGGGCCTTGCCACGAACGACGGCAGCACGCACCTGCGCCGGCGCCGTACCACCAATGTGGTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACGGCGAACCAGATAGTCAGCCAGGTCGGTAGCGGTGGAGAAACCGC  >  NC_002947/235438‑235723
                                                                                                                                         |                                                                                                                                                    
gCCAACAGGGCCTTGCCACGAACGACGGCAGCACGCACCTGCGCCGGCGCCGTACCACCAATGTGGTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTc                                                                                                                                           <  5:345214/149‑1 (MQ=255)
            ttGCCACGAACGACGGCAGCACGCACCTGCGCCGGCGCCGTACCACCAATGTGGTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCAtc                                                                                                                                 >  8:136044/1‑147 (MQ=255)
                 acgaatgaCGGCAGCACGCACCTGCGCCGGCGCCGTACCACCAATGTGGTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGACGCAGTTCGTCCAGGCTCATCTCGGc                                                                                                                            >  6:361487/1‑147 (MQ=255)
                 acgaacgaCGGCAGCACGCACCTGCGCCGGCGCCGTACCACCAATGTGGTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCAGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGc                                                                                                                            >  5:284347/1‑147 (MQ=255)
                                    acCTGCGCCGGCGCCGTACCACCAATGTGGTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGTCTCATCTCGGCCAGGTCCTTGCCAGTGTCca                                                                                                        >  6:254707/1‑148 (MQ=255)
                                    acCTGCGCCGGCGCCGTACCACCAATGTGGTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCAGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCcac                                                                                                       <  6:284347/149‑1 (MQ=255)
                                     ccTGCGCCGGCGCCGTACCACCAATGTGGTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCcaca                                                                                                      >  7:233252/1‑149 (MQ=255)
                                      cTGCGCCGGCGCCGTACCACCAATGTGGTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAg                                                                                                                         >  1:14559/1‑129 (MQ=255)
                                      cTGCGCCGGCGCCGTACCACCAATGTGGTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCacac                                                                                                     >  5:85260/1‑149 (MQ=255)
                                      cTGCGCCGGCGCCGTACCACCAATGTGGTCACGCGCATTCACCGACCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAg                                                                                                                         <  2:14559/129‑1 (MQ=255)
                                                 ccgTACCACCAATGTGGTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCAcc                                                                                          <  7:284214/149‑1 (MQ=255)
                                                  cgTACCACCAATGTGGTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACcg                                                                                         >  7:128165/1‑149 (MQ=255)
                                                   gTACCACCAATGTGGTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACcgc                                                                                        <  5:254707/149‑1 (MQ=255)
                                                    tACCACCAATGTGGTCACGCGGATTAACCGAGCCTACCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTGCAGGCTCATCTCGGGCAGGTCCTTGCCAGTGTGCAGGCCAGACTTCACggcg                                                                                       >  4:19578/1‑149 (MQ=255)
                                                       caccaATGTGGTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTgg                                                                                    >  8:119540/1‑149 (MQ=255)
                                                        accaATGTGGTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCAGGGCCACCACGTCGTGCTCGATCTGGTGGCTGAATTGGCGCAGTTCGTCGAGGCTCATCTCGCCCAGGTCGTTGCGAGTGTCCACACCATACATCACCGCGTGGc                                                                                   >  4:95536/1‑149 (MQ=255)
                                                        accaATGTGGTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGc                                                                                   >  1:330896/1‑149 (MQ=255)
                                                               tgGTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGGCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGAt                                                                            <  5:270443/149‑1 (MQ=255)
                                                                 gTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTgg                                                                                    <  7:220252/139‑1 (MQ=255)
                                                                 gTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTgg                                                                                    >  8:220252/1‑139 (MQ=255)
                                                                 gTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGAtct                                                                          >  2:136765/1‑149 (MQ=255)
                                                                 gTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGAtct                                                                          >  8:173721/1‑149 (MQ=255)
                                                                      gcgcATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGACGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTgg                                                                     >  3:239964/1‑149 (MQ=255)
                                                                              aCCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCa                                                                                                                                                  <  7:284043/64‑1 (MQ=255)
                                                                              aCCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCa                                                                                                                                                  >  8:284043/1‑64 (MQ=255)
                                                                                 gAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGaa                                                          <  7:119540/149‑1 (MQ=255)
                                                                                      ttCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACgg                                                       >  4:215889/1‑147 (MQ=255)
                                                                                        ccAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAg                                                    <  7:136044/148‑1 (MQ=255)
                                                                                         cAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGc                                                  >  7:245941/1‑149 (MQ=255)
                                                                                                 gCACGGCAAACACGTACTGCTCGATCTCGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGACCTTGCCAGTGTCCACACCATAGTTCTCCTCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACGGCg                                           <  7:54487/148‑1 (MQ=255)
                                                                                                     ggCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACg                                              >  8:329409/1‑141 (MQ=255)
                                                                                                     ggCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACg                                              <  7:329409/141‑1 (MQ=255)
                                                                                                             cGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACGGCGAACCAGATAGTCa                              <  8:233252/149‑1 (MQ=255)
                                                                                                              gTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACGGCGAACCAGATAGTCAg                             >  5:66409/1‑149 (MQ=255)
                                                                                                               tCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACGGCGAACCAGATAGTCAGc                            <  6:85260/149‑1 (MQ=255)
                                                                                                               tCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACGGCGAACCAGATAGTCAGc                            <  2:60322/149‑1 (MQ=255)
                                                                                                                     tCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAg                                                    <  1:231381/119‑1 (MQ=255)
                                                                                                                     tCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAg                                                    >  2:231381/1‑119 (MQ=255)
                                                                                                                       gATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACGGCGAACCAGATAGTCAGCCAGGTc                      <  8:245941/147‑1 (MQ=255)
                                                                                                                             gTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACGGCGAACCAGATAGTCAGCCAGGTCGGTAGCgg              >  7:47103/1‑149 (MQ=255)
                                                                                                                             gTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACGGCGAACCAGATACTCAGCCAGGACGCTAGCgg              >  4:308436/1‑149 (MQ=255)
                                                                                                                             gTCGCTGAATTGACGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACGGCGAACCAGATAGTCAGCCAGGTCGGTAGCg               <  4:239964/148‑1 (MQ=255)
                                                                                                                                 cTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAAAGATCTCGTGGCAGTCACGGATCGGCAGGCCACGGCGAACCAGATAGTCAGCCAGGTCGGTAGCGGTGga          <  3:19578/149‑1 (MQ=255)
                                                                                                                                 cTGAATTGACGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACGGCGAACCAGATAGTCAGCCAGGTCGGTAGCGGTGga          <  5:361487/149‑1 (MQ=255)
                                                                                                                                    aaTTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACGGCGAACCAGATAGTCAGCCAGGTCGGTAGCGGTGgaga        >  7:97932/1‑148 (MQ=255)
                                                                                                                                      ttGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACg                                                        <  7:189171/98‑1 (MQ=255)
                                                                                                                                      ttGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACg                                                        >  8:189171/1‑98 (MQ=255)
                                                                                                                                      ttGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACGGCGAACCAGATAGTCAGCCAGGTCGGTAGCGGTGGAGAAAc     <  7:82164/149‑1 (MQ=255)
                                                                                                                                         gcgcAGTTCGTCCAGGCTGATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACGGCGAACCAGTTAGTCAGCCAGGTCGGTAGCGGTGGAGAAACcgc  <  3:308436/149‑1 (MQ=255)
                                                                                                                                         gcgcAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACg                                              <  3:97094/105‑1 (MQ=255)
                                                                                                                                         gcgcAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACg                                              >  4:97094/1‑105 (MQ=255)
                                                                                                                                         gcgcAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACGGCGAACCAGa                                    >  5:152449/1‑115 (MQ=255)
                                                                                                                                         gcgcAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACGGCGAACCAGa                                    <  6:152449/115‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
GCCAACAGGGCCTTGCCACGAACGACGGCAGCACGCACCTGCGCCGGCGCCGTACCACCAATGTGGTCACGCGCATTCACCGAGCCTTCCAGCGTCAGCACGGCAAACACGTCCTGCTCGATCTGGTCGCTGAATTGGCGCAGTTCGTCCAGGCTCATCTCGGCCAGGTCCTTGCCAGTGTCCACACCATACTTCACCGCGTGGCCAACGATCTCGTGGCAGTCACGGAACGGCAGGCCACGGCGAACCAGATAGTCAGCCAGGTCGGTAGCGGTGGAGAAACCGC  >  NC_002947/235438‑235723

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑