Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,759,486 A→G 15.4% intergenic (+84/‑55) fpvA_3 → / → vanA Ferripyoverdine receptor/vanillate demethylase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,759,4860AG15.4% 59.7 / 5.6 26intergenic (+84/‑55)fpvA_3/vanAFerripyoverdine receptor/vanillate demethylase A
Reads supporting (aligned to +/- strand):  ref base A (9/13);  new base G (2/2);  total (11/15)
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.05e-01

TATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATTGTGAGCGGATAACAATTTCACACTCTAGAGAGGAGGACAGCTATGTACCCCAAAAACACCTGGTACGTCGCCTGCACCCCCGATGAGATCGCCACCAAACCCCTGGGCCGG  >  ACB122/4759355‑4759609
                                                                                                                                   |                                                                                                                           
tATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGGTTCTGATCGCCTTAGTTTGAACGCAAAAAACCGCACCGAGGTGCGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAAtgt                                                                                                                         <  1:130531/136‑1 (MQ=255)
          cACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTTGAGCTGTTGACAATTAATCATCGGCTCGTTTAAAGTTTGGaaattt                                                                                                               >  2:94579/1‑132 (MQ=255)
          cACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATtgt                                                                                                               >  3:129589/1‑136 (MQ=255)
          cACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATtgt                                                                                                               >  8:32129/1‑136 (MQ=255)
          cACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATtgt                                                                                                               >  4:102608/1‑136 (MQ=255)
                         tCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATTGTGAGCGGATAACAAtt                                                                                                <  2:160157/136‑1 (MQ=255)
                         tCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGGGCGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATTGTGAGCGGATAACAAtt                                                                                                <  2:35456/136‑1 (MQ=255)
                           aCCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATTGTGAGCGGATAACAATTTc                                                                                              <  5:34771/136‑1 (MQ=255)
                           aCCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATTGTGAGCGGATAACAATTTc                                                                                              <  1:11419/136‑1 (MQ=255)
                           aCCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATTGTGAGCGGATAACAATTTc                                                                                              <  4:111314/136‑1 (MQ=255)
                                            cTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATTGTGAGCGGATAACAATTTCACACTCTAGAgagga                                                                               >  6:21669/1‑134 (MQ=255)
                                            cTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATTGTGAGCGGATAACAATTTCACACTCTAGAgagga                                                                               <  3:103212/134‑1 (MQ=255)
                                            cTGATCGCCTGACGTTGAACGAAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATTGTGAGCGGATAACAATTTCACACTCTAGAgagga                                                                               <  5:21669/134‑1 (MQ=255)
                                             tGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGGGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATTGTGAGCGGATAACAATTTCACACTCTAGAGAggagga                                                                            <  4:129589/136‑1 (MQ=255)
                                               aTCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATTGTGAGCGGATAACAATTTCACACTCTAGAGAggagga                                                                            <  2:10083/134‑1 (MQ=255)
                                                   ccTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATTGTGAGCGGATAACAATTTCACACTCTAGAGAGGAGGACAGCTa                                                                      <  3:187319/136‑1 (MQ=255)
                                                            gaaccaaaaaacCCGCACCCTGGTGCGTTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATTGTGAGCGGATAACAATTTCACACTCTAGAGAGGAGGACAGCTATGTAcccc                                                              <  1:94579/129‑1 (MQ=255)
                                                                                 ggTGCGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATTGTGAGCGGATAACAATTTCACACTCTAGAGAGGAGGACAGCTATGTACCCCAAAAACACCTGGTACGTCGCCt                                        <  6:171515/136‑1 (MQ=255)
                                                                                     cGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATTGTGAGCGGATAACAATTTCACACTCTAGAGAGGAGGACAGCTATGTACCCCAAAAACACCTGGTACGTCGCCTGCAc                                    >  1:165782/1‑136 (MQ=255)
                                                                                        ttttttGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATTGTGAGCGGATAACAATTTCACACTCTAGAGAGGAGGACAGCTATGTACCCCAAAAACACCTGGTACGTCGCCTGCAcccc                                 >  8:23878/1‑136 (MQ=255)
                                                                                                                       atcGGCTCGTATAATGTGTGGAATTGTGAGCGGATAACAATTTCACACTCTAGAGAGGAGGACAGCTATGTACCCCAAAAACACCTGGTACGTCGCCTGCACCCCCGATGAGATCGCCACCAACCCCCTGGGACgg  >  6:107749/1‑136 (MQ=255)
                                                                                                                            cTCGTATGTTGTGTGGAATTGTGAGCGGATAACAAttt                                                                                               >  1:133399/1‑38 (MQ=255)
                                                                                                                            cTCGTATGTTGTGTGGAATTGTGAGCGGATAACAAttt                                                                                               <  2:133399/38‑1 (MQ=255)
                                                                                                                            cTCGTATAATGTGTGGAATTGTGAGCGGATAACAATTTCACACTCTAGAGAGTAGGACAGCTATGTACCCCAAAAACACCTGGTACGTCGCCTGCACCCCCGATGAGATCGCCACCACACCCCTGGGCCgg  >  7:89248/1‑131 (MQ=255)
                                                                                                                              cGTATGTTGTGTGGAATTGTGAGCGGATAACAAttt                                                                                               <  2:123148/36‑1 (MQ=255)
                                                                                                                              cGTATGTTGTGTGGAATTGTGAGCGGATAACAAttt                                                                                               >  1:123148/1‑36 (MQ=255)
                                                                                                                                   |                                                                                                                           
TATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCGAGCTGTTGACAATTAATCATCGGCTCGTATAATGTGTGGAATTGTGAGCGGATAACAATTTCACACTCTAGAGAGGAGGACAGCTATGTACCCCAAAAACACCTGGTACGTCGCCTGCACCCCCGATGAGATCGCCACCAAACCCCTGGGCCGG  >  ACB122/4759355‑4759609

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: