Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,401,844 A→G 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,8440AG94.1% 161.0 / 0.9 51I471T (ATC→ACC) abgBp‑aminobenzoyl‑glutamate hydrolase, B subunit
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (25/23);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

ACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGG  >  NC_000913/1401752‑1401937
                                                                                            |                                                                                             
acgGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCa                                                                                        >  2:754811/1‑100 (MQ=255)
   gTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTg                                                                                     <  2:738802/100‑1 (MQ=255)
   gTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTg                                                                                     >  2:509694/1‑100 (MQ=255)
       ggATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATAc                                                                                 >  1:902961/1‑100 (MQ=255)
                aCTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCgtg                                                                        <  1:642636/100‑1 (MQ=255)
                aCTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCgag                                                                        >  1:1110524/1‑98 (MQ=255)
                aCTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCg                                                                          <  2:1110514/98‑1 (MQ=255)
                   cATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTcc                                                                     >  1:833977/1‑100 (MQ=255)
                   cATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTc                                                                      >  1:994792/1‑99 (MQ=255)
                    aTGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCg                                                                    <  2:1211484/100‑1 (MQ=255)
                     tGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTcc                                                                     <  2:778915/98‑1 (MQ=255)
                      gggCATTCCTCGGTTGTTGTGTTTGCGTTGTTATTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGtt                                                                  >  1:1235673/1‑100 (MQ=255)
                       ggCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTGTTCGGAGTAGGGCAGTGATACGGTTGCGTGTCCGTTa                                                                 >  2:876522/1‑100 (MQ=255)
                         cATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTAGCGTGTCCGTTACt                                                               >  2:887192/1‑100 (MQ=255)
                           ttCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTg                                                             <  2:571944/100‑1 (MQ=255)
                           ttCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTg                                                             >  1:420651/1‑100 (MQ=255)
                           ttCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTg                                                             >  1:1231549/1‑100 (MQ=255)
                           ttCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTg                                                             <  1:865968/100‑1 (MQ=255)
                             ccTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCt                                                           <  2:100062/100‑1 (MQ=255)
                             ccTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCt                                                           >  2:84811/1‑100 (MQ=255)
                              cTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTg                                                          <  2:420647/100‑1 (MQ=255)
                              cTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTg                                                          <  1:104495/100‑1 (MQ=255)
                                cGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGAt                                                        <  1:731980/100‑1 (MQ=255)
                                cGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGAt                                                        <  2:283721/100‑1 (MQ=255)
                                 ggTTGTTGTGTTTGCGATGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATg                                                       >  1:25943/1‑100 (MQ=255)
                                       tgtgTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCttg                                                 <  2:902954/100‑1 (MQ=255)
                                         tgtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCttgtt                                               >  2:456765/1‑100 (MQ=255)
                                         tgtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCttgt                                                >  2:896374/1‑99 (MQ=255)
                                          gtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCttgttg                                              <  2:701826/100‑1 (MQ=255)
                                                 ttgttgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCt                                       <  1:76540/100‑1 (MQ=255)
                                                   gttgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTg                                     <  1:962696/100‑1 (MQ=255)
                                                   gttgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTg                                     >  2:638958/1‑100 (MQ=255)
                                                    ttgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGc                                    <  2:1231538/100‑1 (MQ=255)
                                                     tgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCa                                   <  1:458812/100‑1 (MQ=255)
                                                     tgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCa                                   >  1:63284/1‑100 (MQ=255)
                                                        tgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATa                                >  2:1295364/1‑100 (MQ=255)
                                                         gttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAg                               >  1:620779/1‑100 (MQ=255)
                                                         gttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAg                               >  2:534701/1‑100 (MQ=255)
                                                             ttttAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCt                           >  2:937718/1‑100 (MQ=255)
                                                             ttttAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCt                           >  1:798641/1‑100 (MQ=255)
                                                              tttAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTg                          >  2:543860/1‑100 (MQ=255)
                                                               ttAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGa                         >  1:139016/1‑100 (MQ=255)
                                                                 aaaGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAAt                       >  2:957672/1‑100 (MQ=255)
                                                                   aGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCa                     >  1:1040997/1‑100 (MQ=255)
                                                                     gTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAg                   >  1:59411/1‑100 (MQ=255)
                                                                       gACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCa                     >  2:536622/1‑96 (MQ=255)
                                                                       gACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGa                 <  2:708836/100‑1 (MQ=255)
                                                                       gACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGa                 <  1:1102199/100‑1 (MQ=255)
                                                                         cGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAag               <  1:1144485/100‑1 (MQ=255)
                                                                           gtgtCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAagag             <  1:299555/100‑1 (MQ=255)
                                                                                      ttCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTgg  <  1:896381/100‑1 (MQ=255)
                                                                                            |                                                                                             
ACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGG  >  NC_000913/1401752‑1401937

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