Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,407,293 T→A 6.6% intergenic (+27/‑82) accC → / → yhdT biotin carboxylase/DUF997 domain‑containing protein YhdT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,407,2930TA 6.6% 136.4 / 2.1 56intergenic (+27/‑82)accC/yhdTbiotin carboxylase/DUF997 domain‑containing protein YhdT
Reads supporting (aligned to +/- strand):  ref base T (29/23);  new base A (2/2);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

ACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTT  >  NC_000913/3407151‑3407435
                                                                                                                                              |                                                                                                                                              
aCGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAt                                                                                                                                                >  1:859720/1‑143 (MQ=255)
        aaaaCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAACACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGc                                                                                                                                        >  2:676526/1‑143 (MQ=255)
        aaaaCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAGGGGCGGGTTTGCCGGc                                                                                                                                        >  1:248054/1‑143 (MQ=255)
         aaaCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGcc                                                                                                                                       <  1:876612/143‑1 (MQ=255)
           aCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCtt                                                                                                                                     >  1:663537/1‑143 (MQ=255)
           aCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCtt                                                                                                                                     <  1:265860/143‑1 (MQ=255)
           aCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCAGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCtt                                                                                                                                     <  1:621174/143‑1 (MQ=255)
              aaCGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGACAAACAACTCGGTCTTCAGGAAAAGTAAGACTGCTAAAGCGGCAAAAGGCCGGATTGTCCGGCCtctgt                                                                                                                                  >  2:69321/1‑139 (MQ=255)
                cGTTGATCTGCAGATCCGGATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAATAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTtta                                                                                                                                <  2:612107/143‑1 (MQ=255)
                cGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCGACAAGCAGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTtta                                                                                                                                >  1:124911/1‑143 (MQ=255)
                   tGAGCTGAAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTttatta                                                                                                                             <  2:882940/143‑1 (MQ=255)
                        tGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTggg                                                                                                                        >  2:827587/1‑143 (MQ=255)
                               gcGCATCATGAACTACTAGAACTTCCAGCATGGTTGCACTAGCACGTATGAGTTTGTGACAGGTTTGGGGGTTCAGGAAAAATAAGACTGCTAAATCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGAc                                                                                                                 <  2:660536/142‑1 (MQ=255)
                                               gagaACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTACTACTGGGGATCGACAACCCCCATACGGCAc                                                                                                 >  2:635555/1‑143 (MQ=255)
                                                   aCTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATc                                                                                             >  1:89951/1‑143 (MQ=255)
                                                     ttCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATccc                                                                                           >  1:391195/1‑143 (MQ=255)
                                                     ttCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATccc                                                                                           >  2:240067/1‑143 (MQ=255)
                                                      tCCAGCATGGTGGCACTAACATACACTGTCTGGAGAAAAAGCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTGAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGGACAATcccc                                                                                          <  2:336659/143‑1 (MQ=255)
                                                      tCCAGCATGCTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATcccc                                                                                          <  1:526928/143‑1 (MQ=255)
                                                         aGCATGGTGGCACTAACATCCACTCTCTGGAGAAAACACTCGGTCTTCAGGAAAACTAAGACTGCTAAACCGTCAAAAGGCCGGCTTCTCCGGCCTTTTCTATTACTGGGGATCGACAACCCCCATCAGGTACAATCACCGCt                                                                                       >  2:76672/1‑143 (MQ=255)
                                                           cATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCttt                                                                                     <  2:174219/143‑1 (MQ=255)
                                                                tggCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCCAAAGGCCGGATTTTCCGGCCTTTCTTATTACTAAGGATCGACAACCCCCACACCGTACAATCCCCGCTCACTTCa                                                                                >  2:38957/1‑143 (MQ=255)
                                                                    aCTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCa                                                                            >  1:682058/1‑143 (MQ=255)
                                                                          aTCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCCTAAGGTCCAATCCCCGCTTTCTTCACCCATCAggg                                                                      >  2:244813/1‑143 (MQ=255)
                                                                          aTCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAggg                                                                      >  1:781779/1‑143 (MQ=255)
                                                                           tCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGa                                                                     <  2:859720/143‑1 (MQ=255)
                                                                           tCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCACATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGa                                                                     >  1:683025/1‑143 (MQ=255)
                                                                           tCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAATTTCCGGCCTTTTTTATTACTGGGGATCGCCAACCACCCTAAGGTACAATGCCCGCTTTCTTCACCCACCAGGGa                                                                     >  2:139549/1‑143 (MQ=255)
                                                                               cTATCTGGAGAAAAAACTCGGTCTTCAGAAAAACTAAGGATGCTGAAGTTGCAAGAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGCCaaa                                                                 >  1:372753/1‑143 (MQ=255)
                                                                                   cTGGAGAAAAAACTCGGTCTTCAGGAAAAAGAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAAt                                                             <  2:345505/143‑1 (MQ=255)
                                                                                    tGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAACTg                                                            >  1:138070/1‑143 (MQ=255)
                                                                                      gggaAAGAACTGGGTCTTCAGGAACACTAAGACTGCTGAACCGTCACAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTAGAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGa                                                          <  1:377236/141‑1 (MQ=255)
                                                                                         aaaaaaCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATCAGGTACAATCCCCGCTTTCTTCACCCATCACGGACAAAAAACGGaccc                                                       >  2:479760/1‑141 (MQ=255)
                                                                                            agacgcGGTCTTCAGGAACAAGAAGCCGGCTAAAGCGGCAAAAGGCCGGATGTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACAGTCg                                                    <  2:769753/138‑1 (MQ=255)
                                                                                              aCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGtt                                                  <  1:827587/143‑1 (MQ=255)
                                                                                                ttgttCTTCAGGAAAATTAAGACTGCTAAAGCGTCAAAAGGCCGGATGTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGtttt                                                <  1:267826/139‑1 (MQ=255)
                                                                                                tCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACACTCCCCGCTTTCTTCACCCATCAGGGACAAACCATGGACCCTCGtttt                                                >  1:621720/1‑143 (MQ=255)
                                                                                                tCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGtttt                                                >  1:463837/1‑143 (MQ=255)
                                                                                                 cGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAGAAAAAGGACACACGctttt                                               >  1:362126/1‑142 (MQ=255)
                                                                                                     cTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCa                                           <  1:560995/143‑1 (MQ=255)
                                                                                                      ttCAGGAAAACTAAGCCTGCTAAGGCGTCAAAAGGCCGCATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAg                                          <  2:875203/143‑1 (MQ=255)
                                                                                                      ctCAAGAAAAATAGGACTGCTAAAGCGTCAAAAGGCCGCATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCAGAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAg                                          <  2:608595/142‑1 (MQ=255)
                                                                                                       tCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAgg                                         <  1:185738/143‑1 (MQ=255)
                                                                                                         aGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCCTAAGGTCCAATCCCCGCTCTCTTCGCCCACCAGGGCCAACACACGGACCCCCGTTTCGTCCCGGcc                                       >  2:261317/1‑143 (MQ=255)
                                                                                                          ggAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGAACCTCGTTTTGTTCAGGccc                                      >  1:107458/1‑143 (MQ=255)
                                                                                                            aaaaaTAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCAt                                    >  1:861665/1‑143 (MQ=255)
                                                                                                              aaaTAAGACGGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATaa                                  <  1:300931/143‑1 (MQ=255)
                                                                                                                  aaGACTGCTCACGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAagag                              <  2:559418/143‑1 (MQ=255)
                                                                                                                     aCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGCCAAAAAATGGCCACTCGTTTTGTTCAGTCCCATAAACGGgcg                           >  2:258261/1‑143 (MQ=255)
                                                                                                                         ttAAAGCGTCAAATGGCCGGTTTTTCCGGTCTTTTTTATTACTGGGGATCGACAACCCCCATGAGGTACAATCCCCGCTTTCTTCCCCCATCAGGGACAAAGAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCgcgct                       <  1:783764/142‑1 (MQ=255)
                                                                                                                           aaaGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCgcgctgg                     >  1:96808/1‑143 (MQ=255)
                                                                                                                                  aaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAACAATGGACACTCGTTTTGTTCCGGCTCATGAAGAGGCgcgctgagcgctgg              >  1:319838/2‑143 (MQ=255)
                                                                                                                                   aaaaGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACACCCCCCAGAAGGTACAATCCCCGCTTTCTGCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGcgctgggcgctggg             <  2:598807/143‑1 (MQ=255)
                                                                                                                                        gCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGa        <  1:729566/143‑1 (MQ=255)
                                                                                                                                           ggATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACCCTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGAGCTCAcac     >  1:818647/1‑141 (MQ=255)
                                                                                                                                           ggATTTTCCGGCCTTTTGTATTAGTGGGGCTACATAAGGTCCAGGAGATACAATCCACGCATTCTTCAACCCTTAGGGACAAAAAATGGACACTCGTTTTGATCAGGCCCATAAATAGGCGCGATGGGCGCTGGGGCTgcacc     >  2:271281/1‑139 (MQ=255)
                                                                                                                                            gATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACTAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCt    <  1:676526/143‑1 (MQ=255)
                                                                                                                                              ttttCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCGGGGGCTGACCCttt  >  1:455750/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                              
ACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTT  >  NC_000913/3407151‑3407435

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