Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,139 A→T 11.4% intergenic (+3/‑72) 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,1390AT11.4% 108.8 / 6.4 35intergenic (+3/‑72)yejM/proLessential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base A (15/16);  new base T (2/2);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-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                                                                                                    <  1:937048/126‑1 (MQ=255)
                                              tACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCa                                                                                                    >  2:937048/1‑126 (MQ=255)
                                                     tgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCg                                                                                                >  4:133142/1‑123 (MQ=255)
                                                     tgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCg                                                                                                <  3:133142/123‑1 (MQ=255)
                                                        aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGtt                                                                                                                        >  4:438599/1‑96 (MQ=38)
                                                        aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGtt                                                                                                                        <  3:438599/96‑1 (MQ=39)
                                                              gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                                                            >  1:641418/1‑86 (MQ=255)
                                                              gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                                                            <  2:641418/86‑1 (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)
                                                                           ccACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGatta                                                                                                                                      <  2:994403/63‑1 (MQ=255)
                                                                           ccACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGatta                                                                                                                                      >  1:994403/1‑63 (MQ=255)
                                                                                              tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                                                      >  3:1033199/1‑60 (MQ=255)
                                                                                              tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                                                      <  4:224482/60‑1 (MQ=255)
                                                                                              tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                                                      >  3:224482/1‑60 (MQ=255)
                                                                                              tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                                                      <  4:1033199/60‑1 (MQ=255)
                                                                                                  aGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATACATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                                                  <  1:780539/92‑1 (MQ=255)
                                                                                                  aGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                                                  >  2:780539/1‑92 (MQ=255)
                                                                                                        gACAGACGAGAAGCGTTTTATCGCTAACTGattaa                                                                                                                                     >  4:98671/1‑35 (MQ=255)
                                                                                                        gACAGACGAGAAGCGTTTTATCGCTAACTGattaa                                                                                                                                     <  3:98671/35‑1 (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                                        >  3:572926/1‑103 (MQ=39)
                                                                                                                                 aaCTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAg                                        <  4:572926/103‑1 (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)
                                                                                                                                         |                                                                                                                                    
AAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAAT  >  NC_000913/2286002‑2286271

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: ‑