Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,534,897 A→T 8.5% intergenic (+28/+68) nhoA → / ← yddE N‑hydroxyarylamine O‑acetyltransferase/PhzC‑PhzF family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,534,8970AT 8.5% 201.5 / 8.1 59intergenic (+28/+68)nhoA/yddEN‑hydroxyarylamine O‑acetyltransferase/PhzC‑PhzF family protein
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base T (3/2);  total (30/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 = 9.51e-01

ATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACG  >  NC_000913/1534751‑1535045
                                                                                                                                                  |                                                                                                                                                    
atGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTaaa                                                                                                                                                    <  4:852809/149‑1 (MQ=255)
  gCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACg                                                                                                                                                  <  2:28902/149‑1 (MQ=255)
             aGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAAc                                                                                                                                                   <  4:433983/137‑1 (MQ=255)
             aGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAAc                                                                                                                                                   >  3:433983/1‑137 (MQ=255)
             aGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                       >  1:172029/1‑149 (MQ=255)
              gAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCa                                                                                                                                      <  3:1044196/149‑1 (MQ=255)
                 cAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATAc                                                                                                                                   <  2:172029/149‑1 (MQ=255)
                                   gatgatGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                          >  3:911749/1‑124 (MQ=255)
                                   gatgatGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGACGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                          <  4:911749/124‑1 (MQ=255)
                                     tgatgCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCt                                                                                                               <  3:357036/149‑1 (MQ=255)
                                       atgCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTa                                                                                                                                             <  4:38333/117‑1 (MQ=255)
                                       atgCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTa                                                                                                                                             >  3:38333/1‑117 (MQ=255)
                                           gAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCAGACATTAGCCCGGTTGTAGACCTGACTgg                                                                                                         >  3:385450/1‑149 (MQ=255)
                                            aaaCATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGtt                                                                                                                       >  3:1132798/1‑134 (MQ=255)
                                            aaaCATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGtt                                                                                                                       <  4:1132798/134‑1 (MQ=255)
                                            aaaCATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                        >  3:445570/1‑117 (MQ=38)
                                                    cTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACa                                                                                                                                  >  3:372052/1‑115 (MQ=255)
                                                    cTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACa                                                                                                                                  <  4:372052/115‑1 (MQ=255)
                                                       tACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCa                                                                                             >  3:1248786/1‑149 (MQ=255)
                                                          cGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                          <  3:1040657/101‑1 (MQ=255)
                                                          cGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                          >  4:1040657/1‑101 (MQ=255)
                                                              gacgAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACTCGGAGGCAGGAAAATAATGTATGTCAGCTTGCCGAATGCGGCGTAAACGCCTTACCTGGCATACATTAGCCTGGTTGTAGACCTGACAGGCGAAGCGCATAAGATAACg                                                                                      <  4:385450/146‑1 (MQ=255)
                                                                   gTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAAc                                                                                                                                                   >  2:1249848/1‑83 (MQ=255)
                                                                   gTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAAc                                                                                                                                                   <  1:1249848/83‑1 (MQ=255)
                                                                    ttAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCa                                                                                                                                      <  2:630486/95‑1 (MQ=255)
                                                                    ttAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCa                                                                                                                                      >  1:630486/1‑95 (MQ=255)
                                                                    ttAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCa                                                                                <  4:1248786/149‑1 (MQ=255)
                                                                    ttAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCa                                                                                >  1:355242/1‑149 (MQ=255)
                                                                             gTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCaga                                                                              <  2:634364/142‑1 (MQ=255)
                                                                             gTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCaga                                                                              >  1:634364/1‑142 (MQ=255)
                                                                              tGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTaa                                                                                                                                                     >  4:1383925/1‑70 (MQ=255)
                                                                              tGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTaa                                                                                                                                                     <  3:1383925/70‑1 (MQ=255)
                                                                                   gcggcgTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGc                                                                                                                                                 <  2:268377/69‑1 (MQ=255)
                                                                                   gcggcgTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGc                                                                                                                                                 >  1:268377/1‑69 (MQ=255)
                                                                                      gcgTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTc                                                              <  2:355242/149‑1 (MQ=255)
                                                                                                   aCCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCTGGCATACATTa                                                                                                                               <  1:1087735/71‑1 (MQ=255)
                                                                                                   aCCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTa                                                                                                                               >  2:1087735/1‑71 (MQ=255)
                                                                                                                aaaaTAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCt                                                                                                                                               <  1:1138723/42‑1 (MQ=255)
                                                                                                                aaaaTAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCt                                                                                                                                               >  2:1138723/1‑42 (MQ=255)
                                                                                                                  aataatTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGc                                                                    <  1:480744/115‑1 (MQ=255)
                                                                                                                  aataatTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGc                                                                    >  2:480744/1‑115 (MQ=255)
                                                                                                                       tttATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACaaaa                                                >  3:499059/1‑130 (MQ=255)
                                                                                                                       tttATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACaaaa                                                <  4:499059/130‑1 (MQ=255)
                                                                                                                       tttATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAAt                             >  2:1348249/1‑149 (MQ=255)
                                                                                                                       tttATGTCAGGCTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGGCCATTCGGCATGAATCAAAATCACCGCTGTGCCTGAAAt                             >  3:961408/1‑149 (MQ=255)
                                                                                                                           tGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTg                         >  3:952820/1‑149 (MQ=255)
                                                                                                                            gTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCaga                                                                              >  1:1291911/1‑95 (MQ=255)
                                                                                                                            gTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCaga                                                                              <  2:1291911/95‑1 (MQ=255)
                                                                                                                            gTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGa                        <  4:952820/149‑1 (MQ=255)
                                                                                                                             tCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCa                                                                                >  3:461352/1‑92 (MQ=255)
                                                                                                                             tCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCa                                                                                <  4:461352/92‑1 (MQ=255)
                                                                                                                             tCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGAc                       <  1:1348249/149‑1 (MQ=255)
                                                                                                                                    gCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTcc                <  4:961408/149‑1 (MQ=255)
                                                                                                                                      cGGATAAGGCGTTTACGCCGCATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACaaaa                                                <  2:419047/115‑1 (MQ=38)
                                                                                                                                      cGGATAAGGCGTTTACGCCGCATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACaaaa                                                >  1:419047/1‑115 (MQ=38)
                                                                                                                                        gataaggcgTTTACGCCGCATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCt                                  >  4:116578/6‑127 (MQ=255)
                                                                                                                                        gataaggcgTTTACGCCGCATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCt                                  <  3:116578/122‑1 (MQ=255)
                                                                                                                                        gATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACaaa                                                 >  4:1427413/1‑112 (MQ=255)
                                                                                                                                        gATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACaaa                                                 <  3:1427413/112‑1 (MQ=255)
                                                                                                                                              gcgTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTAt      >  3:1392106/1‑149 (MQ=255)
                                                                                                                                               cgTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATc     >  1:113160/1‑149 (MQ=255)
                                                                                                                                                 tAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCAc   >  3:522723/1‑149 (MQ=255)
                                                                                                                                                  aaaCGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACg  >  4:1314782/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACG  >  NC_000913/1534751‑1535045

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