Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,992,392 G→T 25.6% D414Y (GAT→TAT)  cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,992,3920GT25.6% 111.8 / 36.4 51D414Y (GAT→TAT) cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (22/16);  new base T (7/6);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

CTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTC  >  NC_000913/3992245‑3992537
                                                                                                                                                   |                                                                                                                                                 
cTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGa                                                                                                                                                  <  1:110328/149‑1 (MQ=255)
     aaTGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGTATATCg                                                                                                                                             >  2:153434/1‑149 (MQ=255)
     aaTGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCg                                                                                                                                             >  4:470477/1‑149 (MQ=255)
         cTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGTATATCGGCGt                                                                                                                                         <  4:49159/149‑1 (MQ=255)
            gATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCt                                                                                                                                      >  1:388405/1‑149 (MQ=255)
               aaCCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGAc                                                                                                                                   <  3:470477/149‑1 (MQ=255)
                           tGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGTATATCGGCGTGCTGACGCGTAAGCTGTa                                                                                                                       >  2:391557/1‑149 (MQ=255)
                                gATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCATCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCcg                                                                                                                  >  2:467202/1‑149 (MQ=255)
                                    gATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGtt                                                                                                              <  3:443594/149‑1 (MQ=255)
                                       cAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGa                                                                                                           <  2:566627/149‑1 (MQ=255)
                                        aGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGaa                                                                                                          <  2:446873/149‑1 (MQ=255)
                                         ggTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGTATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAg                                                                                                         >  3:249541/1‑149 (MQ=255)
                                         ggTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAg                                                                                                         <  2:432044/149‑1 (MQ=255)
                                               tGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTAc                                                                                                   <  1:68320/149‑1 (MQ=255)
                                                 aGGCGCACAATGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTGCGCCGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCa                                                                                                 >  4:215176/1‑149 (MQ=255)
                                                   gcgcACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAgg                                                                                               >  3:5942/1‑149 (MQ=255)
                                                        caACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGTATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAg                                                                                          <  4:276098/149‑1 (MQ=255)
                                                         aaCGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGt                                                                                         >  3:204548/1‑149 (MQ=255)
                                                          aCGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTg                                                                                        >  3:414876/1‑149 (MQ=255)
                                                            gAGTTGCTCGACGCGATTATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCTTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATTCCGCGTTTGAAGCATTACCAGGTAAAGTGAc                                                                                      <  2:348582/149‑1 (MQ=255)
                                                                 gCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGTATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTgg                                                                                 <  1:153434/149‑1 (MQ=255)
                                                                 gCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTgg                                                                                 <  1:317385/149‑1 (MQ=255)
                                                                   tCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTa                                                                               >  1:523126/1‑149 (MQ=255)
                                                                     gACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTaaa                                                                             >  1:183467/1‑149 (MQ=255)
                                                                      aCGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAAc                                                                            <  3:534364/149‑1 (MQ=255)
                                                                          gatgatGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGaa                                                                                                          <  4:278677/115‑1 (MQ=255)
                                                                          gatgatGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGaa                                                                                                          >  3:278677/1‑115 (MQ=255)
                                                                             gatgCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGTATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGa                                                                     <  2:4635/149‑1 (MQ=255)
                                                                             gatgCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGTATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGa                                                                     <  2:278452/149‑1 (MQ=255)
                                                                                        aCCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGAt                                                          <  2:544957/149‑1 (MQ=255)
                                                                                         ccGTAATCTGATCCGCTTTGAGCGTCGCAATAACCTTAGCGTCTACGCCAGTCCGCAGTATATCGGCGTGCTGACGCGTAAGCTGTGTGCCGCGTTTGAAGCTTTACTAGGTTAAGTGACGCTGGTAAAGCTGCAGATTTCACCCGAtt                                                         >  2:181641/1‑148 (MQ=255)
                                                                                             aaTCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGAtctctc                                                     <  4:414876/149‑1 (MQ=255)
                                                                                                  gATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAAc                                                <  2:136356/149‑1 (MQ=255)
                                                                                                  gATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTCTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGAATTCACCCGATCTCTCGGAAc                                                >  2:192345/1‑149 (MQ=255)
                                                                                                   aTCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGTATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAAcc                                               >  1:75856/1‑149 (MQ=255)
                                                                                                   aTCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAAcc                                               >  1:118819/1‑149 (MQ=255)
                                                                                                   aTCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAAcc                                               >  3:85629/1‑149 (MQ=255)
                                                                                                   aTCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACACGATCTCTCGGAAcc                                               >  3:200051/1‑149 (MQ=255)
                                                                                                       gCTTTGCGCGTCGCAATAACCTTAGCGTCCCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAAt                                           <  4:232245/149‑1 (MQ=255)
                                                                                                                 tCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGTATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCttta                                 <  4:472676/149‑1 (MQ=255)
                                                                                                                           cTTAGCGTCTCCGCCAGTCCGCAGTATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTgcc                       >  4:174171/1‑149 (MQ=255)
                                                                                                                           cTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCATGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTgcc                       >  2:568662/1‑149 (MQ=255)
                                                                                                                               gCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGccgccg                   >  2:321609/1‑149 (MQ=255)
                                                                                                                                  tctcCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCgccgggc                >  3:196584/1‑149 (MQ=255)
                                                                                                                                    tcCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCtttatt                               >  2:28748/1‑132 (MQ=255)
                                                                                                                                    tcCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCtttatt                               <  1:28748/132‑1 (MQ=255)
                                                                                                                                      cGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGccgggccggg            >  2:243594/1‑149 (MQ=255)
                                                                                                                                      cGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGccgggccggg            >  1:320096/1‑149 (MQ=255)
                                                                                                                                          aGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTaa        >  2:350295/1‑149 (MQ=255)
                                                                                                                                               gCAGTATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGTACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGtt   >  4:473466/1‑149 (MQ=255)
                                                                                                                                                cAGGATATCGGCGTGCTGACGCGTAGGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTc  >  2:153915/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTC  >  NC_000913/3992245‑3992537

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑