Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,145,170 G→T 35.1% intergenic (+43/+22) yegJ → / ← yegK uncharacterized protein/ser/thr phosphatase‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,145,1700GT35.1% 29.3 / 10.9 51intergenic (+43/+22)yegJ/yegKuncharacterized protein/ser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (20/12);  new base T (10/8);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.77e-01

TATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCG  >  CP009273/2145023‑2145316
                                                                                                                                                     |                                                                                                                                                  
tataCGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTc                                                                                                                                                  <  2:143667/150‑1 (MQ=255)
                aaaTTTACGGCGGTTATACAATTGATCTTTTATTGGTCACCTATCAAAAAGAAGAAGCCGACGAGAAAAGAGAAAAATTAGTTCGCTAACCAGCCGTTATTTTTATTGTTAAACCGGGGGGTTGAATTGTTGTCTTCAACCGCCATgaga                                                                                                                                  <  2:140136/150‑1 (MQ=255)
                        ggcggTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTTTTGTCTTTCACCGCCCTGAGAAaatcact                                                                                                                          >  1:1025148/1‑150 (MQ=255)
                          cggTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTGTTTTCAACCGCCATGAGAAACTCActct                                                                                                                        >  2:1140426/1‑150 (MQ=255)
                          cggTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTCAATTGTTGGTTTTATCCCCCCTGAGAAACTCActct                                                                                                                        >  2:834388/1‑150 (MQ=255)
                           ggTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTg                                                                                                                       >  1:476866/1‑150 (MQ=255)
                            gTTATACCATTGATCTTTTATTGGTCACATATCAAAAAGAAAAAGCCAACGAGAAAAAAGAAAAATTAGTTCGAAAACAAGCAGTTATATTAATTTTTAAAACGGGCGGTTGATTTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGc                                                                                                                      <  1:203899/150‑1 (MQ=255)
                            gTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAACCGCCATGAGAAAATCACCCTGc                                                                                                                      >  1:856586/1‑150 (MQ=255)
                                 aCCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTTTTCAACCGCCATGAGAAAATCACTCTGCCcac                                                                                                                  >  2:839504/1‑149 (MQ=255)
                                   cATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAgc                                                                                                               >  2:173378/1‑150 (MQ=255)
                                   cATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTTTTCTACCCCCCTGAGAAAATCACTCTGccccccgc                                                                                                               >  2:346803/1‑145 (MQ=255)
                                      tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTGATTGTTGTCCTCCACCACCATGTGAAACACACTCTGCCCACAgcga                                                                                                             >  2:234404/1‑148 (MQ=255)
                                       gATCTTTTTTTGCTCCCCTATCAAAAAAAAGAAGCCGACGAGCAAAGAAAAAAATTAGTTCTCTAATCAGCCTTTATTTTTATTTTTAAAACGGGCGGTTGATTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGcc                                                                                                            <  2:624299/149‑1 (MQ=255)
                                               aTTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTTTTCAACCGCCCTTAGAAAATCACTCTGCCCCCAGCGCCCGTGCCAAc                                                                                                   >  2:1044760/1‑150 (MQ=255)
                                               aTTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTTAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGGGCCAAc                                                                                                   >  1:971076/1‑150 (MQ=255)
                                                ttgtcccccTAACAAAAAAAAGAGGCCGACGAGCTAAGAGAAAAATTATTTCGCTAATCGCCCTTTATTTTTATTGTTAAACCTGGCGGTTGAATTGTTTTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACg                                                                                                  <  2:731654/143‑1 (MQ=255)
                                                      aCCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTat                                                                                            >  1:1004128/1‑150 (MQ=255)
                                                       ccTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTat                                                                                            >  1:606401/1‑149 (MQ=255)
                                                       ccTATCAAAAAGAAGAAGCAGACGAGATAAGAGAAAAATTATTTCGATAATCGGCTTTTATATTTATTGTTAAACCGGGGGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatc                                                                                           <  2:693074/150‑1 (MQ=255)
                                                                 agaagaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTGTCTTCAACCGCCATGAGAAAATCCCTCTGCCCACAGCGCCAGTGCCAACGCTTTTTCATCCTCCGt                                                                                  >  1:872388/1‑149 (MQ=255)
                                                                 agaagaagCCGACGAGCAAAAAAAAAAATTAGTTCGAAAATCAGCTTTTATATTTATTGTTAAAACGGGGGGTTAAATTGTTGTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTa                                                                                 <  2:774071/150‑1 (MQ=255)
                                                                  gaggaggCCGACGAAAAAAGAGAAAAATTATTTCACTATTCAGCCTTTATATTTTTTGTTTAAACGGGGGGTTGATTTGTTGTTTTCAGCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTAc                                                                                <  2:33163/144‑1 (MQ=255)
                                                                   aagaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGGTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACg                                                                               >  2:14357/1‑150 (MQ=255)
                                                                    agaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCCTGAGAAAATCACTCTGCCCACAGCGCCCGTGCCAACGGTTTATCATCATCCGTACGc                                                                              >  1:446662/1‑150 (MQ=255)
                                                                      aagCCGAGGAGCAAAAAAAAAAATTATTTCGTTAATCACCGTTTATTTTTTTTGTTAAAACGGGCGTTTAAATTTTTGTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTc                                                                            <  2:881234/150‑1 (MQ=255)
                                                                      aagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAACGCCAGTGCCAACGGTTTATCATCATCCGTACGCTc                                                                            >  1:876550/1‑150 (MQ=255)
                                                                         ccGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGtt                                                                         >  2:213716/1‑150 (MQ=255)
                                                                         ccGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTCTCCTCAACCGCCATGAGAAAATCCCTCTCCCCACCGCCCCAATGCCAACGCCTTATCATTTTCCGTACCCTCgat                                                                         >  1:57248/1‑148 (MQ=255)
                                                                             cgaGCTAAGAGAAAATTTATTCCGTTAATCACCCTTTATTTTTATTTTTTAAACGGGGGTTTAATTTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTAcc                                                                     <  1:170678/150‑1 (MQ=255)
                                                                              gaGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCCACGTTTTTTCATCATCCCTACCCCCGTTTAccg                                                                    >  2:425972/1‑150 (MQ=255)
                                                                                gCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTTAACCGCCATTAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCCTACGCTCGTTTAccgcc                                                                  >  1:399957/1‑150 (MQ=255)
                                                                                  taaagaaaaaaTTTGTTTCGATATTCAGCTTTTATTTTTTTTTTTAAAACGGGGGTTTGATTTTTTGTCTCAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTAccgccg                                                                 <  1:253908/143‑1 (MQ=255)
                                                                                  tAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCgg                                                                >  2:400146/1‑150 (MQ=255)
                                                                                    agagCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCCTACGCCCGTTTACCGCCGGGc                                                              >  2:431716/1‑150 (MQ=255)
                                                                                     gagCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCt                                                             >  1:376525/1‑150 (MQ=255)
                                                                                       gCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAAACGCCATGAGAAACTCACTTTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTgg                                                           >  1:1000980/1‑150 (MQ=255)
                                                                                         aaaaTTATTTCGTTAATCCGCCTTTATTTTTTTTGTTAAACCTGGGGTTTAAATTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGAc                                                         <  1:716363/150‑1 (MQ=255)
                                                                                         aaaaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTTTTCCACCGCCCTGAGAAAATCACTCTGCCCACAGCGCCCGTGCCAACGTTTTTTCATCATCCGCACCCCCGTTTACCCCCCGGCCGGa                                                          >  1:487349/1‑149 (MQ=255)
                                                                                          aaaTTATTCCGCTACTCACCCTTTTTATTTTTTTTTAAACCGGGGGTTTGATTTTTTTCCTCCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACa                                                        <  1:997273/150‑1 (MQ=255)
                                                                                               aGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTGTCTTCAAACGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCTTGGCaaaaaa                                                   >  1:660163/1‑149 (MQ=255)
                                                                                                 tttCTCTAATCAGCCTTTATTTTTATTTTTAAAACGGGCGTTTAAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAAc                                                   <  2:164570/148‑1 (MQ=255)
                                                                                                      cTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTTAAACGCCATGAGAAACTCACTCTGCCCCCAGCGCCAGTGCCAAAGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCTGGACAAAAAAAGCttt                                             >  1:1190264/1‑150 (MQ=255)
                                                                                                        aaTCAGCTTTTTTTTTTTTTTTTAAAACGGGGGTTTAATTTTTTTTTTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTa                                            <  1:1044760/149‑1 (MQ=255)
                                                                                                        aaTCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTCTTCAACCCCCATTAGAAAATCACTCTGCCCCCCGCCCCCGTGCCCACGGTTTTTTATCATCCGTACGCCCGTTTTCCGCCCGGCTGGGCAAAAACTGCtttt                                            >  2:587407/1‑148 (MQ=255)
                                                                                                         aTCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATTATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCtttta                                           >  1:1103777/1‑147 (MQ=255)
                                                                                                                  tataTTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCCCCCGGCTGGACAAAAACTGCTTTTACAATTCAgg                                  >  2:294495/1‑149 (MQ=255)
                                                                                                                   attttttttttttAAACGGGCGGTTGATTTTTTGTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGt                                 <  2:476866/146‑1 (MQ=255)
                                                                                                                    tattttttTTTAAAACGGGGGTTTAATTTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGta                                <  2:446662/149‑1 (MQ=255)
                                                                                                                       ttattTTTTAAACGGGCGGTTAAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGtaata                             <  1:511822/149‑1 (MQ=255)
                                                                                                                              taaaaCGGGGGTTTGATTTTTTTTCTCCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAg                      <  2:561505/147‑1 (MQ=255)
                                                                                                                                 accggggggTTGATTTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgttgt                  <  1:1123577/143‑1 (MQ=255)
                                                                                                                                         ggtgtttttgttcTTTTCAGCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCg           <  1:346803/143‑1 (MQ=255)
                                                                                                                                                ttgttgTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtgctgc    >  2:17109/1‑150 (MQ=255)
                                                                                                                                                ttgttgTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtgctgc    >  1:1151904/1‑150 (MQ=255)
                                                                                                                                                  gttgtCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTgccg  >  1:1159955/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCG  >  CP009273/2145023‑2145316

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