Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,162,433 C→T 8.5% A421V (GCC→GTC)  argD_2 → Acetylornithine/succinyldiaminopimelateaminotrans ferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,162,4330CT 8.5% 146.8 / 2.2 47A421V (GCC→GTC) argD_2Acetylornithine/succinyldiaminopimelateaminotrans ferase
Reads supporting (aligned to +/- strand):  ref base C (20/23);  new base T (2/2);  total (22/25)
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

CGATCCGTCCTGCCATATCGCGCTGTCGGCATTGCGCCGGTTGTTCGCCGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCGCTGGTACAGGCGTCGGCCCAGTACAACGATGCCTTGGGTTACTTGATCGC  >  ACB122/3162308‑3162565
                                                                                                                             |                                                                                                                                    
cGATCCGTCCTGCCATATCGCGCTGTCGGCATTGCGCCGGTTGTTCGCCGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGTCTATCCGGa                                                                                                                             >  6:172563/1‑135 (MQ=255)
cGATCCGTCCTGCCATATCGCGCTGTCGGCATTGCGCCGGTTGTTCGCCGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGTCTATCCGGa                                                                                                                             >  6:165770/1‑135 (MQ=255)
cGATCCGTCCTGCCATATCGCGCTGTCGGCATTGCGCCGGTTGTTCGCCGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGa                                                                                                                             >  1:181919/1‑135 (MQ=255)
      gTCCTGCCATATCGCGCTGTCGGCATTGCGCCGGTTGTTCGCCGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGAt                                                                                                                       >  2:158845/1‑135 (MQ=255)
      gTCCTGCCATATCGCGCTGTCGGCATTGCGCCGGTTGTTCGCCGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGAt                                                                                                                       <  3:150229/135‑1 (MQ=255)
      gTCCTGCCATATCGCGCTGTCGGCATTGCGCCGGTTGTTCGCCGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGAt                                                                                                                       >  5:197587/1‑135 (MQ=255)
        ccTGCCATATCGCGCTGGCGGCATTGCGGCGGTTGTTCGCCGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCg                                                                                                                     <  1:62361/135‑1 (MQ=255)
               tatCGCGCTGTCGGCATTGCGCCGGTTGTTCGCCGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGAtgt                                                                                                              <  4:381756/135‑1 (MQ=255)
                    cgcTGTCGGCATTGCGCCGGTTGTTCGCCGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCgg                                                                                                        >  7:192009/1‑136 (MQ=255)
                     gcTGTCGGCATTGCGCCGGTTGTTCGCCGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCgg                                                                                                        >  4:142115/1‑135 (MQ=255)
                         tCGGCATTGCGCCGGTTGTTCGCCGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGt                                                                                                    >  1:299396/1‑135 (MQ=255)
                            gCATTGCGCCGGTTGTTCGCCGATGACGACGCGATGCTCAAAGCCATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCGGAGGTGCGCGGGCGTGGc                                                                                                 >  5:34958/1‑135 (MQ=255)
                              aTTGCGCCGGTTGTTCGCCGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCtt                                                                                               >  1:336987/1‑135 (MQ=255)
                                  cgcCGGTTGTTCGCCGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCt                                                                                            <  1:307430/134‑1 (MQ=255)
                                  cgcCGGTTGTTCGCCGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCt                                                                                            <  6:385608/134‑1 (MQ=255)
                                          gttCGCCGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTg                                                                                  >  8:346680/1‑136 (MQ=255)
                                               ccGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAAtt                                                                             <  7:18792/136‑1 (MQ=255)
                                               ccGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAAt                                                                              >  6:325833/1‑135 (MQ=255)
                                                cGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGTCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAAtt                                                                             <  5:165770/135‑1 (MQ=255)
                                                cGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGTCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAAtt                                                                             <  5:172563/135‑1 (MQ=255)
                                                cGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAAtt                                                                             <  5:302624/135‑1 (MQ=255)
                                                      cgacgCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACg                                                                        <  5:325833/134‑1 (MQ=255)
                                                                cTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAgtg                                                               >  4:238723/1‑133 (MQ=255)
                                                                    aaGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTcag                                                           <  1:158845/133‑1 (MQ=255)
                                                                       aCATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCt                                                      <  6:197587/135‑1 (MQ=255)
                                                                           aCGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCGCTg                                                 <  7:279346/136‑1 (MQ=255)
                                                                             gCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATCGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCTCGCCATGAGTGTCAGCTCCTCGCtc                                                 >  1:330216/1‑133 (MQ=255)
                                                                                   tgggCGCATACCTCAACGCCTCGATAAATTAAGTTAATAGGGCCTATCCGTTTGTGATCGCCGAGGTGCGGGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCGCTGGTACAgg                                          <  6:34958/134‑1 (MQ=255)
                                                                                   aGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTAGCCGGATGTGATCGCCGATGTGCGCGGGCGGGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCGCTGGTACAgg                                          <  2:181919/135‑1 (MQ=255)
                                                                                    gggCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCGCTGGTACAGGc                                         >  4:212634/1‑135 (MQ=255)
                                                                                      gcgcCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCGCTGGTACAGGCGt                                       <  2:207162/135‑1 (MQ=255)
                                                                                       cgcCTACCTCAAAGGCTCTCGCAATGAACTGTAGTTGGCCTATCCGGATGTTATAGCCGATGTGCGCGGGCGTGGCTTGCGGGTGGGGCTGGAATTGAACTACTTTAGTGGCAGCATCTCTCTGGTAGAGGCGTCg                                     <  7:260253/136‑1 (MQ=255)
                                                                                       cgcCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCTCTGGTACAGGCGTCg                                     <  8:101138/136‑1 (MQ=255)
                                                                                        gcCTCCCTCAAAGCCGCGCCCAATGAACTGAAGATGGCCTAGCCGGATGTGAGCGCCGATGGGCGCGGGCGGGGCTGGCTGGTGGGGCTGGAATGGCACGACTTGAGTGTCAGCAGCGCGCTGGTACAGgcggcg                                     <  5:20182/135‑1 (MQ=255)
                                                                                        gcCTACCTCAAAGCCTCGCTCCATGACCTGAAGATGGCCTATCCGGATGGGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCGCTGGTACAGGCGTCg                                     <  3:42789/135‑1 (MQ=255)
                                                                                        gcCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCGCTGGTACAGGCGTCg                                     <  3:60470/135‑1 (MQ=255)
                                                                                        gcCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCGCTGGTACAGGCGTCg                                     <  3:81598/135‑1 (MQ=255)
                                                                                            aCCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCGCTGGTACAGGCGTCGGccc                                 >  5:214283/1‑135 (MQ=255)
                                                                                            aCCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCGCTGGTACAGGCGTCGGc                                   >  3:377829/1‑133 (MQ=255)
                                                                                              cTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCGCTGGTACAGGCGTCGGCCCa                                >  3:27557/1‑134 (MQ=255)
                                                                                                  aaGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCGCTGGTACAGGCGTCGGCCCAGTACa                           <  6:214283/135‑1 (MQ=255)
                                                                                                  aaGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCGCTGGTACAGGCGTCGGCCCAGTACa                           <  4:27557/135‑1 (MQ=255)
                                                                                                           tCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGCATTGCACGACTTGAGTGTCAGCAGCTCGCTGGTACAGGCGTCGGCCCAGTACCACGCTGccg                  >  3:34595/1‑134 (MQ=255)
                                                                                                           tCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTTAGTGTCAGAAGCTCGCTGGTAAAGGCGTCGGCCCAGTAAAACGATGCCt                  >  6:348477/1‑135 (MQ=255)
                                                                                                           tCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCGCTGGTACAGGCGTCGGCCCAGTACAACGATGCCt                  >  2:17811/1‑135 (MQ=255)
                                                                                                                  aCTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCGCTGGTACAGGCGTCGGCCCAGTACAACGATGCCTTGGGTTAc          >  7:62525/1‑136 (MQ=255)
                                                                                                                          tGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCGTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCGCTGGTACAGGCGTCGGCCCAGTACAACGATGCCTTGGGTTACTTGATCGc  <  8:192009/136‑1 (MQ=255)
                                                                                                                             |                                                                                                                                    
CGATCCGTCCTGCCATATCGCGCTGTCGGCATTGCGCCGGTTGTTCGCCGATGACGACGCGATGCTCAAAGACATACGCACCAAGGGCGCCTACCTCAAAGCCTCGCTCAATGAACTGAAGATGGCCTATCCGGATGTGATCGCCGATGTGCGCGGGCGTGGCTTGCTGGTGGGGCTGGAATTGCACGACTTGAGTGTCAGCAGCTCGCTGGTACAGGCGTCGGCCCAGTACAACGATGCCTTGGGTTACTTGATCGC  >  ACB122/3162308‑3162565

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