Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB143 4,437,644 C→T D125N (GAC→AAC)  nicD ← N‑formylmaleamate deformylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,437,6440CT97.2% 114.1 / ‑4.0 36D125N (GAC→AAC) nicDN‑formylmaleamate deformylase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (18/17);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

TCGCTCCAGGTGGCACAGAAGGCGCGCATGTCATCGCCGCTCATGCCAACCGTGGCTTGGCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTCGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTGCGCGGATGGCAAAACGCGCCCCCATCGAGTGACCGAGCAGATGGTAGCTATCCAGGCCCAGCGCGGCGGCAAAGGCGGGAATATCCG  >  ACB143/4437514‑4437775
                                                                                                                                  |                                                                                                                                   
tcgCTCCAGGTGGCACAGAAGGCGCGCAAGTCATCGCCGCTCATGCCAACCGGGGCTTGGCGGCGCGAATCGACGTACCACGGCAACTTGCTCGGGTCGGCGCGGCGGCCCGGCCCGGCCACCGGGGGGTTGac                                                                                                                                  >  3:15801/1‑134 (MQ=255)
 cgCTCCAGGTGGCACAGAAGGCGCGCATGTCATCGCCGCTCATGCCAACCGTGGCTTGGCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGaccag                                                                                                                               >  2:271957/1‑136 (MQ=255)
      cAGGTGGCACAGAAGGCGCGCATGTCATCGCCGCTCATGCCAACCGTGGCTTGGCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGaccagcac                                                                                                                            >  2:2421/1‑134 (MQ=255)
       aGGTGGCACAGAAGGCGCGCATGTCATCGCCGCTCATGCCAACCGTGGCTTGGCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGaccagcaccag                                                                                                                         >  4:171459/1‑136 (MQ=255)
                aGAAGGCGCGCATGTCATCGCCGCTCATGCCAACCGTGGCTTGGCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAg                                                                                                                >  4:363030/1‑136 (MQ=255)
                   aGGCGCGCATGTCATCGCCGCTCATGCCAACCGAGGCTTGGCGGATCGAATCGACGTACCACGGCAACTTGATCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTcc                                                                                                             >  5:347336/1‑136 (MQ=255)
                    ggCGCGCATGTCATCGCCGCTCATGCCAACCGTGGCTTGGCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTccc                                                                                                            <  8:5209/136‑1 (MQ=255)
                    ggCGCGCATGTCATCGCCGCTCATGCCAACCGTGGCTTGGCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTccc                                                                                                            <  6:364821/136‑1 (MQ=255)
                       gcgcATGTCATCGCCGCTCATGCCAACCGTGGCTTGGCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGt                                                                                                         >  5:110249/1‑136 (MQ=255)
                              tCATCGCCGCTCATGCCAACCGTGGCTTGGCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCtt                                                                                                  <  6:363281/136‑1 (MQ=255)
                                 tCGCCGCTCATGCCAACCGTGGCTTGGCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTggg                                                                                               <  2:353765/136‑1 (MQ=255)
                                          aTGCCAACCGTGGCTTGGCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGCCACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTgcgc                                                                                      >  2:254237/1‑136 (MQ=255)
                                          aTGCCAACCGTGGCTTGGCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTgcgc                                                                                      >  1:133429/1‑136 (MQ=255)
                                          aTGCCAACCGTGGCTTGGCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTgcgc                                                                                      >  5:365008/1‑136 (MQ=255)
                                          aTGCCAACCGTGGCTTGGCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTgcgc                                                                                      >  5:89693/1‑136 (MQ=255)
                                          aTGCCAACCGTGGCTTGGCGGATCGAAACGCCGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGACCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGAGCCCCTTTGGCGGCTgcgc                                                                                      >  8:321204/1‑136 (MQ=255)
                                                      gCTTGGCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTGCGCGGATGGCAAAAc                                                                          <  3:205899/136‑1 (MQ=255)
                                                      gCTTGGCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTGCGCGGATGGCAAAAc                                                                          <  8:21198/136‑1 (MQ=255)
                                                        ttGGCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTGCGCGGATGGCAAAAcgc                                                                        <  5:76766/136‑1 (MQ=255)
                                                          ggCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTGCGCGGATGGCAAAAcgc                                                                        <  1:92998/134‑1 (MQ=255)
                                                                 cgaatcgaGGTACCACGGCAACTTGGTCGGGTAGGCGCGGCGGCCCGGCCCGGACACGGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTGGGGGGATGGCAAAACGCGCCCCCATc                                                               <  1:112292/136‑1 (MQ=255)
                                                                   aatcgaCGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTGCGCGGATGGCAAAACGCGCCCCCATc                                                               >  5:364402/1‑134 (MQ=255)
                                                                         cGTACCACGGCAACTTGCTCGGGTCGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTGCGCGGATGGCAAAACGCGCCCCCATCGAGTGAcc                                                       <  4:29017/136‑1 (MQ=255)
                                                                         cGTACCACGGCAACTTGCTCGGGTAGGCGCGGGGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTGCGCGGATGGCAAAACGCGCCCCCATCGAGTGAcc                                                       <  7:255436/136‑1 (MQ=255)
                                                                         cGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTGCGCGGATGGCAAAACGCGCCCCCATCGAGTGAcc                                                       <  4:370365/136‑1 (MQ=255)
                                                                           aaCCACGGCACCTTGCGCGGGAAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTCGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTGCGCGGATGGCAAAACGCGCCCCCATCGAGTGACCGa                                                     <  4:54087/135‑1 (MQ=255)
                                                                                    aaCTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGGGGCTGCGCGGATGGCAAAACGCGCCCCCATCGAGTGACCGAGCAGATGGt                                            <  6:312257/136‑1 (MQ=255)
                                                                                           tCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTGCGCGGATGGCAAAACGCGCCCCCATCGAGTGACCGAGCAGATGGTAGCTATc                                     >  2:23153/1‑136 (MQ=255)
                                                                                              ggTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTGCGCGGATGGCAAAACGCGCCCCCATCGAGTGACCGAGCAGATGGTAGCTATCCAg                                  <  8:67523/136‑1 (MQ=255)
                                                                                                           gcccggcccgGACACCGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTGCGCGGATGGCAAAACGCGCCCCCATCGAGTGACCGAGCAGATGGTAGCTATCCAGGCCCAGCgcggcg                     >  7:214872/1‑136 (MQ=255)
                                                                                                               ggcccggCCACCGGGGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGGTGCGCGGATGGCAACACGCGCCCCCATCGAGTGACCGAGCAGATGGTAGCTATCCAGGCCCAGCGCGGCGGCaa                 <  6:25191/136‑1 (MQ=255)
                                                                                                                         ccGGTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTGCGCGGATGGCAAAACGCGCCCCCATCGAGTGACCGAGCAGATGGTAGCTATCCAGGCCCAGCGCGGCGGCAAAGGCGGGAat       <  4:198044/136‑1 (MQ=255)
                                                                                                                           ggTGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGGGCCCCTTGGGCGGCTGCGCGGATGGCAAAACGCGCCCCCATCGAGTGACCGAGCAGATGGTAGCTATCCAGGCCCAGCGCGGCGGCAAAGGCGGGAatat     <  1:365640/136‑1 (MQ=255)
                                                                                                                             tGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTTCCACTTGGGCGGCTGCGCGGATGCCAAAACGCGCCCCCATCGAGTGACCGAGCAGCTGGTAGCTATACAGGCACAGCGCGGCGGCAAAGGCGGGAATATcc   >  7:332326/1‑136 (MQ=255)
                                                                                                                             tGGGTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTGCGCGGATGGCAAAACGCGCCCCCATCGAGTGACCGAGCAGATGGTAGCTATCCAGGCCCAGCGCGGCGGCAAAGGCGGGAATATcc   >  2:39602/1‑136 (MQ=255)
                                                                                                                              gggTTGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTGCGCGGATGGCAAAACGCGCCCCCATCGAGTGACCGAGCAGATGGTAGCTATCCAGGCCCAGCGCGGCGGCAAAGGCGGGAATATCCg  >  1:18778/1‑136 (MQ=255)
                                                                                                                                  |                                                                                                                                   
TCGCTCCAGGTGGCACAGAAGGCGCGCATGTCATCGCCGCTCATGCCAACCGTGGCTTGGCGGATCGAATCGACGTACCACGGCAACTTGCTCGGGTAGGCGCGGCGGCCCGGCCCGGACACCGGTGGGTCGACCAGCACCAGCCGCTGCAGTCCCGGTGCCCCTTGGGCGGCTGCGCGGATGGCAAAACGCGCCCCCATCGAGTGACCGAGCAGATGGTAGCTATCCAGGCCCAGCGCGGCGGCAAAGGCGGGAATATCCG  >  ACB143/4437514‑4437775

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: