Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,165,628 A→T 7.5% intergenic (+107/‑40) yegP → / → yegQ UPF0339 family protein/putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,165,6280AT 7.5% 177.0 / 3.5 52intergenic (+107/‑40)yegP/yegQUPF0339 family protein/putative peptidase
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base T (2/2);  total (26/26)
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.97e-01

AGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTA  >  NC_000913/2165482‑2165768
                                                                                                                                                  |                                                                                                                                            
aGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAg                                                                                                                                             >  1:275298/1‑148 (MQ=255)
            aaGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAgg                                                                                                                                            >  4:304836/1‑137 (MQ=255)
            aaGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAgg                                                                                                                                            <  3:304836/137‑1 (MQ=255)
              gCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTg                                                                                                                               >  4:816583/1‑148 (MQ=255)
               ccAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCgg                                                                                                                                    >  4:1137812/1‑142 (MQ=255)
               ccAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCgg                                                                                                                                    <  3:1137812/142‑1 (MQ=255)
               ccAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGa                                                                                                                              <  2:275298/148‑1 (MQ=255)
                 aGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGAcc                                                                                                                            <  3:816583/148‑1 (MQ=255)
                 aGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGAcc                                                                                                                            <  3:1035753/148‑1 (MQ=255)
                            gACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTg                                                                                                                          >  4:487054/1‑139 (MQ=255)
                            gACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTg                                                                                                                          <  3:487054/139‑1 (MQ=255)
                                      aaTTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGa                                                                                                                                         >  1:1091808/1‑114 (MQ=255)
                                      aaTTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGa                                                                                                                                         <  2:1091808/114‑1 (MQ=255)
                                             aTGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTtcatca                                                                                                                                              <  3:168038/102‑1 (MQ=255)
                                             aTGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTtcatca                                                                                                                                              >  4:168038/1‑102 (MQ=255)
                                             aTGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTa                                                                                                >  2:517589/1‑148 (MQ=255)
                                                          gttgttGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCgg                                                                                                                                    >  1:778351/1‑99 (MQ=255)
                                                          gttgttGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCgg                                                                                                                                    <  2:778351/99‑1 (MQ=255)
                                                           ttgttgGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCtt                                                                                  <  1:517589/148‑1 (MQ=255)
                                                            tgttgGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGa                                                                                                                  <  1:261551/115‑1 (MQ=255)
                                                            tgttgGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGa                                                                                                                  >  2:261551/1‑115 (MQ=255)
                                                                           cTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGc                                                                  >  1:584195/1‑148 (MQ=255)
                                                                             ttACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACgagaga                                                                                                          >  4:501832/1‑106 (MQ=255)
                                                                             ttACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACgagaga                                                                                                          <  3:501832/106‑1 (MQ=255)
                                                                                      ccccccACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATg                                                       <  2:584195/148‑1 (MQ=255)
                                                                                        ccccACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTaa                                                                                               <  2:770763/106‑1 (MQ=255)
                                                                                        ccccACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTaa                                                                                               >  1:770763/1‑106 (MQ=255)
                                                                                         cccACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGa                                                                                                                                   <  2:457371/69‑1 (MQ=255)
                                                                                         cccACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGa                                                                                                                                   >  1:457371/1‑69 (MQ=255)
                                                                                           cacaGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTaa                                                                                               <  3:62951/103‑1 (MQ=255)
                                                                                           cacaGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTaa                                                                                               >  4:62951/1‑103 (MQ=255)
                                                                                                    aCAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAAcc                                                                                            <  2:655571/97‑1 (MQ=255)
                                                                                                    aCAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAAcc                                                                                            >  1:655571/1‑97 (MQ=255)
                                                                                                                  tAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATg                                    >  2:691527/1‑139 (MQ=255)
                                                                                                                  tAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATg                                    <  1:691527/139‑1 (MQ=255)
                                                                                                                                   ccTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGtcat                                                                                                                      <  2:790333/40‑1 (MQ=255)
                                                                                                                                   ccTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGtcat                                                                                                                      >  1:790333/1‑40 (MQ=255)
                                                                                                                                    cTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCttt                                                                                 <  3:1019987/76‑1 (MQ=255)
                                                                                                                                    cTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCttt                                                                                 >  4:1019987/1‑76 (MQ=255)
                                                                                                                                    cTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTcc                                                                               <  3:595037/78‑1 (MQ=255)
                                                                                                                                    cTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTcc                                                                               >  4:595037/1‑78 (MQ=255)
                                                                                                                                     tAACCGCTTCATCAGGTGAAGCGGATCTGACCTGt                                                                                                                         <  2:440217/35‑1 (MQ=255)
                                                                                                                                     tAACCGCTTCATCAGGTGAAGCGGATCTGACCTGt                                                                                                                         >  1:440217/1‑35 (MQ=255)
                                                                                                                                       aCCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTcc                                                                               <  2:1019400/75‑1 (MQ=255)
                                                                                                                                       aCCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTcc                                                                               >  1:1019400/1‑75 (MQ=255)
                                                                                                                                       aCCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCcgc      >  2:204086/1‑148 (MQ=255)
                                                                                                                                        ccGCTTCACCTGATGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGaa                                                                                        >  3:591275/1‑65 (MQ=38)
                                                                                                                                        ccGCTTCACCTGATGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGaa                                                                                        <  4:591275/65‑1 (MQ=38)
                                                                                                                                        ccGCTTCACCTGATGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGa                             <  2:257951/124‑1 (MQ=255)
                                                                                                                                        ccGCTTCACCTGATGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGa                             >  1:257951/1‑124 (MQ=255)
                                                                                                                                            ttCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTAt                                     >  2:751332/1‑112 (MQ=255)
                                                                                                                                            ttCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTAt                                     <  1:751332/112‑1 (MQ=255)
                                                                                                                                            ttCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTa  >  2:159876/1‑147 (MQ=255)
                                                                                                                                            ttCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTa  <  1:159876/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
AGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTA  >  NC_000913/2165482‑2165768

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