Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,407,293 T→A 9.0% 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 9.0% 97.4 / 2.4 41intergenic (+27/‑82)accC/yhdTbiotin carboxylase/DUF997 domain‑containing protein YhdT
Reads supporting (aligned to +/- strand):  ref base T (13/24);  new base A (2/2);  total (15/26)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01

TCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGG  >  NC_000913/3407157‑3407424
                                                                                                                                        |                                                                                                                                   
tCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTCAAGCGTCAAAAGGCCGGATTTTCCg                                                                                                                               >  1:362231/1‑143 (MQ=255)
 cAAAACAAACGTTGATCTGGAGATCCGCATCGTGGGTGACGACAACATCCAGCATGGTGGCACTAACATCCACGATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCgg                                                                                                                              <  2:407190/143‑1 (MQ=255)
   aaaCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGcc                                                                                                                            >  2:141601/1‑143 (MQ=255)
      cgaaCGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCttt                                                                                                                         <  1:541349/141‑1 (MQ=255)
        aaCGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCttttg                                                                                                                       <  2:305543/143‑2 (MQ=255)
          cGTTGGTCGTCAGATCCGCAGCATGAATGACGAGAACTTCCGGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGTCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTtta                                                                                                                     <  2:41031/143‑1 (MQ=255)
          cGTTGATCTGCAGATCCGCATCATGAATCACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTtta                                                                                                                     <  1:146742/143‑1 (MQ=255)
           gTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAATTTCCGTCCTTTTttat                                                                                                                    >  2:95037/1‑143 (MQ=255)
                 cTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTgg                                                                                                              <  2:442954/143‑1 (MQ=255)
                   gCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTgggg                                                                                                            >  2:265875/1‑143 (MQ=255)
                    cAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAGCTCGGTCTTCAGGAAAAATAAGACTGCTAAGGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGa                                                                                                           <  2:124510/143‑1 (MQ=255)
                     cgctCCGCACCATGAATCACGAGAACTTCCAGCATGGGGGCACTAACATCCACTATCCGGAGAAAAAACGCGGTCTTCAGGAAAAAGAAGACTGCCAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGAt                                                                                                          <  1:411352/140‑1 (MQ=255)
                     aGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGAt                                                                                                          <  2:329446/143‑1 (MQ=255)
                           gCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACaa                                                                                                    <  1:334907/143‑1 (MQ=255)
                             atcagGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAAcc                                                                                                  <  2:470804/143‑1 (MQ=255)
                               catGAATGACGAGAACCTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAAcccc                                                                                                <  1:273622/143‑1 (MQ=255)
                                     tgaCGAGAACTTCCAGCATGGTGGCACTAACATCCAATATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTCCTGGGGGTCGACAACCCCCATacg                                                                                          >  1:375039/1‑141 (MQ=255)
                                              cTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGCAAAAACTCGGGCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAGAGGCCGGATCTTCCGGCCTTTTTTATCACTGAGGACAGACGACCCCACTAACGGACAATct                                                                                 >  2:371403/1‑142 (MQ=255)
                                              cTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACGGCTAAAGCGTCAAAAGGCCGGCTTTTGCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATcc                                                                                 >  1:253537/1‑143 (MQ=255)
                                                                        aCTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaa                                                       >  1:141322/1‑143 (MQ=255)
                                                                        aCTATCTGAAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaa                                                       <  1:440510/143‑1 (MQ=255)
                                                                              tGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATg                                                 <  1:142365/143‑1 (MQ=255)
                                                                                gagaAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGa                                               >  1:498933/1‑143 (MQ=255)
                                                                                        aCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTATTTATTACTGGTGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACGCGtt                                       >  1:364318/1‑143 (MQ=255)
                                                                                          tCGCTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGtttt                                     <  1:445522/143‑1 (MQ=255)
                                                                                            ggacTTCAGGAAAAATAAAAATGAAAATGCGTAAAAGATTCGGATGATCCGGCGAGTAGTATTACGGAAGATCGACAACGCCCATAATGTACAATCCCCGCTTTCTTTACAAAAAAAAAAAAAAAAATGGACACTCGTTTTGt                                   <  2:107474/140‑1 (MQ=255)
                                                                                               cTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGTGATCGACCACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCa                                <  2:427970/143‑1 (MQ=255)
                                                                                               cTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCa                                <  2:141322/143‑1 (MQ=255)
                                                                                               cTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCCGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCa                                <  1:27278/143‑1 (MQ=255)
                                                                                                 tCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTCCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATACGGTACAACCCCCGCTTTCTTCACCCATCAGGCACAAAAAATGGACACTCGTTTTGTTCAgg                              <  1:251917/143‑1 (MQ=255)
                                                                                                    ggAAAGATAAGACTGCTAAAGCCTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGccc                           <  2:355131/143‑1 (MQ=255)
                                                                                                         aaTAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATaaa                      >  1:548197/1‑143 (MQ=255)
                                                                                                            aaGACTGCTAAAGCGGCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAagag                   <  1:227286/143‑1 (MQ=255)
                                                                                                            aaGACTGCGAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAATTGTACACTCGTTTTGTTCAGGCCCATAAagag                   <  1:531404/143‑1 (MQ=255)
                                                                                                              gACTGCTAATCCGTCAAAAGGCCGGATTTGCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACGATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGgc                 <  2:540951/143‑1 (MQ=255)
                                                                                                               aCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCTCCGCTTTCGTCACCCATCAGGGACACAAAATGGACACACGTTTTGTTCAGCCCCATAAAGAGgcg                >  2:467407/1‑143 (MQ=255)
                                                                                                               aCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGgcg                <  2:498933/143‑1 (MQ=255)
                                                                                                                cTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGgcgc               >  2:113523/1‑143 (MQ=255)
                                                                                                                      aaGCGTGAAAGGGCCTGATTTTCCGGCCCGTTGTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCGTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCgcgctggg         <  2:276627/143‑1 (MQ=255)
                                                                                                                        gCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACCCGTTTTGTTCAGCCCCATAAAGAGGCgcgctgggcg       >  1:234938/1‑143 (MQ=255)
                                                                                                                             aaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACACAAAATGGACACTCGTTTTGTTCACGCCCATCAACACGTgcgctgggcgctgcg  >  1:257020/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                   
TCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGG  >  NC_000913/3407157‑3407424

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