Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,142 T→A 12.1% intergenic (+6/‑69) yejM → / → proL essential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,286,1420TA12.1% 99.7 / 6.7 33intergenic (+6/‑69)yejM/proLessential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base T (14/15);  new base A (2/2);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

AAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAAT  >  NC_000913/2286002‑2286271
                                                                                                                                            |                                                                                                                                 
aaGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                                                            <  4:704167/148‑1 (MQ=255)
   aCGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGGAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGt                                                                                                                         >  3:835892/1‑148 (MQ=255)
               aaCAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATc                                                                                                             <  4:835892/148‑1 (MQ=255)
                                   cACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                                                <  1:660840/109‑1 (MQ=255)
                                   cACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                                                >  2:660840/1‑109 (MQ=255)
                                              tACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCa                                                                                                    >  2:937048/1‑126 (MQ=255)
                                              tACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCa                                                                                                    <  1:937048/126‑1 (MQ=255)
                                                     tgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCg                                                                                                <  3:133142/123‑1 (MQ=255)
                                                     tgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCg                                                                                                >  4:133142/1‑123 (MQ=255)
                                                        aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGtt                                                                                                                        >  4:438599/1‑96 (MQ=38)
                                                        aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGtt                                                                                                                        <  3:438599/96‑1 (MQ=39)
                                                              gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                                                            <  2:641418/86‑1 (MQ=255)
                                                              gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                                                            >  1:641418/1‑86 (MQ=255)
                                                                  gaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGtt                                                                                                                        >  4:521689/1‑86 (MQ=255)
                                                                  gaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGtt                                                                                                                        <  3:521689/86‑1 (MQ=255)
                                                                        aaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaaa                                                                                                                              >  4:224132/1‑74 (MQ=255)
                                                                        aaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaaa                                                                                                                              <  3:224132/74‑1 (MQ=255)
                                                                                              tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                                                      <  4:224482/60‑1 (MQ=255)
                                                                                              tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                                                      <  4:1033199/60‑1 (MQ=255)
                                                                                              tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                                                      >  3:224482/1‑60 (MQ=255)
                                                                                              tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                                                      >  3:1033199/1‑60 (MQ=255)
                                                                                                  aGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATACATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                                                  <  1:780539/92‑1 (MQ=255)
                                                                                                  aGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                                                  >  2:780539/1‑92 (MQ=255)
                                                                                                            gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACg                                                                                    <  1:715576/80‑1 (MQ=255)
                                                                                                            gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACg                                                                                    >  2:715576/1‑80 (MQ=255)
                                                                                                                   aGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGa                                                                 <  2:435757/92‑1 (MQ=255)
                                                                                                                   aGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGa                                                                 >  1:435757/1‑92 (MQ=255)
                                                                                                                               cTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                                                  >  3:665135/1‑63 (MQ=255)
                                                                                                                               cTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                                                  <  4:665135/63‑1 (MQ=255)
                                                                                                                               cTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCGTGGGGTGTCGGGGGTCGGAGGTTCAAAt  >  1:869534/1‑143 (MQ=255)
                                                                                                                               cTAACTGATTAATTATAAATCAGTTAGCGAAATATCTGACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCGTGGGGTGTCGGGGGTCGGAGGTTCAAAt  <  2:869534/143‑1 (MQ=255)
                                                                                                                                 aaCTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAg                                        <  4:572926/103‑1 (MQ=39)
                                                                                                                                 aaCTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAg                                        >  3:572926/1‑103 (MQ=39)
                                                                                                                                 aaCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCg             >  4:121894/1‑130 (MQ=255)
                                                                                                                                 aaCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCg             <  3:121894/130‑1 (MQ=255)
                                                                                                                                        taattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt                                             >  1:916069/1‑91 (MQ=255)
                                                                                                                                        taattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt                                             <  2:916069/91‑1 (MQ=255)
                                                                                                                                            tataAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGt                                                               >  1:226020/1‑69 (MQ=255)
                                                                                                                                            tataAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGt                                                               <  2:226020/69‑1 (MQ=255)
                                                                                                                                            tataAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGc                                               <  2:406350/85‑1 (MQ=255)
                                                                                                                                            tataAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGc                                               >  1:406350/1‑85 (MQ=255)
                                                                                                                                            |                                                                                                                                 
AAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAAT  >  NC_000913/2286002‑2286271

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