Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,482,560 T→A 8.2% I101F (ATC→TTC)  azoR ← NADH‑azoreductase, FMN‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,482,5600TA 8.2% 143.9 / 4.2 49I101F (ATC→TTC) azoRNADH‑azoreductase, FMN‑dependent
Reads supporting (aligned to +/- strand):  ref base T (20/25);  new base A (2/2);  total (22/27)
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.91e-01

TCCATCTTTGTGGATCCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCGGAAAGTGCCAGAGCTTCCTGCTGACGCGGAGTCAGCGGCGCATCGCTCGGACGCAGAGCGCCAACCAGTTCGCC  >  NC_000913/1482417‑1482704
                                                                                                                                               |                                                                                                                                                
tCCATCTTTGTGGATCCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGa                                                                                                                                             <  1:236974/149‑1 (MQ=255)
 ccATCTTTGTGGATCCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTg                                                                                                                                              >  8:101282/1‑147 (MQ=255)
   aTCTTTGTGGATCCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGt                                                                                                                                          <  5:54884/149‑1 (MQ=255)
     cTTTGTGGATCCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGTTGTTGAAGTta                                                                                                                                        >  5:172555/1‑149 (MQ=255)
            gATCCCGCCGCGGCTGGTAATAACGATGGCTTTTATACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCg                                                                                                                                 >  1:127702/1‑149 (MQ=255)
             aTCCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCgg                                                                                                                                >  8:103374/1‑149 (MQ=255)
                ccgccgCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGAAGTTGAAGTTATACATCGGTGc                                                                                                                             >  1:113789/1‑149 (MQ=255)
                ccgccgCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGAAGTTGAAGTTATACATCGGTGc                                                                                                                             >  1:110337/1‑149 (MQ=255)
                        gCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAataa                                                                                                                     <  1:23153/149‑1 (MQ=255)
                            gTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAataacg                                                                                                                   <  5:49510/147‑1 (MQ=255)
                            gTAATAACGATGGC‑TTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTtat                                                                                                                                       >  8:117065/1‑126 (MQ=255)
                            gTAATAACGATGGC‑TTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTtat                                                                                                                                       <  7:117065/126‑1 (MQ=255)
                             taataaCGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCTGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAataacgata                                                                                                                <  2:235686/149‑1 (MQ=255)
                                         cTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCt                                                                                                    >  4:180649/1‑149 (MQ=255)
                                          ttttttACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAataacgat                                                                                                                 >  8:245021/1‑135 (MQ=255)
                                          ttttttACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAataacgat                                                                                                                 <  7:245021/135‑1 (MQ=255)
                                            ttttACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCttt                                                                                                  <  1:140090/148‑1 (MQ=255)
                                            ttttACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCtt                                                                                                   <  1:131872/147‑1 (MQ=255)
                                             tttACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCttt                                                                                                  >  6:58402/1‑147 (MQ=255)
                                             tttACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGAAGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTc                                                                                                 <  2:110337/148‑1 (MQ=255)
                                             tttACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGAAGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTc                                                                                                 <  2:113789/148‑1 (MQ=255)
                                              ttACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAg                                                                                               >  3:44885/1‑149 (MQ=255)
                                                  ccGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCt                                                                                             >  5:24354/1‑147 (MQ=255)
                                                           gACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAAc                                                                                  <  7:103374/149‑1 (MQ=255)
                                                               ttCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTg                                                                                                        >  1:74945/1‑123 (MQ=255)
                                                               ttCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTg                                                                                                        <  2:74945/123‑1 (MQ=255)
                                                                  cGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCTACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCg                                                                            <  8:89156/148‑1 (MQ=255)
                                                                   ggACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCTATCAACTCATCgg                                                                           >  5:220111/1‑148 (MQ=255)
                                                                   ggACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCgg                                                                           >  3:104221/1‑148 (MQ=255)
                                                                   ggACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCgg                                                                           >  7:245492/1‑148 (MQ=255)
                                                                       cGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCGGAAAGt                                                                      <  6:24354/149‑1 (MQ=255)
                                                                         ttCTAGGTATAGCGGAAAGTAACGCCTGCGCGGGAAACCAGGTCAAAATAGTTTTTCAACTGAGTCGATATGTTGAAGTTATACATCTGTTCCGCAATAACGATAACGTCGTTGGCTTTCAGCGCGGCAATCAACTCATCGGAAAGTGc                                                                    <  3:34097/149‑1 (MQ=255)
                                                                            tcGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCGGAAAGTGc                                                                    >  4:196110/1‑146 (MQ=255)
                                                                                  taGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCGGAAAGTGCCAGAGCtt                                                            <  1:13858/148‑1 (MQ=255)
                                                                                     cGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCGGAAAGTGCCAGAGCTTCct                                                         <  3:180649/148‑1 (MQ=255)
                                                                                      ggAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCGGAAAGTGCCAGAGCTTCctgc                                                       <  4:44885/149‑1 (MQ=255)
                                                                                             aaCGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCGGAAAGTGCCAGAGCTTCCTGCTGAcgcg                                                >  4:132204/1‑149 (MQ=255)
                                                                                              aCGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCGGAAAGTGCCAGAGCTTCCTGCTGAcgcg                                                >  7:195486/1‑148 (MQ=255)
                                                                                               cGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCGGAAAGTGCCAGAGCTTCCTGCTGACGCGGa                                              <  4:28953/149‑1 (MQ=255)
                                                                                                 ccTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCGGAAAGTGCCAGAGCTTCCTGCTGACGCGGAGt                                            <  2:127702/149‑1 (MQ=255)
                                                                                                       cgGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCGGAAAGTGCCAGAGCTTCCTGCTGACTCCGCGTCAgcag                                      <  4:104221/149‑3 (MQ=255)
                                                                                                               cAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATACCGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCGGAAAGTGCCAGAGCTTCCTGCTGACGCGGAGTCAGCGGCGCATCGc                              <  3:132204/149‑1 (MQ=255)
                                                                                                                    cAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCGGAAAGTGCCAGAGCTTCCTGCTGACGCGGAGTCAGCGGCGCATCGCTCGGa                         >  5:142952/1‑149 (MQ=255)
                                                                                                                     aaaaTAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAAc                                                                                  <  6:212654/91‑1 (MQ=255)
                                                                                                                     aaaaTAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAAc                                                                                  >  5:212654/1‑91 (MQ=255)
                                                                                                                                tCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCGGAAAGTGCCAGAGCTTCCTGCTGACGCGGAGTCAGCGGCGCATCGCTCGGACGCAGAGCGcc              <  1:150176/148‑1 (MQ=255)
                                                                                                                                         tCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCGGAAAGTGCCAGAGCTTCCTGCTGACGCGGAGTCAGCGGCGCATCGCTCGGACGCAGAGCGCCAACCAGTTCg    >  5:47293/1‑149 (MQ=255)
                                                                                                                                          cGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCGTTCAGCGCGGCAATCAACTCATCGGAAAGTGCCAGAGCTTCCTGCTGACGCGGAGTCAGCGGCGCATCGCTCGGACGCAGAGCGCCAACCAGTTCGc   <  3:164435/149‑1 (MQ=255)
                                                                                                                                           gagaTGTTGAAGTTATACATCGGTGCCTCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCGGAAAGTGCCAGAGCTTCCTGCTGACGCGGAGTCAGCGGCGCATCGCTCGGACGCAGAGCGCCAACCAGTTCGcc  <  2:86553/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
TCCATCTTTGTGGATCCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAGCGGAAAGTAACGCCTGCGCGGGCAACCAGGTCAAAATAGTTTTTCAACTGAGTCGAGATGTTGAAGTTATACATCGGTGCCGCAATAACGATAACGTCGTGGGCTTTCAGCTCGGCAATCAACTCATCGGAAAGTGCCAGAGCTTCCTGCTGACGCGGAGTCAGCGGCGCATCGCTCGGACGCAGAGCGCCAACCAGTTCGCC  >  NC_000913/1482417‑1482704

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