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,4750GA96.3% 159.9 / ‑1.6 54*110* (TAG→TAAygaMputative membrane‑anchored DUF883 family ribosome‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (26/26);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

CCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGC  >  NC_000913/2800346‑2800606
                                                                                                                                 |                                                                                                                                    
ccAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAAt                                                                                                                             <  1:420791/139‑1 (MQ=255)
                    ccGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGCTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATACTAAATACGCGTCttt                                                                                                        <  1:1198500/139‑1 (MQ=255)
                                                                aGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATc                                                            >  2:1584019/1‑139 (MQ=255)
                                                                 gCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTgg                                                                  >  2:129242/1‑132 (MQ=255)
                                                                 gCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTgg                                                                  <  1:129242/132‑1 (MQ=255)
                                                                  cGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcac                                                          >  1:1625843/1‑139 (MQ=255)
                                                                  cGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcac                                                          >  2:1283165/1‑139 (MQ=255)
                                                                      ggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGata                                                                                                                         >  1:1012303/1‑73 (MQ=255)
                                                                      ggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGata                                                                                                                         <  2:1012303/73‑1 (MQ=255)
                                                                          cagcagCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATcc                                                  <  2:66774/139‑1 (MQ=255)
                                                                           agcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAg                                                                                              >  1:1225046/1‑94 (MQ=255)
                                                                           agcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAg                                                                                              <  2:1225046/94‑1 (MQ=255)
                                                                           agcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATccc                                                 >  1:877215/1‑139 (MQ=255)
                                                                                 tGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAcc                                           >  2:1193053/1‑139 (MQ=255)
                                                                                  gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGGaaa                                                                                                                                <  2:1221252/54‑1 (MQ=39)
                                                                                  gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCaaa                                                                                                                                >  1:1221252/1‑54 (MQ=255)
                                                                                       aggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                     <  1:1554934/139‑1 (MQ=255)
                                                                                         gtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACg                                   <  2:901261/139‑1 (MQ=255)
                                                                                             ttttATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCa                               >  1:1750458/1‑139 (MQ=255)
                                                                                                  tCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCaaaa                                                                                                                               >  2:725597/1‑39 (MQ=255)
                                                                                                  tCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCaaaa                                                                                                                               <  1:725597/39‑1 (MQ=255)
                                                                                                    ggCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCtt                                                                           <  1:212260/88‑1 (MQ=255)
                                                                                                    ggCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCtt                                                                           >  2:212260/1‑88 (MQ=255)
                                                                                                    ggCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCt                                 <  1:106246/130‑1 (MQ=255)
                                                                                                    ggCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCt                                 >  2:106246/1‑130 (MQ=255)
                                                                                                    ggCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAg                        >  2:653865/1‑139 (MQ=255)
                                                                                                      cgcACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                     <  1:463606/124‑1 (MQ=255)
                                                                                                      cgcACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                     >  2:463606/1‑124 (MQ=255)
                                                                                                         aCTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTaa                   >  1:1152900/1‑139 (MQ=255)
                                                                                                             ttAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATcc                                                  <  1:72712/104‑1 (MQ=255)
                                                                                                             ttAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATcc                                                  >  2:72712/1‑104 (MQ=255)
                                                                                                                aGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCtt                                                                                                         >  1:13312/1‑46 (MQ=255)
                                                                                                                aGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCtt                                                                                                         <  2:13312/46‑1 (MQ=255)
                                                                                                                aGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAc            <  1:1193053/139‑1 (MQ=255)
                                                                                                                aGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAc            <  2:1625843/139‑1 (MQ=255)
                                                                                                                  cATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcac                                                          <  2:159403/91‑1 (MQ=255)
                                                                                                                  cATGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcac                                                          >  1:159403/1‑91 (MQ=255)
                                                                                                                    tggccAAATCCTTATGGAAAAATGATTATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACt        <  1:1523839/135‑1 (MQ=255)
                                                                                                                    tGCGCAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACt        >  1:90122/1‑139 (MQ=255)
                                                                                                                     gcgcAAATCGTAATGCAAAAATGATAATAAATAcgc                                                                                                               <  2:302283/36‑1 (MQ=255)
                                                                                                                     gcgcAAATCGTAATGCAAAAATGATAATAAATAcgc                                                                                                               >  1:302283/1‑36 (MQ=38)
                                                                                                                     gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTg                                                                       >  1:1372727/1‑75 (MQ=255)
                                                                                                                     gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTg                                                                       <  2:1372727/75‑1 (MQ=255)
                                                                                                                     gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACg                                   >  2:425297/1‑111 (MQ=255)
                                                                                                                     gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACg                                   <  1:425297/111‑1 (MQ=255)
                                                                                                                     gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTa       >  2:1481652/1‑139 (MQ=255)
                                                                                                                     gcgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTa       >  1:683882/1‑139 (MQ=255)
                                                                                                                     gcgcAAATCGTAATGCAAAAATGATAATAATAAGCGGCCTTTAACCCGAAAGCCTGTCTTGGGGGTTTTTTTTTCCTTGTGGAATAAAAAAATTCCCCCTCCCCCTTCACCCTCTTATCAAGGGTATTTTGGACAAcat       >  2:1523839/1‑137 (MQ=255)
                                                                                                                      cgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt      <  1:1584019/139‑1 (MQ=255)
                                                                                                                      cgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt      <  2:683882/139‑1 (MQ=255)
                                                                                                                      cgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt      <  2:1152900/139‑1 (MQ=255)
                                                                                                                      cgcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGAGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt      <  1:1283165/139‑1 (MQ=255)
                                                                                                                       gcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGa                                                                                                >  2:289552/1‑48 (MQ=255)
                                                                                                                       gcAAATCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGa                                                                                                <  1:289552/48‑1 (MQ=255)
                                                                                                                          aaTCGTAATGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGc  <  2:1750458/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                    
CCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGC  >  NC_000913/2800346‑2800606

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