Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,028,434 G→A 14.8% D23N (GAC→AAC)  yccU → predicted CoA‑binding protein with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,028,4340GA14.8% 68.8 / 4.7 27D23N (GAC→AAC) yccUpredicted CoA‑binding protein with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base G (9/14);  new base A (2/2);  total (11/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

AATAATGAAAGAAACCGATATTGCAGGCATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCAGGGGTATCACGTCATTCCG  >  AP009048/1028364‑1028502
                                                                      |                                                                    
aataatGAAAGAAACCGATATTGCAGGCATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCg                                                                      <  1:2480497/71‑1 (MQ=255)
  taatGAAAGAAACCGATATTGCAGGCATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGAc                                                                    <  1:4603151/71‑1 (MQ=255)
  taatGAAAGAAACCGATATTGCAGGCATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGAc                                                                    <  1:2233743/71‑1 (MQ=255)
    atGAAAGAAACCGATATTGCAGGCATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGACa                                                                   >  1:3243702/1‑70 (MQ=255)
                  tatTGCAGGCATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTcc                                                      >  1:1061041/1‑69 (MQ=255)
                  tatTGCAGGCATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTcc                                                      >  1:55711/1‑69 (MQ=255)
                          gCATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATcg                                             <  1:149767/70‑1 (MQ=255)
                           cATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATcgcg                                           <  1:633768/71‑1 (MQ=255)
                           cATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATcgcg                                           >  1:1913111/1‑71 (MQ=255)
                            aTTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGt                                          <  1:768623/71‑1 (MQ=255)
                                    gTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAat                                  <  1:3849536/71‑1 (MQ=255)
                                     tCCACTCAAACCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGATATCGCGTGATGaaa                                   >  1:2273415/1‑69 (MQ=255)
                                           cacaCCATTGCGCTGGTGGGGGCAAGCAACAAACCCGAt                                                           >  1:54487/1‑39 (MQ=255)
                                            acacCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTAt                                               >  1:384943/1‑50 (MQ=255)
                                            acacCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTc                          >  1:951001/1‑71 (MQ=255)
                                            acacCATTGCGCTGGTGGGGGCAAGCAACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTc                          >  1:4099232/1‑71 (MQ=255)
                                             cacCATTGCGCTGGTGGGGGCAAGCGACAAACCCg                                                             <  1:1815953/35‑1 (MQ=255)
                                                cATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGAc                       <  1:1702420/70‑1 (MQ=255)
                                                 aTTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCa                     <  1:2707330/71‑1 (MQ=255)
                                                 aTTGCGCTGGTGGGGGCAAGCAACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGAc                       <  1:4031094/69‑1 (MQ=255)
                                                          gtggGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCAGGGGTATc             <  1:4477549/70‑1 (MQ=255)
                                                          gtggGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCAGGGGTATc             <  1:915838/70‑1 (MQ=255)
                                                          gtggGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTAATGAAATATCTACTCGACCAGGGGTATCa            >  1:1173333/1‑71 (MQ=255)
                                                           tggGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCAGGGGTATCa            <  1:166402/70‑1 (MQ=255)
                                                           tggGGGCAAGCAACAAACCCGATCGTCCCAGCTATCGCGtgatg                                      <  1:3738402/44‑1 (MQ=255)
                                                                gCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGAcc                      <  1:1469616/55‑1 (MQ=255)
                                                                 cAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCAGGGGTATCACGTCAt      <  1:3324385/70‑1 (MQ=255)
                                                                    gCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCAGGGGTATCACGTCATTCCg  >  1:508739/1‑71 (MQ=255)
                                                                      |                                                                    
AATAATGAAAGAAACCGATATTGCAGGCATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCAGGGGTATCACGTCATTCCG  >  AP009048/1028364‑1028502

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