Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,630,224 C→G 10.8% R442G (CGA→GGA)  pitA_2 → Low‑affinity inorganic phosphate transporter 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,630,2240CG10.8% 103.6 / 7.5 37R442G (CGA→GGA) pitA_2Low‑affinity inorganic phosphate transporter 1
Reads supporting (aligned to +/- strand):  ref base C (18/14);  new base G (2/2);  total (21/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 = 2.56e-01

ATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAGTGGTGGCGATGTGCACCATTGGCGCGGCGGATATGTTCGGGCTGCCGGTGTC  >  ACB122/4630099‑4630351
                                                                                                                             |                                                                                                                               
aTGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTAtt                                                                                                                         >  8:35820/1‑134 (MQ=255)
aTGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTAtt                                                                                                                         >  8:32790/1‑134 (MQ=255)
                            ggtCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCaa                                                                                            <  6:281449/133‑1 (MQ=255)
                             ttCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCa                                                                                             <  4:52336/133‑1 (MQ=255)
                             ttCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCa                                                                                             >  3:52336/1‑133 (MQ=255)
                                  ggCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCa                                                                                      <  5:276776/135‑1 (MQ=255)
                                  ggCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCa                                                                                      <  7:35820/135‑1 (MQ=255)
                                       aGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGGGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCa                                                                                 <  3:264777/135‑1 (MQ=255)
                                               gacgCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTGGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGc                                                                          >  4:132709/1‑134 (MQ=255)
                                                acgCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGGGACGTATTGTAGTCACCGTTGGCgag                                                                                                       >  4:60480/1‑104 (MQ=255)
                                                acgCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGGGACGTATTGTAGTCACCGTTGGCgag                                                                                                       >  4:56437/1‑104 (MQ=255)
                                                acgCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGGGACGTATTGTAGTCACCGTTGGCgag                                                                                                       <  3:60480/104‑1 (MQ=255)
                                                acgCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGGGACGTATTGTAGTCACCGTTGGCgag                                                                                                       <  3:56437/104‑1 (MQ=255)
                                                    cACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAgg                                                                    >  5:65891/1‑135 (MQ=255)
                                                           tataTTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCg                                                             <  8:11452/135‑1 (MQ=255)
                                                                   gTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAgt                                                     >  2:269012/1‑135 (MQ=255)
                                                                   gTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAgt                                                     >  1:88869/1‑135 (MQ=255)
                                                                       tgGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAgtggtg                                                 >  6:158105/1‑135 (MQ=255)
                                                                        gggTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAgtggtg                                                 <  3:164165/134‑1 (MQ=255)
                                                                        gggTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAgtggtg                                                 <  7:32790/134‑1 (MQ=255)
                                                                          gTGAAGGTGGCTGTAGCCATCGCGCTGGGGCCGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGACGAAGTGGTGGct                                              >  1:121803/1‑134 (MQ=255)
                                                                            gAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAGTGGTGGCGAt                                            >  4:11105/1‑135 (MQ=255)
                                                                            gAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACACACCTCAGCTATGCCCAAGGCGCATCGGCCGAAGTGGTGGCGAt                                            >  1:175152/1‑135 (MQ=255)
                                                                            gAAGGTGGCTGTAGCCATCGCGCGGGGGCTGGGGACCATGGCTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAGTGGGGGCGAt                                            <  3:291458/135‑1 (MQ=255)
                                                                               ggTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAGTGGTGGCGAtgt                                          >  4:153899/1‑134 (MQ=255)
                                                                                      gTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAGTGGTGGCGATGTGCACCAtt                                  >  2:244142/1‑135 (MQ=255)
                                                                                      gTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAGTGGTGGCGATATGCACCAtt                                  >  1:136392/1‑135 (MQ=255)
                                                                                        atccATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCTTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAGTGGTGGCGATGTGCACCATTg                                 <  8:198632/132‑1 (MQ=255)
                                                                                        aGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAGTGGTGGCGATGTGCACCATTg                                 <  3:153206/134‑1 (MQ=255)
                                                                                            aTCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAGTGGTGGCGATGTGCACCATTGGcgcg                            >  8:83567/1‑135 (MQ=255)
                                                                                                    ggggctggggACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAGTGGTGGCGATGTGCACCATTGGCGCGGCGGatat                    >  5:279712/1‑135 (MQ=255)
                                                                                                            ggACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAGTGGTGGCGATGTGCACCATTGGCGCGGCGGATATGTTCGGGc            <  3:11105/135‑1 (MQ=255)
                                                                                                             gACCATGGTTGGCTGGCGACGTATTGTAGTCACCGtt                                                                                                             >  5:301872/1‑37 (MQ=255)
                                                                                                             gACCATGGTTGGCTGGCGACGTATTGTAGTCACCGtt                                                                                                             <  6:301872/37‑1 (MQ=255)
                                                                                                              aCCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAGTGGTGGCGATGTGCACCATTGGCGCGGCGGATATGTTCGGGCTg          >  8:94766/1‑135 (MQ=255)
                                                                                                                    gTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCgg                                                            <  1:70255/79‑1 (MQ=255)
                                                                                                                    gTTGGCTGGAGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCgg                                                            >  2:70255/1‑79 (MQ=255)
                                                                                                                      tggctggcGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAGTGGTGGCGATGTGCACCATTGGCGCGGCGGATATGTTCGGGCTGCCGGTGTc  <  5:203057/135‑1 (MQ=255)
                                                                                                                          tggcGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAgtg                                                    <  2:151306/81‑1 (MQ=255)
                                                                                                                          tggcGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAgtg                                                    >  1:151306/1‑81 (MQ=255)
                                                                                                                             |                                                                                                                               
ATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGCGACGTATTGTAGTCACCGTTGGCGAGAAGATCGGCAAGACCCACCTCAGCTATGCCCAAGGCGCATCGGCCGAAGTGGTGGCGATGTGCACCATTGGCGCGGCGGATATGTTCGGGCTGCCGGTGTC  >  ACB122/4630099‑4630351

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