Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,366,854 T→G 9.7% intergenic (‑81/+37) dcuA ← / ← aspA C4‑dicarboxylate transporter DcuA/aspartate ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,366,8540TG 9.7% 112.1 / 4.4 55intergenic (‑81/+37)dcuA/aspAC4‑dicarboxylate transporter DcuA/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base T (24/24);  new base G (3/3);  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 = 4.57e-01

CTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCG  >  NC_000913/4366721‑4366995
                                                                                                                                     |                                                                                                                                             
cTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATTAGTTTTACAATTAGCATATTAGCCTTCCTTGTTTTTAAACAAGATGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGc                                                                                                                                      <  1:267876/143‑1 (MQ=255)
 tGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTAACAAGAAAAAAGGCACGTCATCTGACGTGcc                                                                                                                                     <  1:16277/143‑1 (MQ=255)
   gcgcCCAATAAGATCGCCAGGGAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCCtt                                                                                                                                   <  2:447597/143‑1 (MQ=255)
     ttccAAGAACATCGCCAGCAAAACTAGGGTTAGTTCTTCACCTCTCCGATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCtttt                                                                                                                                 <  2:337479/141‑1 (MQ=255)
        cAAGAAGATCGCCAGCAAAACTATGATTAGTTCTACAACTAGCATATTATTTTTTCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTa                                                                                                                              <  1:65338/143‑1 (MQ=255)
           gaagaTCTCGAGCAAATCTATGATGATTTCTATAACTAGCATATTATCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTGACAACAAACCACACTCTTTATCTGACGTGCCTTTTTTCttt                                                                                                                           <  2:427244/143‑1 (MQ=255)
                   ccAGGAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTTTTTGTTTTTTAACAAGTTGATATTAGGTTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCAGGTCATCTGACGTGCCTTTTTTATTTGTACTAcc                                                                                                                   <  1:466281/143‑1 (MQ=255)
                        aaaaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTAGTACACTGTa                                                                                                              >  2:106683/1‑143 (MQ=255)
                          aaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACg                                                                                                            >  2:28753/1‑143 (MQ=255)
                           aCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGCCGTGCCTTTTTTATTTGTACTACCCTGTAcca                                                                                                           >  2:103021/1‑141 (MQ=255)
                            cTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTAGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCGGACGTGCCTTTGTGACTTGTACTACCCTGTACGAt                                                                                                          >  2:353422/1‑143 (MQ=255)
                             tATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGAtt                                                                                                         <  2:25131/143‑1 (MQ=255)
                                gatgaGTTCGACAACGATCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATAGTCGATTGTCATTGTCAAGTTACTGCTCACAAGAAAAAAGGCAAGTCATCTGACGTGCCTGTTTTATTTGTACTACCCTGTACGATTACt                                                                                                      <  2:203516/143‑1 (MQ=255)
                                      ttCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGCTCTTTGATTGTTATTTTTATGTTACTGCTCACAAGAAAAAAGACACGTCATCTAAAGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGc                                                                                                >  1:151256/1‑143 (MQ=255)
                                      ttCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGc                                                                                                <  2:463281/143‑1 (MQ=255)
                                              cTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTTTTTTTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGGCATCTGACGGGCCTTGTTGATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatca                                                                                        <  1:389774/143‑1 (MQ=255)
                                               tAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTCCGATTACTGTTCGCTTCCATCAg                                                                                       >  2:460413/1‑143 (MQ=255)
                                                   atatTAGCCTTGTTTGTTTTTTAACTAGTTTTTTTTTGTTTGTTTTTTTTAAGGTACTGCGCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGtata                                                                                   <  1:441734/143‑1 (MQ=255)
                                                         gCCTTCCTTGTTTTTTAATTTTTTTATAATAGATTGTTGTTTTTAAGTTAATGCTTACCAGACAACAGGCACGTCAGCTGTCGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTTTTCGCGTTCATCACTCTATCGttc                                                                             >  2:59176/1‑142 (MQ=255)
                                                           cttccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTTCTGCTCACAAGAAAAAAGGCACGTCATCTGAAGTGCCTTTTTTAGTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTg                                                                           <  1:249690/143‑1 (MQ=255)
                                                               cttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCACCTGACGTGCCTTTTCTATTTGTAATACCCTGTACGAATACAATTCGCATTCATCAGTATAGCGTTTTGCttt                                                                       >  1:263271/1‑143 (MQ=255)
                                                               cgtGTTTTGTAACAAGTTGATATTTGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCttt                                                                       <  1:117112/141‑1 (MQ=255)
                                                                ttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGTCGTGCCTTTTTTATTTGTACTACACTGTACAATTACTGTTCGCTTTCATCAGTATAGAGTTTTACTTTg                                                                      >  1:77626/1‑143 (MQ=255)
                                                                 tGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGAACGTCATCTGACGTGCCTTATATAATTGTAATACCCTGTACGATTACTGTTCGCTTTCTTCAGTATACAGTTTTAGTTTGt                                                                     >  1:371791/1‑143 (MQ=255)
                                                                  gTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGCCGTGCCTTTTTTATTTGTACTAGCCTGTACGATGACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTa                                                                    >  2:220264/1‑143 (MQ=255)
                                                                       ttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTTCTTTGTAAACCg                                                               >  1:291076/1‑143 (MQ=255)
                                                                        tAACAAGTTGATATTAGATTGTTATTTTTGAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACACTGTAAGATTACTGTTCACTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCgg                                                              >  1:68937/1‑143 (MQ=255)
                                                                           cAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTAACGTGCCTTTTTAATATGTACTACCCTGTACGATTACTGTTCGCTTTCATCACTATAGCGTTTTGCTTTGTAATCCGGGTg                                                           >  2:229741/1‑143 (MQ=255)
                                                                            agcttGATATTAGATTGTTATTTTTCAGTTACTGCTTACAATATAATCGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCTTGTACGATTACTGTTCGCTTTCATCTGTATAGCTTTTTGCTTTGTAAGCTGGGTGc                                                          <  2:177123/140‑1 (MQ=255)
                                                                                tGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAACAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCa                                                      <  2:249520/143‑1 (MQ=255)
                                                                                   tatTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCAGGTGCATCAGAt                                                   >  1:415055/1‑143 (MQ=255)
                                                                                    atTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTTCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCAGGGTGCCTCAGAtt                                                  >  1:352931/1‑143 (MQ=255)
                                                                                         aTTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTa                                             <  2:292138/143‑1 (MQ=255)
                                                                                         aTTGTGTTTTTTACGTGACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTa                                             <  2:112853/143‑1 (MQ=255)
                                                                                          tttttATTTTTATGTTTCTGCTCACAAGAAACAATGCACGTCATCTGTCGTGTCTTTTTTATTTGTACTACCCTGTACGATTACGGTTCGCTTTCAGCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTAc                                            <  1:268467/140‑1 (MQ=255)
                                                                                          ttGTTATTTTTAAGTTACTGCTCTCTTGAATTTATGCACGTCATCTTTCTTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTAc                                            <  1:134691/143‑1 (MQ=255)
                                                                                             ttATTTTTAAGTTACTGCTCACTTGAAAAAATTCATGTTTTCTGACATACCTTATTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGga                                         <  1:59217/143‑1 (MQ=255)
                                                                                              tATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTAATCAGTATAGCGTTTTACTTTGTAAGCCAGGTGCATCAGATTCTGTACGgag                                        >  1:337960/1‑143 (MQ=255)
                                                                                                    tAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTTCATCAGATTCTGTACGGAGAAAatc                                  >  1:68261/1‑142 (MQ=255)
                                                                                                          aCTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCCTCAGATCCTGTACGGAGCACATACCGTCa                            >  2:485251/1‑143 (MQ=255)
                                                                                                          aCTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCa                            >  2:204068/1‑143 (MQ=255)
                                                                                                              cTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATTAGTTTGGCGCAGACCAAAGTAAGCCGGATGCATCAGATTCTGTACGGAGAAAATATCGTAAAGtt                        >  1:349906/1‑143 (MQ=255)
                                                                                                               tCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTc                       >  2:347624/1‑143 (MQ=255)
                                                                                                                 actaggaaaaaGGCACGTCATGTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCAGCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCg                     <  2:184866/137‑1 (MQ=255)
                                                                                                                 acaAGAAAAAAGGCACGTCATCTGACGTGCCGTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCg                     <  2:313020/143‑1 (MQ=255)
                                                                                                                  caAGAACAGATGCACTTCATCTGACGTGCCTTTTTTATTTGTTCTACCCTGTTCTATTACTGTTCGCTTTCATCTGTATAGCGTTTTGTTTTGTAAGCTGTGGGTATCATTTTCTGTACGGAGAAAATATCCCCAAATTCCGc                    <  2:288758/143‑1 (MQ=255)
                                                                                                                   aaGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTATCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCt                   >  1:260951/1‑143 (MQ=255)
                                                                                                                    ataaaaaaGGCACGTCATCTGACGTGTCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTCCCGCAGTATATCGTTTTGCTCTACAAGCCGAGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCtt                  <  1:352360/141‑1 (MQ=255)
                                                                                                                     gagagagggCACGGCAGGTGTCGGGCCTTTTTTATTTGTACTTCCCGGTACGATTACTGTGCGCTTTCATCAGTGTAGCGTTTTGCTTTGTAAGCCGGGTGCAGCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTc                 <  2:125566/136‑1 (MQ=255)
                                                                                                                     gAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACACTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCACCAGCCTCTGTACGGAGAAAATATCGTCAAGGTCCGCTTc                 >  2:440197/1‑143 (MQ=255)
                                                                                                                           aGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGGCAAGTTCCGCGTGAGTCaa           >  2:168012/1‑143 (MQ=255)
                                                                                                                           aGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCCTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGTAAATATCGTCCAGTTCCGCTTCAGTCaa           >  1:232728/1‑143 (MQ=255)
                                                                                                                             gcgcggcATGTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACa         <  2:371612/137‑1 (MQ=255)
                                                                                                                                cGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCCTCAGATTCTGTACGGAGGAAATATCGTCAAGTTCCACTTCAGTCGACAGAc      >  2:283321/1‑143 (MQ=255)
                                                                                                                                 ttcttCTGACTTGCCTTTTTTATTTGTACTATCCTGTACGATTACTGTTCGCTTTCAGCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc     <  1:40490/139‑1 (MQ=255)
                                                                                                                                    aTCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCGAGTTGCGCTTCAGTCAACAGACcgcg  >  2:209745/1‑143 (MQ=255)
                                                                                                                                     |                                                                                                                                             
CTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCG  >  NC_000913/4366721‑4366995

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑