Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,685,094 A→C 19.4% Q451H (CAA→CAC)  argD_4 → Acetylornithine/succinyldiaminopimelateaminotrans ferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,685,0940AC19.4% 96.3 / 4.4 50Q451H (CAA→CAC) argD_4Acetylornithine/succinyldiaminopimelateaminotrans ferase
Reads supporting (aligned to +/- strand):  ref base A (26/13);  new base C (2/9);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 6.11e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTCGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCGTGGAAAAAGGCCAGG  >  ACB122/4684960‑4685224
                                                                                                                                      |                                                                                                                                  
aCTTCGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACaa                                                                                                                                    >  2:130957/1‑135 (MQ=255)
    cGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGCGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTCCAAgcgc                                                                                                                                >  3:104598/1‑135 (MQ=255)
     gCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAgcgc                                                                                                                                <  1:271665/134‑1 (MQ=255)
        aGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAACCCCACATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCACCGGCAGCCCGGTCCAAGCgccg                                                                                                                              >  8:81326/1‑133 (MQ=255)
             ccccGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTGTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGCTACACGCgccggcagccg                                                                                                                       >  3:57973/1‑135 (MQ=255)
             ccccGGCCTGCAGCCGGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTATCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCgccggcagccg                                                                                                                       <  8:142367/135‑1 (MQ=255)
               ccGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCgccggcagccggc                                                                                                                     >  8:134437/1‑135 (MQ=255)
               ccGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCgccggcagccggc                                                                                                                     >  2:154464/1‑135 (MQ=255)
                       cagcctgcaCGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGgca                                                                                                             <  4:192581/135‑1 (MQ=255)
                         gcctgcaCGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGgcagc                                                                                                           <  1:195170/135‑1 (MQ=255)
                               aCGGGTAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAgcg                                                                                                     >  4:48763/1‑135 (MQ=255)
                               aCGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCACCGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAgcg                                                                                                     >  5:280112/1‑135 (MQ=255)
                               aCGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACACGCGCCGGCAGCCGGCCCCTGGCAGCACAgcg                                                                                                     >  1:251171/1‑135 (MQ=255)
                               aCGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAgcg                                                                                                     >  1:51232/1‑135 (MQ=255)
                                        aTCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGcc                                                                                            <  4:117575/135‑1 (MQ=255)
                                          cGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCgg                                                                                          >  3:174223/1‑135 (MQ=255)
                                          cGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCgg                                                                                          >  1:42084/1‑135 (MQ=255)
                                            ccTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCgggg                                                                                        >  4:170636/1‑135 (MQ=255)
                                               gTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGctg                                                                                     >  3:190370/1‑135 (MQ=255)
                                                  ccAAACCCATACCGCGCGCCCCGCGGCGCCCGCCCCCGGCGCCCTTTTCGGCGCGCTGCAGGCCCCCACCGGCAGCCCGGAACACGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGctgcct                                                                                  <  1:130957/134‑1 (MQ=255)
                                                    aaaCCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGctgcctgc                                                                                >  8:6719/1‑135 (MQ=255)
                                                    aaaCCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGctgcctgc                                                                                >  3:115023/1‑135 (MQ=255)
                                                    aaaCCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGctgcctgc                                                                                >  3:68878/1‑135 (MQ=255)
                                                      aCCCATATCGATCGCCCCGCGTAGCCCGCCACCTGCGCCCTTTTCGTCGCGCTGCATGTGCCCACCGGCAGCCCGGTACACGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTgcgc                                                                              <  7:134437/135‑1 (MQ=255)
                                                      aCCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCGCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCGGGCAGCCGGCACCTGGCAGCACTGCGGCGACGGCCGGGGCTGCCTgcgc                                                                              <  2:8445/135‑1 (MQ=255)
                                                      aCCCATATCGAGCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTgcgc                                                                              <  5:69220/135‑1 (MQ=255)
                                                        ccATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCAc                                                                            >  1:37808/1‑135 (MQ=255)
                                                              cgatcgCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCa                                                                      >  6:257309/1‑135 (MQ=255)
                                                              cgatcgCCCCGAGGAGCCCGCCACCTGCGCCCTGTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCa                                                                      <  6:250273/135‑1 (MQ=255)
                                                                   gCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTggg                                                                 >  4:41893/1‑135 (MQ=255)
                                                                            agCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGc                                                        <  4:45925/135‑1 (MQ=255)
                                                                                 gCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCaa                                                   >  3:113069/1‑135 (MQ=255)
                                                                                 gCCACCTGCGCCCTTTTCGTCGAGCTGCATCTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGCCGGCCGGGGAGGCGTCCCCCCCCCGCTTGGGGGTTGGTGGATGCaa                                                   >  2:178939/1‑135 (MQ=255)
                                                                                 gCCACCTGCGCCCGTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCCCCCCGCATTGGGCTTTGTGGCTGCaa                                                   >  7:47735/1‑135 (MQ=255)
                                                                                                  cgtcgAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACg                                      <  4:55464/131‑1 (MQ=255)
                                                                                                  cgtcgAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACg                                      >  3:55464/1‑131 (MQ=255)
                                                                                                   gtcgAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCg                                 >  5:285339/1‑135 (MQ=255)
                                                                                                   gtcgAGCTGCATGGGCCCAAGGGCTGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCg                                 <  4:139786/135‑1 (MQ=255)
                                                                                                   ggcgcgCTGCATGTGCCCCACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCg                                 <  5:117478/130‑1 (MQ=255)
                                                                                                                      ccGGCCGCCCGGTACCCGCGCCGGCAGCCGGCACGGGGAAGCACCGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCGTgg              <  3:21154/134‑1 (MQ=255)
                                                                                                                       cGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCGTGGa             <  8:64289/135‑1 (MQ=255)
                                                                                                                         gCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCTGCAGACGGCCAGGCCGTGGaa            >  4:92661/1‑134 (MQ=255)
                                                                                                                         gCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCGTGGaa            >  4:261384/1‑134 (MQ=255)
                                                                                                                            gCCCGGTCCACGCGCCGGCCGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCGTGGAAAAag        <  4:190370/135‑1 (MQ=255)
                                                                                                                            gCCCGGTACACGGGCCGGCTGCGGGCACCTGGCAGCACTGCGGCGACGGCCGGGGCGGCCTGCGCACCCCGCATTGGGGTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCGTGGAAAAag        <  3:48763/135‑1 (MQ=255)
                                                                                                                            gCCCGGTACACGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCGTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCGTGGAAAAag        <  3:88921/135‑1 (MQ=255)
                                                                                                                            gCCCGGTACACCCGCCGGCAGCCGGCACCGGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCGTGGAAAAag        <  3:191544/135‑1 (MQ=255)
                                                                                                                            gCCCGGCCCACGCGCCGGCCGCGGGCCCCGGGCCGCCCCGGGGCGAGGGCCGGGGCGGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCGGGAAGACGCACCGGCAGACGGCCAGGCCGTGGAAAAag        <  5:233575/135‑1 (MQ=255)
                                                                                                                            gCCCGGCCCACGCGCCGGCAGCCGGCACCTGGCAGCACCGGGGCGACGGCGGGGGCTGCCGGCGCACCCGGCTTTGGGCTTTGGGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCGTGGAAAAag        <  3:172513/135‑1 (MQ=255)
                                                                                                                                ggTACAATCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCGTGGAAAAaggaca    >  6:264930/1‑135 (MQ=255)
                                                                                                                                  tACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCGTGGAAAAaggccagg  >  8:259016/1‑135 (MQ=255)
                                                                                                                                      |                                                                                                                                  
ACTTCGCAAGCCGCCCCGGCCTGCAGCCTGCACGGGCAGTATCGCCTGTCGCAAACCCATATCGATCGCCCCGAGGAGCCCGCCACCTGCGCCCTTTTCGTCGAGCTGCATGTGCCCAACGGCAGCCCGGTACAAGCGCCGGCAGCCGGCACCTGGCAGCACAGCGGCGACGGCCGGGGCTGCCTGCGCACCCCGCATTGGGCTTTGTGGCTGCAAGGCCTGGAAGACGCACCGGCAGACGGCCAGGCCGTGGAAAAAGGCCAGG  >  ACB122/4684960‑4685224

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