Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,560,309 T→A P124P (CCA→CCT yjiH ← uncharacterized protein YjiH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,560,3090TA93.5% 77.0 / ‑0.5 31P124P (CCA→CCTyjiHuncharacterized protein YjiH
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base A (15/14);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-01

GCTCCGACAAATTTAAACACCAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGCC  >  NC_000913/4560064‑4560553
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
gCTCCGACAAATTTAAACACCAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAAt                                                                                                                                                                                                                                                 <  1:132547/251‑1 (MQ=255)
        aaaTTTAAACACCAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCcaaca                                                                                                                                                                                                                                         <  1:383385/251‑1 (MQ=255)
                  acCAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGtttt                                                                                                                                                                                                                                 >  1:295012/1‑249 (MQ=255)
                        aTGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAgg                                                                                                                                                                                                                         >  1:108465/1‑251 (MQ=255)
                         tGACGGCTAAAGGTAAGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGCGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGACAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTGAGCGATCAGCGCCAGCGAGCAGGTCCCAGGAATGCCAAACAACGGTTTTAGCACAgt                                                                                                                                                                                                                         >  2:371621/1‑249 (MQ=255)
                            cGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCTAAAATCACTTTAGCGCGTTCGGTAATTGCGCCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTc                                                                                                                                                                                                                     <  1:356915/251‑1 (MQ=255)
                                 aaaGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTGCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCAtt                                                                                                                                                                                                                  <  1:332220/249‑1 (MQ=255)
                                         gattacgggtgtgtcaaGAAAGGCAAACACCGCGACTCCGGAAGAGAAGTAGTTGGGGATGATGGCGCGACCGCTGGTTTGATTTGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCGTTGGTCATACCGGCGGCGGCAGCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGAGCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTgcg                                                                                                                                                                                                        <  2:298357/235‑1 (MQ=255)
                                           atACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTgcgcg                                                                                                                                                                                                      <  1:10544/250‑1 (MQ=255)
                                                     gCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATc                                                                                                                                                                                            >  1:40630/1‑251 (MQ=255)
                                                                  ccgagcgtgACGCAGGAAGAGAATTAGTTTGTGATGATGTCGTTACCGCGTGTTTGTTATGCTGCGAAAATCACTTTATCGCGTTCGGTAATTTCTCCTTCCTTCTCCATCTCGTTGGTCATACCGGCGGCGGCATCTGTGTTTTGCAGGTTAGCGATCAGCGCCCGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTATCTGCTGCGCGTCGCGTAATCCTCCAAGCCCATc                                                                                                                                                                               <  2:212786/243‑1 (MQ=255)
                                                                    cacCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCgg                                                                                                                                                                             >  1:264249/1‑251 (MQ=255)
                                                                         cgacgaCGGAAGAGAAGTATATGGTGATGATGGCGCTACCTCTGGATTGACATGCGTCGAAAATACCTTTATCGCGTTCGGTAATTTCGCCTTCCTGCGCCAGCTCTTTGATCATTCCCTCGGCGGCATCGGTGTTTTGCAGGTGTGCTAGCAGCGCCAGTGGGCATATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCTCGTAATCCACCCAGACCATCGGCAAtg                                                                                                                                                                        <  2:30144/251‑2 (MQ=255)
                                                                         cgACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATc                                                                                                                                                                               >  1:212786/1‑244 (MQ=255)
                                                                         cgACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAAtt                                                                                                                                                                        >  1:371299/1‑251 (MQ=255)
                                                                                  aaGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATc                                                                                                                                                               >  1:325754/1‑251 (MQ=255)
                                                                                               ggtgtgcatgGCGCTCGCGCTGGTTTGATATGCGGCGATAAGCACTTTATCGCGTGCGGTAATTTCTCCTTCCTTCGCCAGCTCTTTGGTCATACCGGCTGCTGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCTCAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTTATACAATCCCCagagagata                                                                                                                                                   <  2:325754/243‑3 (MQ=255)
                                                                                                            cTACCGCTGGTTTGATATGCGGCGCAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCCGCGCCAGCGGGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGACGCGAAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGc                                                                                                                                     >  2:67938/1‑251 (MQ=255)
                                                                                                                   tGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCATGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGCTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCAACCATGACCATCAGTAATTGAAATAATCACCAGAGAGAGAATGACTGACGGTGCCAGCTcc                                                                                                                              >  2:85703/1‑251 (MQ=255)
                                                                                                                                    aaaaTCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCACGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGcc                                                                                                             >  1:285001/1‑251 (MQ=255)
                                                                                                                                                          tAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGCGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCAccc                                                                                       >  1:264613/1‑251 (MQ=255)
                                                                                                                                                                ccccTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTGTGCAGGTTAGCGATCAGCGCCAGCGAGCAGAGCCCAGGAATGCCCAACAACGGTTTGAGCACAGGGGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCATCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACccgcgcc                                                                                  <  2:40630/248‑1 (MQ=255)
                                                                                                                                                                                            aCCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCTGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCAc                                                     >  1:171684/1‑251 (MQ=255)
                                                                                                                                                                                                gcggcggcATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAAAAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAgg                                                 <  1:165255/251‑1 (MQ=255)
                                                                                                                                                                                                    cggcATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGt                                             >  1:199217/1‑251 (MQ=255)
                                                                                                                                                                                                     ggcATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACATGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGtt                                            <  1:27119/251‑1 (MQ=255)
                                                                                                                                                                                                          cggggttgtGCAGGGTAGCGATCAGGGCCAGCGAGCAGAGCCCAGGAGTGCCCAACAGCGGGGTTAGCACAGGTGTCATTAACTGCGGCGCGGCGCGTAAGCCACCAAGACCATCGGTAATTGAAATAAGCCCCAGAGAGAGAATGACTGACGGTGCCAGCGCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCGCCAggg                                                <  2:82418/234‑2 (MQ=255)
                                                                                                                                                                                                          cGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGc                                                >  1:82418/1‑242 (MQ=255)
                                                                                                                                                                                                                   gCAGGTTAGCGATCAGCGCCAGCGAGCATATCCCAGGATGGCCCAACAACGGTGGTAGCACAGGTGTCATTAACTGCTGCTCGGCACGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGATAGAGATTGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCAGCTTTAGCCCTCGCGCCACCCGCGCCTCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAAAGATcc                               <  2:363741/250‑1 (MQ=255)
                                                                                                                                                                                                                                cAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGaa                  <  1:68194/250‑1 (MQ=255)
                                                                                                                                                                                                                                               gATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGcc  <  1:85703/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
GCTCCGACAAATTTAAACACCAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGCC  >  NC_000913/4560064‑4560553

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