Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,534,743 G→A 7.9% A240A (GCG→GCA)  nhoA → N‑hydroxyarylamine O‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,534,7430GA 7.9% 150.7 / 2.9 51A240A (GCG→GCA) nhoAN‑hydroxyarylamine O‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (20/26);  new base A (2/2);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

GCACTGGCCGCAGTCGCATTTTCGCCATCATTTGCTGATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTG  >  NC_000913/1534599‑1534883
                                                                                                                                                |                                                                                                                                            
gcACTGGCCGCAGTCGCATTTTCGCCATCATTTGCTGATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCAt                                                                                                                                          >  1:198947/1‑149 (MQ=255)
 cACTGGACGCAGTGGCATTTTCGCCATCATTTGCTTATGTGCCGCCATTTGCCTGACGGCGGAAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCAtt                                                                                                                                         <  2:183535/149‑1 (MQ=255)
      gCCGCAGTCGCATTTTCGCCATCATTTGCTGATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATg                                                                                                                                    >  4:200302/1‑149 (MQ=255)
        cGCAGTCGCATTTTCGCCATCATTTGCTGATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCt                                                                                                                                  <  2:26286/149‑1 (MQ=255)
        cGCAGTCCCATTTTCGCCATCATTTGCTGCAGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCt                                                                                                                                  <  4:35871/149‑1 (MQ=255)
               gCATTTTCGCCATCATTTGCTGATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGc                                                                                                                           <  7:206491/149‑1 (MQ=255)
                 aTTTTCGCCATCATTTGCTGATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCaa                                                                                                                         <  5:69710/149‑1 (MQ=255)
                     tCGCCATCATTTGCTGATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGaa                                                                                                                      <  8:28158/148‑1 (MQ=255)
                       gCCATCATTTGCTGATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAg                                                                                                                   >  1:1402/1‑149 (MQ=255)
                        ccATCATTTGCTGATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGt                                                                                                                  >  8:111380/1‑149 (MQ=255)
                         aataatTTGCTGATGTGCCGCCATTTGCCGGACTGAGGCATGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATTGGCACGCGGTTGAGCAGCGAAATCTACCGGATGTGGCGTCTTTATATGCTGTGATGCAAGAACAGtt                                                                                                                 <  1:154709/145‑1 (MQ=255)
                          atcatTTGCTGATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGttt                                                                                                                <  7:79187/149‑1 (MQ=255)
                           tcatTTGCTGATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTg                                                                                                               <  7:111380/149‑1 (MQ=255)
                           tcatTTGCTGATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTg                                                                                                               >  6:151878/1‑149 (MQ=255)
                             atTTGCTGATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGt                                                                                                             <  2:1402/149‑1 (MQ=255)
                                  cTGATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTgg                                                                                                         <  2:198947/148‑1 (MQ=255)
                                    gATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCATCATTATATGTTGTGATGc                                                                                                                           >  6:8524/1‑128 (MQ=255)
                                    gATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCATCATTATATGTTGTGATGc                                                                                                                           >  6:11666/1‑128 (MQ=255)
                                    gATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCATCATTATATGTTGTGATGc                                                                                                                           <  5:8524/128‑1 (MQ=255)
                                    gATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCATCATTATATGTTGTGATGc                                                                                                                           <  5:11666/128‑1 (MQ=255)
                                       gtgCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTg                                                                                                    <  2:86675/148‑1 (MQ=255)
                                        tgCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGga                                                                                                  <  6:209893/149‑1 (MQ=255)
                                           cgccATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGTAACGCGGTGGAGCAGCTAAATCTACCGGATGTTGCTTCATTATATGCTGTGATTCAAGAACAGTTTGGTCTGGGAGTGgatga                                                                                               >  6:66334/1‑149 (MQ=255)
                                              cATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGc                                                                                            <  5:151878/149‑1 (MQ=255)
                                              cATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGc                                                                                            >  1:205592/1‑149 (MQ=255)
                                                tttGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCg                                                                                           >  7:38535/1‑148 (MQ=255)
                                                   gCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGaaa                                                                                        <  4:170607/148‑1 (MQ=255)
                                                    ccGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAAc                                                                                       <  3:200302/148‑1 (MQ=255)
                                                     cGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGgatgat                                                                                              >  3:125080/1‑140 (MQ=255)
                                                     cGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGgatgat                                                                                              <  4:125080/140‑1 (MQ=255)
                                                     cGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACAt                                                                                     >  4:6335/1‑149 (MQ=255)
                                                       gACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATgg                                                                                   <  4:153501/149‑1 (MQ=255)
                                                              gcAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTAcc                                                                            <  3:133922/149‑1 (MQ=255)
                                                                   cTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGa                                                                       <  7:51681/149‑1 (MQ=255)
                                                                    tGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGAt                                                                      >  3:158932/1‑149 (MQ=255)
                                                                    tGACACTGACCAATTATCATTTTACCCATAATGAAAATGGGCACGCGGTGGAGCACAGAAATCTACCGGATCTGGCGTCAAAATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGTATGATGCGAAACATGGCATTACCGTGCAt                                                                      >  4:37987/1‑149 (MQ=255)
                                                                       cacTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAg                                                                   <  3:127268/149‑1 (MQ=255)
                                                                              cAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAgcgc                                                            >  8:229185/1‑149 (MQ=255)
                                                                               aaTTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGACCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCTAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCt                                                           <  5:232986/149‑1 (MQ=255)
                                                                                   ttCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTg                                                       >  5:97918/1‑149 (MQ=255)
                                                                                      aTTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATg                                                    <  2:205592/149‑1 (MQ=255)
                                                                                             ccATTATGAAAATGGGCACGCGGTTGAGCAGCGAAATCTACCGTATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTTGTCTGGGCGGGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCTGCGt                                             <  5:66334/149‑1 (MQ=255)
                                                                                             ccATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGt                                             >  5:91927/1‑149 (MQ=255)
                                                                                                 tATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGt                                                                          >  6:203610/1‑116 (MQ=255)
                                                                                                 tATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGt                                                                          <  5:203610/116‑1 (MQ=255)
                                                                                                     aaaaTGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACg                                     >  7:139828/1‑149 (MQ=255)
                                                                                                              aCGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATTAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAg                            >  8:239532/1‑149 (MQ=255)
                                                                                                                          agcGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAaataa                <  8:38535/149‑1 (MQ=255)
                                                                                                                                 tCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGt         >  5:202367/1‑149 (MQ=255)
                                                                                                                                   tACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGt         >  4:69496/1‑147 (MQ=255)
                                                                                                                                        gATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTg  >  3:202366/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                            
GCACTGGCCGCAGTCGCATTTTCGCCATCATTTGCTGATGTGCCGCCATTTGCCGGACGGCGGCAAGCTGACACTGACCAATTTTCATTTTACCCATTATGAAAATGGGCACGCGGTGGAGCAGCGAAATCTACCGGATGTGGCGTCATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTG  >  NC_000913/1534599‑1534883

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