Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,464,032 A→G N165S (AAC→AGC)  paaX → DNA‑binding transcriptional repressor PaaX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,464,0320AG96.5% 82.3 / NA 29N165S (AAC→AGC) paaXDNA‑binding transcriptional repressor PaaX
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (8/20);  total (9/20)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTGGCAAGGTTTTGGCGCACTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTACGAAACCTTTATTCAGTCATTCCGCCCGCTGGTG  >  NC_000913/1463898‑1464168
                                                                                                                                      |                                                                                                                                        
aTCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTGGCAAGGTTTTGGCGCACTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTgg                                                                                                                                     >  2:75901/1‑139 (MQ=255)
     cGCTGGCTGATGTCAAAAAGCAGTTGATCTGGCAAGGTTTTGGCGCACTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGt                                                                                                                                >  2:447863/1‑140 (MQ=255)
         ggctgATGTCAAAAAGCAGTTGATCTGGCAAGGTTTTGGCGCACTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGTTTTg                                                                                                                            <  1:4641/140‑1 (MQ=255)
          gctgATGTCAAAAAGCAGTTGATCTGGCAAGGTTTTGGCGCACTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGttttta                                                                                                                           >  2:467220/1‑138 (MQ=255)
          gctgATGTCAAAAAGCAGTTGATCTGGCAAGGTTTTGGCGCACTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGTTTTGa                                                                                                                           >  1:535014/1‑140 (MQ=255)
          gctgATGTCAAAAAGCAGTTGATCTGGCAAGGTTTTGGCGCACTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAAAGTGATTTGttttta                                                                                                                           >  1:585755/1‑138 (MQ=255)
             gATGTCAAAAAGCAGTTGATCTGGCAAGGTTTTGGCGCACTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGTTTTGAAgc                                                                                                                        <  1:503020/140‑1 (MQ=255)
             gATGTCAAAAAGCAGTTGATCTGGCAAGGTTTTGGCGCACTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGTTTTGAAgc                                                                                                                        <  1:499063/140‑1 (MQ=255)
                    aaaaGCAGTTGATCTGGCAAGGTTTTGGCGCACTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGGGATTTGTTTTGAAGCGCaaaaa                                                                                                                 >  1:402228/1‑138 (MQ=255)
                    aaaaGCAGTTGATCTGGCAAGGTTTTGGCGCACTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGCGTTTTGTTTTGAAGGGCAAATa                                                                                                                 >  1:650930/1‑140 (MQ=255)
                            ttGATCTGGCAAGGTTTTGGCGCACTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGTATAGCGTTATTGGTTTTGAAGGGCAAATACCACTGgc                                                                                                         >  1:190614/1‑140 (MQ=255)
                                    gCAAGGTTTTGGCGCACTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGGATTCGAGGGCTGTTTTTAAGGGCGAAAAGCACCGGCGCTTtctc                                                                                                 >  2:442860/1‑140 (MQ=255)
                                    gCAAGGTTTTGGCGCACTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTtctc                                                                                                 <  2:488601/140‑1 (MQ=255)
                                              ggCGCACTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACt                                                                                       <  2:650033/140‑1 (MQ=255)
                                                    cTGGCTCCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGACTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGc                                                                                 <  1:442860/140‑1 (MQ=255)
                                                    cTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGc                                                                                 <  2:446828/140‑1 (MQ=255)
                                                    cTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGc                                                                                 <  1:426305/140‑1 (MQ=255)
                                                      ggCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCa                                                                               <  1:605994/140‑1 (MQ=255)
                                                             aGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTaga                                                                        <  1:444182/140‑1 (MQ=255)
                                                             aGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTaga                                                                        <  2:535014/140‑1 (MQ=255)
                                                                   aTGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTg                                                                  >  2:39405/1‑140 (MQ=255)
                                                                    tGGCATCGCCGTCGCAATAACTCGCCGATGTCCAGCCACCTTTGCATGTAGCGGGTGTGGCGGATAGCGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGc                                                                 <  2:190614/140‑1 (MQ=255)
                                                                    tGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGc                                                                 <  2:661422/140‑1 (MQ=255)
                                                                    tGGCATCGCCGTCGCAAAAACTCGCGGATGTACAGACTCTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGc                                                                 <  1:250156/140‑1 (MQ=255)
                                                                     ggCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCt                                                                <  1:447863/140‑1 (MQ=255)
                                                                       cATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGTTTTTAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGGAGAAGAGTGCTgg                                                              >  2:442279/1‑140 (MQ=255)
                                                                          cGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGAGAGGGGGTTTTGTTTTTAGGCGCAAAAAACACAGGCGCCTTCTCGCGCAGCGCTGCGGGTCAGAGGAGAAGAGGGCTGGcct                                                           >  2:510403/1‑138 (MQ=255)
                                                                                                    cAGACACTTTTGCATGAAGCGGTTGTGGCGGTTAGCGTGTTTTGTTTTTAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTACGa                                 <  2:222113/140‑1 (MQ=255)
                                                                                                        cacTTTTGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTACGAAAcc                             <  1:122043/140‑1 (MQ=255)
                                                                                                              tGCATGAAGCGGGTGTGGCGGATAGCGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTACGAAACCTTTAtt                       <  2:585755/140‑1 (MQ=255)
                                                                                                                     aGCGGGTGTGGCGGATAGCGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCCGAGTAGAAGAGTGCCGGCATTTAACTGAACAAAATGCCATGTACGCAACCTTTATTCAGTCAt                >  2:194192/1‑140 (MQ=255)
                                                                                                                             tgGCGGATAGCGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTACGAAACCTTTATTCAGTCATTccgcccg        <  1:75901/140‑1 (MQ=255)
                                                                                                                               gCGGATAGCGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAACTGCGATGAACGAAACCGTTATTTAGTCATTCCGCCCGCt      >  2:237505/1‑140 (MQ=255)
                                                                                                                                   aTAGCGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTACGAAACCTTTATTCAGTCATTCCGCCCGCTGGTg  <  1:466984/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                        
ATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTGGCAAGGTTTTGGCGCACTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTACGAAACCTTTATTCAGTCATTCCGCCCGCTGGTG  >  NC_000913/1463898‑1464168

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