Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,800,475 G→A *110* (TAG→TAA ygaM → putative membrane‑anchored DUF883 family ribosome‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,800,4750GA98.4% 203.0 / ‑4.5 64*110* (TAG→TAAygaMputative membrane‑anchored DUF883 family ribosome‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (31/32);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTT  >  NC_000913/2800341‑2800604
                                                                                                                                      |                                                                                                                                 
cgcgTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCa                                                                                                                               <  2:968885/139‑1 (MQ=255)
                              ggCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGAccc                                                                                                 <  1:714815/139‑1 (MQ=255)
                                          tCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCt                                                                                     <  1:1406688/139‑1 (MQ=255)
                                                    gTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATAcgc                                                                                                            >  2:164221/1‑106 (MQ=255)
                                                    gTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATAcgc                                                                                                            <  1:164221/106‑1 (MQ=255)
                                                      gAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCtt                                                                         >  1:1271276/1‑139 (MQ=255)
                                                               tGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTgg                                                                >  1:1254447/1‑139 (MQ=255)
                                                                ggTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCttt                                                                        <  2:650079/130‑1 (MQ=255)
                                                                ggTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCttt                                                                        >  1:650079/1‑130 (MQ=255)
                                                                      gCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATca                                                         >  1:36913/1‑139 (MQ=255)
                                                                         tGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACaa                                                      >  2:595707/1‑139 (MQ=255)
                                                                           ggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATAcgc                                                                                                            <  2:721171/83‑1 (MQ=255)
                                                                           ggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATAcgc                                                                                                            >  1:721171/1‑83 (MQ=255)
                                                                           ggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACaaaa                                                    >  2:1023596/1‑139 (MQ=255)
                                                                                agcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATccc                                               >  2:585030/1‑139 (MQ=255)
                                                                                     cTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCaaaaa                                                                                                                           >  2:716706/1‑58 (MQ=255)
                                                                                     cTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCaaaaa                                                                                                                           <  1:716706/58‑1 (MQ=255)
                                                                                       gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAccc                                        >  1:212054/1‑139 (MQ=255)
                                                                                           taggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTc                                    <  1:1074150/139‑1 (MQ=255)
                                                                                            aggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                   <  1:1023596/139‑1 (MQ=255)
                                                                                                aTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCa                                                                                                                               <  1:253669/43‑1 (MQ=255)
                                                                                                aTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCa                                                                                                                               >  2:253669/1‑43 (MQ=255)
                                                                                                aTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCt                               <  2:1254447/139‑1 (MQ=255)
                                                                                                  ttttATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGt                                                                              <  1:594062/90‑1 (MQ=255)
                                                                                                  ttttATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGt                                                                              >  2:594062/1‑90 (MQ=255)
                                                                                                       tCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCa                                                                                                                               <  1:510880/36‑1 (MQ=255)
                                                                                                       tCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCa                                                                                                                               >  2:510880/1‑36 (MQ=255)
                                                                                                         ggCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCt                                                                                         <  2:1480635/72‑1 (MQ=255)
                                                                                                         ggCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAg                      >  2:1102953/1‑139 (MQ=255)
                                                                                                         agcgcACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCt                                                                                         >  1:1480635/2‑72 (MQ=255)
                                                                                                                     aGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                   <  1:1612725/114‑1 (MQ=255)
                                                                                                                     aGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                   >  2:1612725/1‑114 (MQ=255)
                                                                                                                          gcgcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTa     >  2:905887/1‑139 (MQ=255)
                                                                                                                          gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATcccc                                              >  1:561696/1‑98 (MQ=255)
                                                                                                                          gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATcccc                                              <  2:561696/98‑1 (MQ=255)
                                                                                                                          gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATc                                                 <  2:1417921/95‑1 (MQ=255)
                                                                                                                          gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATc                                                 >  1:1417921/1‑95 (MQ=255)
                                                                                                                          gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                   >  1:649652/1‑109 (MQ=255)
                                                                                                                          gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                   >  2:1378911/1‑109 (MQ=255)
                                                                                                                          gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                   <  2:649652/109‑1 (MQ=255)
                                                                                                                          gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                   <  1:1378911/109‑1 (MQ=255)
                                                                                                                          gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACt      >  2:1543876/1‑138 (MQ=255)
                                                                                                                          gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACt      <  1:1543876/138‑1 (MQ=255)
                                                                                                                          gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTa     >  1:235962/1‑139 (MQ=255)
                                                                                                                          gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAACCCCCCTACCCCGTCACGCTCATATCCAGGGTAAtt               <  2:1370648/129‑1 (MQ=255)
                                                                                                                          gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAACCCCCCTACCCCGTCACGCTCATATCCAGGGTAAtt               >  1:1370648/1‑129 (MQ=255)
                                                                                                                           cgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACaaaaa                                                   >  2:1179130/1‑92 (MQ=255)
                                                                                                                           cgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACaaaaa                                                   <  1:1179130/92‑1 (MQ=255)
                                                                                                                           cgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCt                               <  2:633535/112‑1 (MQ=255)
                                                                                                                           cgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCt                               >  1:633535/1‑112 (MQ=255)
                                                                                                                           cgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt    <  1:1019432/139‑1 (MQ=255)
                                                                                                                           cgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt    <  2:671279/139‑1 (MQ=255)
                                                                                                                           cgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt    <  2:1271276/139‑1 (MQ=255)
                                                                                                                            gcAAATCGTAATGCAAAAGTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAgg                     >  1:897212/1‑121 (MQ=255)
                                                                                                                            gcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTg                                                                                                     <  2:545620/41‑1 (MQ=255)
                                                                                                                            gcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTg                                                                                                     >  1:545620/1‑41 (MQ=255)
                                                                                                                            gcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGt                                     <  1:1186592/105‑1 (MQ=255)
                                                                                                                            gcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGt                                     >  2:1186592/1‑105 (MQ=255)
                                                                                                                            gcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACg                                 <  1:796158/109‑1 (MQ=255)
                                                                                                                            gcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACg                                 >  2:796158/1‑109 (MQ=255)
                                                                                                                            gcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAgg                     <  2:897212/121‑1 (MQ=255)
                                                                                                                            gcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt    <  1:828284/138‑1 (MQ=255)
                                                                                                                            gcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt    >  2:828284/1‑138 (MQ=255)
                                                                                                                             cAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAttt  <  2:93255/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                 
CGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTT  >  NC_000913/2800341‑2800604

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