Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,987,442 Δ1 bp 9.8% intergenic (+14/+65) ftnB → / ← yecJ putative ferritin‑like protein/DUF2766 domain‑containing protein YecJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,987,4420T. 9.8% 143.1 / 11.4 51intergenic (+14/+65)ftnB/yecJputative ferritin‑like protein/DUF2766 domain‑containing protein YecJ
Reads supporting (aligned to +/- strand):  ref base T (24/22);  new base . (2/3);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.88e-01

GAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTG  >  NC_000913/1987300‑1987579
                                                                                                                                              |                                                                                                                                          
gAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATAGCACATCATCGGCGCt                                                                                                                                            >  1:106586/1‑143 (MQ=255)
 aaCAGCCGCATGATGGTCTCTTGCTGCAACCCACTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTa                                                                                                                                           <  1:596027/143‑1 (MQ=255)
     gcagcaTGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTACTGc                                                                                                                                       >  1:574414/1‑143 (MQ=255)
            gatgGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGACACGTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCAATGCAttagcgc                                                                                                                               <  1:232388/143‑6 (MQ=255)
                 tCTGTTGGTGCAAACCAGTCTTGATGAAGTGCGGAGTGCGATACTTGCGGGTATGTGCACTGTGAAGCCCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCAATGCATTAGCGCCGATg                                                                                                                          <  2:53822/143‑1 (MQ=255)
                  cggtggctgcAAACCATTCTTGATGTAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCAATGCATTAGCGCCGATGa                                                                                                                         <  2:519271/138‑1 (MQ=255)
                   tGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGaa                                                                                                                        <  2:490076/143‑1 (MQ=255)
                    gTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGACGAAg                                                                                                                       >  2:516794/1‑143 (MQ=255)
                         tgcAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGCAGGtttt                                                                                                                  >  1:733963/1‑143 (MQ=255)
                         tgcAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGtttt                                                                                                                  <  2:98066/143‑1 (MQ=255)
                           cAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTga                                                                                                                <  2:20411/143‑1 (MQ=255)
                           cAAACCATTCTTGATGAAGCGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTga                                                                                                                <  1:235728/143‑1 (MQ=255)
                                 aTTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAAcc                                                                                                          >  2:333902/1‑143 (MQ=255)
                                  ttCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCg                                                                                                         >  2:309045/1‑143 (MQ=255)
                                    cTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCt                                                                                                       <  2:410154/143‑1 (MQ=255)
                                           aaGTGTGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCAt                                                                                                <  1:198732/143‑1 (MQ=255)
                                                cgcAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACGGCGCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAAGCGCTGCCTCATCGGtt                                                                                           >  2:442274/1‑143 (MQ=255)
                                                  cAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGCGACCCCGCTGCCTCATCTGTTTg                                                                                         >  2:126354/1‑143 (MQ=255)
                                                   aGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGa                                                                                        <  1:352035/143‑1 (MQ=255)
                                                       cGAAACTTGCGGGTATGTGCCCTGTGCAGACAGACCAACATGTACTGAATGTCGTGTCACACCAGCTGGATGGATCATCATCGGCGCTGATGCTTTGCGTCGCTGAAAGGTCTGAGCAAACGCAACCTCATCTGTTTGAagca                                                                                    >  2:46175/1‑143 (MQ=255)
                                                           aCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCgg                                                                                <  1:593325/143‑1 (MQ=255)
                                                               gCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCCTTGATCAGCAGCGGCGCTAATGCATTGCGCCGATGAAGGGTTTGAGAAACCGCTGCCTCACCTGTTTGAAGCAGCGGtttt                                                                            >  2:744510/1‑143 (MQ=255)
                                                               gCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGtttt                                                                            <  1:651266/143‑1 (MQ=255)
                                                                cGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCGTCTGTTTGAAGCAGCGGttttt                                                                           >  2:616855/1‑143 (MQ=255)
                                                                     aTGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt                                                                      <  1:333902/143‑1 (MQ=255)
                                                                                  aGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCtg                                                         >  1:633797/1‑143 (MQ=255)
                                                                                    actgaccAACACGTTCTTGATGTCGTGTCACACCAGTTGCTGGGCTCATCATCGGCGCTAATGCATTGCTCCGATGAAGGTTTTGAGAAACCGCTGCCTCAGCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCtgtg                                                       <  2:320217/140‑1 (MQ=255)
                                                                                       gaccAACATGTTCTGAATGGCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTgggg                                                    <  1:83649/143‑1 (MQ=255)
                                                                                           aaCATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGAATTCCCCCTGTGGGGTaca                                                >  1:556982/1‑141 (MQ=255)
                                                                                                  tCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGGAGCAGCGGTTTTTTTAATGGGATTCCCCCTGTGGGGTAAGCGTGAGt                                         >  1:197561/1‑143 (MQ=255)
                                                                                                   cTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGtt                                        <  2:238198/143‑1 (MQ=255)
                                                                                                    tGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTc                                       <  1:233909/143‑1 (MQ=255)
                                                                                                     gAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCa                                      <  1:445924/143‑1 (MQ=255)
                                                                                                          tCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTGTTTTGAGGGGAATCACCCAGGGGGGGGAGCTGGAGTTGaaggag                                 >  1:256189/1‑139 (MQ=255)
                                                                                                           cGTGTCACACCAGCTGCATTGCTTATTATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAgc                                <  2:128221/143‑1 (MQ=255)
                                                                                                           cGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTGTTAATGGGATTCACCCTGTGGGGTAAACGCGAGGTCAATAAgc                                >  2:30248/1‑143 (MQ=255)
                                                                                                                cacaGCAGCTGCATTGATCATCATCGGCGCTAATGCATGGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGAt                           <  1:207196/143‑1 (MQ=255)
                                                                                                                  cacCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGGTTCACCCTGTGGGGTCAACTTGAGTTCAATGAGCGCGATgg                         >  1:132628/1‑143 (MQ=255)
                                                                                                                       gCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCtttt                    <  1:126354/143‑1 (MQ=255)
                                                                                                                       gCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCtttt                    >  1:794487/1‑143 (MQ=255)
                                                                                                                           cATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTCTATCATCGGCTTCTCCCGGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGa                <  1:766439/143‑1 (MQ=255)
                                                                                                                            aTTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCGGTGTGAAGCAGCGGGTTTTTTAATGGGATTCACCCTGTGGGGTAAACGTGAGTTCAATAAGCGCGATGGCTTTTGGGAt               >  2:643043/1‑143 (MQ=255)
                                                                                                                                atcatcatcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTGGAGTTCAATAAGCGCGCGGGCTTTTTGGAtgcca           >  2:84437/1‑139 (MQ=255)
                                                                                                                                 tcatcatcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGccc          >  2:454197/1‑143 (MQ=255)
                                                                                                                                  cttcatcGGCGCAATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCg         >  1:39625/3‑143 (MQ=255)
                                                                                                                                   atcatcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTAACTAAAGGAGCGGCTTTTTTATTGGGATTCACCCTGTGGGGCAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGc        >  1:373596/1‑143 (MQ=255)
                                                                                                                                     catcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCAt      <  1:406344/143‑1 (MQ=255)
                                                                                                                                      atcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGCAGCAGCGGTGTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAGGCGCGATGGCTTTTTGGCTTGCCCGCGTg     >  2:552731/1‑143 (MQ=255)
                                                                                                                                      atcGGCGCAATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATg     >  1:733310/1‑143 (MQ=255)
                                                                                                                                        cGGCGCTAATGCATTGCGCGGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGt   <  2:384307/143‑1 (MQ=255)
                                                                                                                                         ggCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTg  >  1:523831/1‑143 (MQ=255)
                                                                                                                                         ggCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTg  >  1:202125/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                          
GAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTG  >  NC_000913/1987300‑1987579

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: