Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,901,029 Δ1 bp 10.2% coding (557/1557 nt) yoaE ← putative membrane protein/conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,901,0290A.10.2% 156.8 / 12.2 49coding (557/1557 nt)yoaEputative membrane protein/conserved protein
Reads supporting (aligned to +/- strand):  ref base A (26/18);  new base . (2/3);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

AGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGAC  >  NC_000913/1900883‑1901174
                                                                                                                                                  |                                                                                                                                                  
aGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGcaccacc                                                                                                                                                  >  3:260058/1‑149 (MQ=255)
       aCGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgg                                                                                                                                           >  4:82733/1‑149 (MQ=255)
        cgcgCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgg‑g                                                                                                                                         <  1:421551/149‑1 (MQ=255)
             aaTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGCt                                                                                                                                    <  2:275671/149‑1 (MQ=38)
              aTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTg                                                                                                                                   <  4:260058/149‑1 (MQ=255)
              aTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGCtc                                                                                                                                   <  4:531722/149‑2 (MQ=255)
                    ttAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAAc                                                                                                                             <  2:241982/149‑1 (MQ=255)
                            ttCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATcgcg                                                                                                                     >  4:29473/1‑149 (MQ=255)
                            ttCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATAATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATcgcg                                                                                                                     >  1:323359/1‑149 (MQ=255)
                                 tgatgatCGAGAAGCCAATCGCGGCATACAGGTAACCTATCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGc                                                                                                                <  4:389813/149‑1 (MQ=255)
                                  gatgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGCtc                                                                                                                                   >  1:402653/1‑128 (MQ=255)
                                  gatgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGCtc                                                                                                                                   <  2:402653/129‑2 (MQ=37)
                                      atCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGTATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGttt                                                                                                           >  4:236496/1‑149 (MQ=255)
                                      atCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGttt                                                                                                           <  4:412843/149‑1 (MQ=255)
                                       tCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTg                                                                                                          >  2:495680/1‑149 (MQ=255)
                                         gagaAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGa                                                                                                        <  3:328835/149‑1 (MQ=255)
                                              gCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                   >  2:93123/1‑149 (MQ=255)
                                                 aaTCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCa                                                                                                >  1:250883/1‑149 (MQ=255)
                                                  aTCGCGGCATACAGGTAACCTTTCGGTATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGa                                                                                                        <  4:381737/140‑1 (MQ=255)
                                                  aTCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGa                                                                                                        >  3:381737/1‑140 (MQ=255)
                                                               agTAACCTTTCGGAATGTGGGAACCGAAACCTTCTGCCACCAGAATCAGACCACTCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCCGCAACATAACCGCGATc                                                                                  >  4:118632/2‑149 (MQ=39)
                                                                   aCCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCaa                                                                              <  1:57243/149‑1 (MQ=255)
                                                                      tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCa                                                                           >  4:503132/1‑149 (MQ=255)
                                                                           gAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTAcc                                                                      >  1:316163/1‑149 (MQ=255)
                                                                                  gaaaccgaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCca                                                               >  2:196409/1‑149 (MQ=255)
                                                                                       cgaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCTGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCTGCAACATAGCCGCCATCGCAATCACTACCGCCGCCatcatc                                                          >  4:350672/1‑149 (MQ=255)
                                                                                       cgaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAAcc                                                                                        >  1:364275/1‑119 (MQ=255)
                                                                                       cgaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAAcc                                                                                        <  2:364275/119‑1 (MQ=255)
                                                                                        gaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCAtcatca                                                         >  1:525036/1‑149 (MQ=255)
                                                                                         aaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCAc                                                        >  4:539283/1‑149 (MQ=255)
                                                                                              ttCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCa                                                   <  1:93123/149‑1 (MQ=255)
                                                                                              ttCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCa                                                   >  2:237833/1‑149 (MQ=255)
                                                                                                     accaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGtt                                            >  2:247380/1‑149 (MQ=255)
                                                                                                      ccaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTa                                           <  2:333628/149‑1 (MQ=255)
                                                                                                         gACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAAcc                                        <  4:246929/149‑1 (MQ=255)
                                                                                                            tCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATc                                     >  4:71483/1‑149 (MQ=255)
                                                                                                             cAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATcc                                    <  1:495680/149‑1 (MQ=255)
                                                                                                                  cAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGcc                                                                                                    >  1:79489/1‑80 (MQ=255)
                                                                                                                  cAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGcc                                                                                                    <  2:79489/80‑1 (MQ=255)
                                                                                                                           aCAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAAcc                                                                                        <  3:397078/83‑1 (MQ=255)
                                                                                                                           aCAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAAcc                                                                                        >  4:397078/1‑83 (MQ=255)
                                                                                                                           aCAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAtt                      <  2:127810/149‑1 (MQ=255)
                                                                                                                              ggAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTAcc                   >  1:330698/1‑149 (MQ=255)
                                                                                                                                    tCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTAc                    >  1:175548/1‑142 (MQ=255)
                                                                                                                                    tCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTAc                    <  2:175548/142‑1 (MQ=255)
                                                                                                                                     cagacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCa            >  3:247541/1‑149 (MQ=255)
                                                                                                                                      ttaaccagCACC‑CCACGGTGGTGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCaa           >  2:570520/7‑149 (MQ=255)
                                                                                                                                          agagCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACgag       >  3:149011/1‑149 (MQ=255)
                                                                                                                                          agagCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACgag       >  2:268304/1‑149 (MQ=255)
                                                                                                                                          agagCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCAACCAACgag       <  3:539283/149‑1 (MQ=255)
                                                                                                                                               accaccaccGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGAc  >  3:218548/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGAC  >  NC_000913/1900883‑1901174

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑