Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,269,395 T→C 7.0% intergenic (‑383/‑19) yjbS ← / → aphA uncharacterized protein YjbS/acid phosphatase/phosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,269,3950TC 7.0% 145.8 / 2.2 58intergenic (‑383/‑19)yjbS/aphAuncharacterized protein YjbS/acid phosphatase/phosphotransferase
Reads supporting (aligned to +/- strand):  ref base T (20/33);  new base C (2/2);  total (22/36)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

AGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGC  >  NC_000913/4269255‑4269531
                                                                                                                                            |                                                                                                                                        
aGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGACGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAACGTTTCAACATACTGa                                                                                                                                        >  2:259224/1‑143 (MQ=255)
 ggAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATACAAAGTTTCAACATACTGAc                                                                                                                                       >  2:253326/1‑143 (MQ=255)
      gTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTAttt                                                                                                                                  <  2:740165/143‑1 (MQ=255)
       tATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCACCATACTGACTATTTa                                                                                                                                 >  1:599714/1‑143 (MQ=255)
       tATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTa                                                                                                                                 <  2:117823/143‑1 (MQ=255)
        agttagggggTGAAATATAATCATGCTGCTTGTGATCACAAAAGCCGCAGATAAGCCTTCCGTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTGTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAg                                                                                                                                <  1:227047/136‑1 (MQ=255)
        aTTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAg                                                                                                                                <  2:580863/143‑1 (MQ=255)
        aTTTACTGGGTTAAATATAATCACCCTGCTTTTCATCACAAAAACCGCATATAATCCTTCCTTTCCCCGGCAGCTGGCGTCATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAg                                                                                                                                <  1:146766/143‑1 (MQ=255)
         ttgccTGGCTTACATGTAACCCTCCTGCTTGTCATCACCAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACCAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAgg                                                                                                                               <  1:830490/139‑1 (MQ=255)
                    tattataAACATCCTGTTTTTCAGCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATgc                                                                                                                    <  1:901960/140‑1 (MQ=255)
                          aaTCCTCCTGCCTCTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAAGAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGa                                                                                                              <  2:351631/143‑1 (MQ=255)
                                   gCTCTTCATCACAAAATCCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATTGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAgg                                                                                                     <  2:565592/143‑1 (MQ=255)
                                    cTTGTCATCACAAAAAGCGCTGATAATCCTTGCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGc                                                                                                    <  2:547841/143‑1 (MQ=255)
                                     ttttCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCa                                                                                                   >  2:345658/1‑143 (MQ=255)
                                      tttCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCGCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCaa                                                                                                  <  1:129798/143‑1 (MQ=255)
                                       ttCACCACAAAACCCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAAt                                                                                                 <  1:552959/143‑1 (MQ=255)
                                           tcaCAAAAACCGCAGATAATCCTTCATTTCCCCGGCAGCTGGCGGTATGGTCAGATGGTTTTTGCAGCAACCCACACAATAACAAGTTTCAACATACTGACTATTTCGGGAGTAATCTGCGGAAGATGGCACGGGcagtcagt                                                                                             >  1:85307/1‑143 (MQ=255)
                                              cacaaaCCGCAGATAATCCTTCCTATACACGGCAGCTGGCGATAGGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGcc                                                                                          >  1:173301/4‑143 (MQ=255)
                                                    ccGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGc                                                                                    >  1:451415/1‑143 (MQ=255)
                                                         gATAATCCTTCCTTGCCCCGGCAGCGGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTAtt                                                                               <  2:212614/143‑1 (MQ=255)
                                                            aaTCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGtt                                                                            <  1:566532/143‑1 (MQ=255)
                                                             aTCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGACAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATCGTTc                                                                           >  1:325267/1‑143 (MQ=255)
                                                                   gcttCCCCGGGCAACTGGCGTTATGGGCAGATGGGTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATTCTTACTATTTAGTGCCGAATATGCGCAAGACCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTa                                                                     <  2:902237/142‑1 (MQ=255)
                                                                   ccttTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTa                                                                     <  2:709873/143‑1 (MQ=255)
                                                                       tCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACa                                                                 <  2:496227/143‑1 (MQ=255)
                                                                           cGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACCGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTc                                                             >  2:466423/1‑143 (MQ=255)
                                                                           cGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACCGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTc                                                             >  2:783719/1‑143 (MQ=255)
                                                                               aGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCt                                                         <  1:800926/143‑1 (MQ=255)
                                                                                     cGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTCTTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGcc                                                   >  1:155791/1‑143 (MQ=255)
                                                                                       ttATGGTCAGATGTTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCt                                                 <  1:299395/143‑1 (MQ=255)
                                                                                       atATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCt                                                 <  1:637941/142‑1 (MQ=255)
                                                                                        tATGTTTAGATGGTTTTTGCAACAAAGCACACAATAAAAAGTTCCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTg                                                <  1:350932/143‑1 (MQ=255)
                                                                                              cAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCa                                          <  1:536460/143‑1 (MQ=255)
                                                                                                gATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCAGTTTGATCAGCGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCAtc                                        >  2:541215/1‑143 (MQ=255)
                                                                                                     tttttGCAACGAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCtt                                   <  2:5270/143‑1 (MQ=255)
                                                                                                     tttttGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCtt                                   <  1:399259/143‑1 (MQ=255)
                                                                                                        ttGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTAGTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCAc                                <  2:849038/143‑1 (MQ=255)
                                                                                                         tGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTCAGCAGTTCCGCTGTTGCGCTTGTCGCCTGTCCCTCAcc                               >  2:487386/1‑143 (MQ=255)
                                                                                                            aacaaATCTCACAATAAAAAGTTTCAACATACTGACTGTTGAGTGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTCGCCTCATCTCCTTCACCGCt                            >  1:673731/1‑143 (MQ=255)
                                                                                                             acaaATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCtt                           >  1:105533/1‑143 (MQ=255)
                                                                                                                   ctcACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCt                     <  1:547162/143‑1 (MQ=255)
                                                                                                                    tcACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTg                    <  2:241439/143‑1 (MQ=255)
                                                                                                                     cacaATAAAAAGTTTCCACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTgg                   <  1:333436/143‑1 (MQ=255)
                                                                                                                     cacaATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGGCGTTGCCCTGGCCTCATCCCCTTCACCGCTTAACCCTgg                   >  2:443959/1‑143 (MQ=255)
                                                                                                                      acaATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTggg                  <  2:130982/143‑1 (MQ=255)
                                                                                                                      acaAGAAATAGGTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTggg                  <  2:583418/143‑1 (MQ=255)
                                                                                                                          tAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAAGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTa              >  1:475007/1‑143 (MQ=255)
                                                                                                                            aaaaGTTTCAACATACCGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAAc            <  1:783719/143‑1 (MQ=255)
                                                                                                                            aaaaGTTTCAACATACCGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAAc            <  1:466423/143‑1 (MQ=255)
                                                                                                                            aaaaGGTTCCACATACGGACTATGTAGGGAAAAATATGCGCAAGATCACACAGGCAAGCAGTGCCGTTTGCTTATTGTTCGCTCGAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAAc            <  1:92082/143‑1 (MQ=255)
                                                                                                                             aaaGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTCACg           >  1:269621/1‑143 (MQ=255)
                                                                                                                              aattttCCACATACTGACTATTTAGGGAAAACTATGCGCAAGATCACACCGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGTCTCATCTCCTTCTCCGCTTAACCCTGGGACTAACGt          <  2:251201/140‑1 (MQ=255)
                                                                                                                              aaGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTCACGt          >  1:409776/1‑143 (MQ=255)
                                                                                                                                  ttCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGCCTAACCTTGcc      >  2:365125/1‑143 (MQ=255)
                                                                                                                                   tCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCCTCTCCTTCACCGCTTAACCCCGGGACGACCGCTGCCa     >  2:678379/1‑143 (MQ=255)
                                                                                                                                     gaCATCCTGACTATTTAGGGAAAAATATGCGCAAGCGCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAgg   <  2:412323/142‑1 (MQ=255)
                                                                                                                                      aCATACTGACTATTTAGGGAAAACTATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGc  <  1:259224/143‑1 (MQ=255)
                                                                                                                                      aCATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTCACGTTGCCAGGc  >  2:72654/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                        
AGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGC  >  NC_000913/4269255‑4269531

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