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,8440AG96.7% 199.7 / ‑2.1 61I471T (ATC→ACC) abgBp‑aminobenzoyl‑glutamate hydrolase, B subunit
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (31/28);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTG  >  NC_000913/1401747‑1401940
                                                                                                 |                                                                                                
ggaCGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGGTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTa                                                                                                >  2:70579/1‑100 (MQ=255)
   cgacgGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTGGGGGTAggg                                                                                             >  1:965105/1‑100 (MQ=255)
    gacgGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGc                                                                                            <  2:40208/100‑1 (MQ=255)
      cgGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCa                                                                                           >  2:1307759/1‑99 (MQ=255)
        gTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTg                                                                                        <  2:891473/100‑1 (MQ=255)
            ggATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATAc                                                                                    >  2:329203/1‑100 (MQ=255)
                 aCATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTg                                                                               >  1:289296/1‑100 (MQ=255)
                  cATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGc                                                                              >  2:932532/1‑100 (MQ=255)
                   aTACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCg                                                                             <  2:477350/100‑1 (MQ=255)
                   aTACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCg                                                                             >  2:157211/1‑100 (MQ=255)
                        cATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGCCGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTcc                                                                        >  2:540453/1‑100 (MQ=255)
                        cATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTcc                                                                        >  2:795750/1‑100 (MQ=255)
                         aTGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCg                                                                       <  1:232309/100‑1 (MQ=255)
                         aTGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCg                                                                       <  1:1433911/100‑1 (MQ=255)
                          tGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGt                                                                      <  2:1201966/100‑1 (MQ=255)
                           gggCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGtt                                                                     >  2:526313/1‑100 (MQ=255)
                            ggCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGtt                                                                     >  2:1179323/1‑99 (MQ=255)
                            ggCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTa                                                                    <  2:1257039/100‑1 (MQ=255)
                            ggCATTCCTCGGCTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTa                                                                    >  1:1032309/1‑100 (MQ=255)
                             gCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTAc                                                                   >  2:1099994/1‑100 (MQ=255)
                                  ccTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCt                                                              <  1:24382/100‑1 (MQ=255)
                                   cTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTg                                                             >  1:172886/1‑100 (MQ=255)
                                          gttgtGTTTGCGTTGTTGTTGTTATTTTAAATGTGACGGTGTCACGTTTTTTGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGctct                                                      >  2:1353506/1‑100 (MQ=255)
                                          gttgtGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTATGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGctct                                                      <  2:860244/100‑1 (MQ=255)
                                            tgtgTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCttg                                                    >  2:1066423/1‑100 (MQ=255)
                                               gtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCttgttg                                                 >  1:282395/1‑100 (MQ=255)
                                               gtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCttgttg                                                 <  1:227586/100‑1 (MQ=255)
                                               gtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGGGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCttgttg                                                 <  2:136906/100‑1 (MQ=255)
                                                   gCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCAt                                             >  2:1244696/1‑100 (MQ=255)
                                                         ttgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGc                                       <  2:1223168/100‑1 (MQ=255)
                                                           gttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCaa                                     >  2:102749/1‑100 (MQ=255)
                                                               ttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGc                                 <  1:982763/100‑1 (MQ=255)
                                                                tATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGcc                                >  1:956333/1‑100 (MQ=255)
                                                                tATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGcc                                <  1:271589/100‑1 (MQ=255)
                                                                tATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGcc                                <  2:58891/100‑1 (MQ=255)
                                                                   tttAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTg                             <  2:1509705/100‑1 (MQ=255)
                                                                   tttAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTg                             <  1:1402503/100‑1 (MQ=255)
                                                                   tttAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTg                             >  2:34535/1‑100 (MQ=255)
                                                                   tttAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTg                             <  1:523591/100‑1 (MQ=255)
                                                                   tttAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTg                             >  2:1373133/1‑100 (MQ=255)
                                                                   tttAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTg                             <  1:1419006/100‑1 (MQ=255)
                                                                   tttAAAGGCGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTg                             >  2:477378/1‑100 (MQ=255)
                                                                    ttAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGa                            >  1:54800/1‑100 (MQ=255)
                                                                     tAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTCCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGaa                           <  2:843852/100‑1 (MQ=255)
                                                                      aaaGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAAt                          <  2:289295/100‑1 (MQ=255)
                                                                       aaGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATc                         >  2:217526/1‑100 (MQ=255)
                                                                         ggTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCaa                       >  1:354487/1‑100 (MQ=255)
                                                                          gTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAg                      >  2:991211/1‑100 (MQ=255)
                                                                          gTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAg                      <  2:896112/100‑1 (MQ=255)
                                                                            gACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGa                    >  1:452839/1‑100 (MQ=255)
                                                                                gtgtCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAagag                <  1:369581/100‑1 (MQ=255)
                                                                                  gttaCGTTATTATGGGTAGGGCAGTGATATGGTGGCGTGTGTTTTACTTTTTGATGCTCTTGTTGGCTTTCGTGCAATAGCTCTGAATCATGGAaga                 >  2:681203/4‑97 (MQ=255)
                                                                                  gtCACGTTTTTCGGGGTAGGGCAGTGCTACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGAt              <  1:70579/100‑1 (MQ=255)
                                                                                  gtCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGAt              <  1:157211/100‑1 (MQ=255)
                                                                                    cACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTCGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATc                         >  2:573743/1‑87 (MQ=255)
                                                                                           ttCGGGGTATGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTgg     <  1:448381/100‑1 (MQ=255)
                                                                                           ttCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTgg     <  1:217527/100‑1 (MQ=255)
                                                                                           ttCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTgg     <  1:448310/100‑1 (MQ=255)
                                                                                            tCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGt    >  1:178723/1‑100 (MQ=255)
                                                                                            tCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGt    <  1:985800/100‑1 (MQ=255)
                                                                                              ggggTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGTAAGAGATTGACTGTGGTtg  >  2:690464/1‑100 (MQ=255)
                                                                                              gggATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTtg  >  1:48032/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTG  >  NC_000913/1401747‑1401940

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