Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,269,058 G→A intergenic (‑46/‑356) 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,0580GA97.1% 114.8 / NA 35intergenic (‑46/‑356)yjbS/aphAuncharacterized protein YjbS/acid phosphatase/phosphotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (12/22);  minor base C (1/0);  total (13/22)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGGTAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGCCGAGGATTCATCATAATAAACGTAAAATTAATGTATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCA  >  NC_000913/4268817‑4269300
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
aTAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCt                                                                                                                                                                                                                                           <  1:112133/251‑1 (MQ=255)
   aaTGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGt                                                                                                                                                                                                                                        <  2:196596/251‑1 (MQ=255)
       cGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAAt                                                                                                                                                                                                                                    <  1:229842/251‑1 (MQ=255)
       cGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAAt                                                                                                                                                                                                                                    <  2:208226/251‑1 (MQ=255)
         tGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATc                                                                                                                                                                                                                                   >  2:317351/1‑250 (MQ=255)
                    aCTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGCGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCtt                                                                                                                                                                                                                       >  1:589625/1‑251 (MQ=255)
                       tGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTtat                                                                                                                                                                                                                    >  2:566103/1‑251 (MQ=255)
                                    gATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGt                                                                                                                                                                                                       <  2:515825/251‑1 (MQ=255)
                                                     tCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATc                                                                                                                                                                                      >  2:105043/1‑251 (MQ=255)
                                                                 tttGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTTTTTTTGTAATGCTGTATCATTAAGTTc                                                                                                                                                                             <  1:256183/248‑1 (MQ=255)
                                                                 tttGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTc                                                                                                                                                                             <  1:117999/248‑1 (MQ=255)
                                                                 tttGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTc                                                                                                                                                                             <  2:507340/248‑1 (MQ=255)
                                                                            tAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTAc                                                                                                                                                                 >  1:441634/1‑249 (MQ=255)
                                                                                                           ttGGCATAAGAGATATCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCt                                                                                                                                    <  1:247086/247‑1 (MQ=255)
                                                                                                             ggCATAAGAGATATCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTg                                                                                                                              >  1:341040/1‑251 (MQ=255)
                                                                                                               cATAAGAGATATCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGcc                                                                                                                            <  2:68929/251‑1 (MQ=255)
                                                                                                                       atatCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGCCGAGGATtc                                                                                                                    <  1:250294/251‑1 (MQ=255)
                                                                                                                         atCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGCCGAGGATtcat                                                                                                                  >  2:195609/1‑251 (MQ=255)
                                                                                                                              tattatTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGCCGAGGATTCATCata                                                                                                              >  1:140111/1‑250 (MQ=255)
                                                                                                                               attattTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCCCCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGCCGAGGATTCATCataat                                                                                                            <  1:428295/251‑1 (MQ=255)
                                                                                                                                                  tCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGCCGAGGATTCATCATAATAAACGTAAAATTAATGTAt                                                                                         >  1:10460/1‑251 (MQ=255)
                                                                                                                                                   cAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGCCGAGGATTCATCATAATAAACGTAAAATTAATGTa                                                                                          <  1:82242/249‑1 (MQ=255)
                                                                                                                                                    aaCAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGCCGAGGATTCATCATAATAAACGTAAAATTAATGTATcc                                                                                       <  2:441634/251‑1 (MQ=255)
                                                                                                                                                                                 tCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGCCGAGGATTCATCATAATAAACGTAAAATTAATGTATCCTTACATCGAGTAATAAACATTTTTtata                                                           <  2:175477/250‑1 (MQ=255)
                                                                                                                                                                                   acaTATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGCCGAGGATTCATCATAATAAACGTAAAATTAATGTATCCTTACATCGAGTAATAAACATTTTTTATACaa                                                        <  2:250345/251‑1 (MQ=255)
                                                                                                                                                                                   acaTATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGCCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTAATTTGCCGAGGATTCATCATAATAAACGTAAAATTAATGTATCCTTACATCGAGTAATAAACATTTTTTATACa                                                         >  1:398166/1‑250 (MQ=255)
                                                                                                                                                                                       atGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGCTAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGCCGAGGATTCATCATAATAAACGTAAACTTAATGTATCCTTACATCGAGTAATAAACATTTTTTATACaaaaa                                                     >  1:415830/1‑250 (MQ=255)
                                                                                                                                                                                       atGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGCCGAGGATTCATCATAATAAACGTAAAATTAATGTATCCTTACATCGAGTAATAAACATTTTTTATACaaaaa                                                     <  1:459075/250‑1 (MQ=255)
                                                                                                                                                                                          atatatTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGCCGAGGATTCATCATAATAAACGTAAAATTAATGTATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGa                                                 <  2:348033/251‑1 (MQ=255)
                                                                                                                                                                                                  atatCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTc                                                                                                                                                                             <  2:632539/119‑1 (MQ=255)
                                                                                                                                                                                                  atatCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTc                                                                                                                                                                             >  1:632539/1‑119 (MQ=255)
                                                                                                                                                                                                                  cgcgAGTAATAAAATTAATCACAATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGCCGAGGATTCATCATAATAAACGTAAAATTAATGTATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTaa                          >  2:265310/1‑250 (MQ=255)
                                                                                                                                                                                                                               aTTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGCCGAGGATTCATCATAATAAACGTAAAATTAATGTATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTg            <  2:180448/251‑1 (MQ=255)
                                                                                                                                                                                                                               aTTAATCACCATTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGCCGAGGATTCATCATAATAAACGTAAAATTAATGTATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTg            <  1:541840/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aTTGTAGGATAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGCCGAGGATTCATCATAATAAACGTAAAATTAATGTATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTtcatca  <  2:629608/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
ATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCTTTATTATTTAGTGCATCTAGCTCAACAGGTAAAGATGTCGTTTTGTTGCTATTCACATATGATATATTCATATCATCATCTCCATCGCGAGTAATAAAATTAATCACCATTGTAGGGTAGGGGGCTGGTCAATCAGAAATCATCTTTATAAACTTCGATTGTTTTTGTAATGCTGTATCATTAAGTTCATTAAATCGTACAGCAGATAAATGTTCTATCAAATTTCGCTCATTTGCCGAGGATTCATCATAATAAACGTAAAATTAATGTATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCA  >  NC_000913/4268817‑4269300

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