Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,401,844 A→G 100% I471T (ATC→ACC)  abgB ← p‑aminobenzoyl‑glutamate hydrolase, B subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,401,8440AG97.8% 140.0 / ‑4.6 45I471T (ATC→ACC) abgBp‑aminobenzoyl‑glutamate hydrolase, B subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (20/23);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

GACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTT  >  NC_000913/1401751‑1401939
                                                                                             |                                                                                               
gacgGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGc                                                                                           <  1:1772779/100‑1 (MQ=255)
   ggTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGt                                                                                        <  2:1231985/100‑1 (MQ=255)
    gTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTg                                                                                       <  2:4209/100‑1 (MQ=255)
       tGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATa                                                                                    >  2:1016794/1‑100 (MQ=255)
       tGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATa                                                                                    <  2:293146/100‑1 (MQ=255)
        ggATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATAc                                                                                   >  2:1883525/1‑100 (MQ=255)
        ggATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATAc                                                                                   >  2:1096419/1‑100 (MQ=255)
            gACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGtt                                                                               >  2:1865002/1‑100 (MQ=255)
            gACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGtt                                                                               <  1:351798/100‑1 (MQ=255)
              cATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGt                                                                                >  1:1858537/1‑97 (MQ=255)
                tACTCATGGGCATTCCTCGGTTGGTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGTGTAGGGCAGTGATACGGTTGCg                                                                            <  2:1863346/99‑1 (MQ=255)
                       gggCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGtt                                                                    <  2:53395/100‑1 (MQ=255)
                              ccTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCt                                                             >  2:1235155/1‑100 (MQ=255)
                               cTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGTTAGGGCAGTGATACGGATGCGTGTCCGTTACTTGCTg                                                            >  1:149361/1‑100 (MQ=255)
                                 cGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGTGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGAt                                                          <  2:491133/100‑1 (MQ=255)
                                     tgttgtGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGAAGGGCAGTGATACGGCTGCGTGTCCGTTACTTGCTGATGctc                                                      >  2:1113210/1‑100 (MQ=255)
                                      gttgtGTTTGCGTTTTTGTTGTTTTTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGctct                                                     <  2:1071201/100‑1 (MQ=255)
                                           gtTTGCGTGGTTGTTGTTATTTTAAAGGTGACTGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCttgtt                                                 <  2:420338/99‑1 (MQ=255)
                                                     ttgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGc                                      <  1:106277/100‑1 (MQ=255)
                                                      tgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGcg                                     >  1:262252/1‑99 (MQ=255)
                                                        ttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAAt                                   <  1:1392572/100‑1 (MQ=255)
                                                        ttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAAt                                   >  1:1993292/1‑100 (MQ=255)
                                                          gttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCGTTCTTGCAATAg                                 >  1:967507/1‑100 (MQ=255)
                                                          gttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAAt                                   >  2:1269050/1‑98 (MQ=255)
                                                           ttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGc                                <  1:860363/100‑1 (MQ=255)
                                                               tttAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTg                            >  1:240470/1‑100 (MQ=255)
                                                               tttAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTg                            <  1:1086936/100‑1 (MQ=255)
                                                                  aaaGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAAt                         >  2:681245/1‑100 (MQ=255)
                                                                   aaGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAAt                         >  1:1870857/1‑99 (MQ=255)
                                                                     ggTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTTGCCTTCTTGCAATAGTCCTGAATCaa                      >  1:235959/1‑100 (MQ=255)
                                                                      gTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAg                     <  2:2381878/100‑1 (MQ=255)
                                                                      gTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAg                     <  1:1311559/100‑1 (MQ=255)
                                                                       tGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAgg                    <  1:1914908/100‑1 (MQ=255)
                                                                        gACGGTGTCACGTTTTTCGGGGTAGGGCAGTTATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAgg                    <  2:2055138/99‑1 (MQ=255)
                                                                            gtgtCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTGGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGaa                  <  1:2042439/97‑1 (MQ=255)
                                                                              gtCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGAt             >  2:1350058/1‑100 (MQ=255)
                                                                              gtCACGTTTTTCGGGGTAGAGCAGAGATACGGTTACATGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAagag               >  2:560544/1‑98 (MQ=255)
                                                                               tCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGAt             <  1:111904/99‑1 (MQ=255)
                                                                                cACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTg           >  1:2041212/1‑100 (MQ=255)
                                                                                cACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTg           <  1:1016805/100‑1 (MQ=255)
                                                                                     ttttCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACtgt      <  2:2044536/100‑1 (MQ=255)
                                                                                       ttCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCGATAGCCCTGAATCAAGGAAGAGATTGACTGTgg    <  2:262249/100‑1 (MQ=255)
                                                                                       ttCGGGATAGGGCAGTGATACGGCTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGGCtgtg     >  2:905872/1‑99 (MQ=255)
                                                                                        tCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGt   >  2:2016228/1‑100 (MQ=255)
                                                                                         cGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGtt  >  2:163685/1‑100 (MQ=255)
                                                                                             |                                                                                               
GACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTT  >  NC_000913/1401751‑1401939

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: