Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,401,040 C→T 11.7% G245E (GGA→GAA)  abgT ← p‑aminobenzoyl‑glutamate transporter; membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,401,0400CT11.7% 87.6 / 3.8 34G245E (GGA→GAA) abgTp‑aminobenzoyl‑glutamate transporter; membrane protein
Reads supporting (aligned to +/- strand):  ref base C (15/15);  new base T (2/2);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

TCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAG  >  NC_000913/1400903‑1401175
                                                                                                                                         |                                                                                                                                       
tCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTTCCTGCCATTGa                                                                                                                              >  4:222747/1‑149 (MQ=255)
tCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTTCCTGCCATTGa                                                                                                                              >  4:218599/1‑149 (MQ=255)
  ggATCGCGCAATATCCCGTTTTGCGGTATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGat                                                                                                                             <  1:243006/148‑2 (MQ=255)
         tcAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGGCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGt                                                                                                                     <  1:265128/148‑1 (MQ=255)
               tCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTc                                                                                                                >  6:160873/1‑148 (MQ=255)
                  cGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTTCCTGCCATTGACCTAACCGTGGCTCgatg                                                                                                            <  3:218599/149‑1 (MQ=255)
                  cGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTTCCTGCCATTGACCTAACCGTGGCTCgatg                                                                                                            <  3:222747/149‑1 (MQ=255)
                         cGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGt                                                                                                     <  6:33567/149‑1 (MQ=255)
                            gATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCgg                                                                                                  <  5:160873/149‑1 (MQ=255)
                            gATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCgg                                                                                                  <  3:269961/149‑1 (MQ=255)
                                 ccATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTAt                                                                                              >  5:132141/1‑148 (MQ=255)
                                   aTCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAg                                                                                           <  8:14232/149‑1 (MQ=255)
                                                cagcaATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGgccgcc                                                                                     >  3:231153/1‑142 (MQ=255)
                                                cagcaATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGgccgcc                                                                                     <  4:231153/142‑1 (MQ=255)
                                                           aGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGAc                                                                   <  6:291968/149‑1 (MQ=255)
                                                                           cTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATc                                                                                             <  3:236310/107‑1 (MQ=255)
                                                                           cTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATc                                                                                             >  4:236310/1‑107 (MQ=255)
                                                                           cTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCAt                                                    >  6:92212/1‑148 (MQ=255)
                                                                           cTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCAt                                                    >  6:317260/1‑148 (MQ=255)
                                                                             gCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATaa                                                  >  3:82216/1‑148 (MQ=255)
                                                                                 tGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATAc                                             >  5:304120/1‑149 (MQ=255)
                                                                                  gCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATAc                                             >  8:15875/1‑148 (MQ=255)
                                                                                   cGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATAc                                             >  2:259418/1‑147 (MQ=255)
                                                                                       aaCCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTa                                       >  8:223319/1‑149 (MQ=255)
                                                                                         ccAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTCACCTAACCGTGGCTCGATGATTTTGTCGGATATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATc                                     <  5:317260/149‑1 (MQ=255)
                                                                                                 cTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTacac                             <  2:136163/149‑1 (MQ=255)
                                                                                                 cTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTacac                             <  7:141276/149‑1 (MQ=255)
                                                                                                  tGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACt                            >  8:327777/1‑149 (MQ=255)
                                                                                                            tCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCAtt                  >  8:147587/1‑149 (MQ=255)
                                                                                                                  tCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGa            <  1:259418/149‑1 (MQ=255)
                                                                                                                       aGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGaa       <  3:14778/149‑1 (MQ=255)
                                                                                                                        gTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAAc      >  2:299803/1‑149 (MQ=255)
                                                                                                                            ctcATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGc                     <  6:27372/130‑1 (MQ=255)
                                                                                                                            ctcATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAg  >  4:90649/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                       
TCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAG  >  NC_000913/1400903‑1401175

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