Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,212,703 C→T *135* (TAG→TAA ycgX ← DUF1398 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,212,7030CT96.5% 171.7 / ‑0.6 57coding (405/405 nt)ycgXDUF1398 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/27);  minor base . (1/1);  total (29/28)
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.22e-01

TAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGCTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGCCAGTGTTGAACAGTAT  >  NC_000913/1212569‑1212838
                                                                                                                                      |                                                                                                                                       
tAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGa                                                                                                                                     >  2:149193/1‑139 (MQ=255)
   tCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTa                                                                                                                                  <  1:149193/139‑1 (MQ=255)
   tCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTa                                                                                                                                  <  2:374018/139‑1 (MQ=255)
       gTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATgg                                                                                                                              >  1:951512/1‑139 (MQ=255)
            tatttatGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTttat                                                                                                                         <  2:827314/139‑1 (MQ=255)
            tatttatGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCAATTATGCTAATGTCTTACTGAATTATGTAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTttat                                                                                                                         >  2:743218/1‑139 (MQ=255)
                 atGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTAttt                                                                                                                    >  1:253189/1‑139 (MQ=255)
                    ggTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTTCACTTTTATTTCCATTATGCTAATTCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTaa                                                                                                                                 <  2:619679/123‑1 (MQ=255)
                    ggTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACGGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTaa                                                                                                                                 >  1:619679/1‑123 (MQ=255)
                     gTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTa                                                                                                                <  2:951512/139‑1 (MQ=255)
                        aCAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACAt                                                                                                             >  2:853729/1‑139 (MQ=255)
                           aTAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCtt                                                                                                          <  1:853729/139‑1 (MQ=255)
                             aaaaaGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCa                                                                                                        <  2:253189/139‑1 (MQ=255)
                                      gTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATaaa                                                                                               <  2:954217/139‑1 (MQ=255)
                                       tATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATaaaa                                                                                              >  1:818793/1‑139 (MQ=255)
                                                 ccATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTAc                                                                                                               <  1:630416/112‑1 (MQ=255)
                                                 ccATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATc                                                                                    >  2:502591/1‑139 (MQ=255)
                                                 ccATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCGTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTAc                                                                                                               >  2:630416/1‑112 (MQ=255)
                                                   aTAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCaa                                                                                  >  2:630137/1‑139 (MQ=255)
                                                   aTAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCaa                                                                                  >  2:101683/1‑139 (MQ=255)
                                                      gagaTAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATgg                                                                               <  1:101683/139‑1 (MQ=255)
                                                         aTAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTaa                                                                                                                                 <  2:444678/86‑1 (MQ=255)
                                                         aTAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTaa                                                                                                                                 >  1:444678/1‑86 (MQ=255)
                                                           aGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCgtag                                                                          >  1:661733/1‑139 (MQ=255)
                                                              ctctGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTttatta                                                                                                                       <  1:746519/91‑1 (MQ=255)
                                                              ctctGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTttatta                                                                                                                       >  2:746519/1‑91 (MQ=255)
                                                               tctGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACAt                                                                                                             <  1:772314/100‑1 (MQ=255)
                                                               tctGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACAt                                                                                                             >  2:772314/1‑100 (MQ=255)
                                                               tctGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCtt                                                                                                          >  2:323321/1‑103 (MQ=255)
                                                               tctGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCtt                                                                                                          <  1:323321/103‑1 (MQ=255)
                                                                        ttATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGt                                                             >  2:484083/1‑139 (MQ=255)
                                                                          aTTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTa                                                           <  1:502591/139‑1 (MQ=255)
                                                                                aTTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTa                                                     <  1:484083/139‑1 (MQ=255)
                                                                                         aTGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGc                                                                                            <  1:901409/91‑1 (MQ=255)
                                                                                         aTGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGc                                                                                            >  2:901409/1‑91 (MQ=255)
                                                                                         aTGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTa                                            >  2:520298/1‑139 (MQ=255)
                                                                                               tACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCa                                                                                   <  2:505143/94‑1 (MQ=255)
                                                                                               tACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCa                                                                                   >  1:505143/1‑94 (MQ=255)
                                                                                               tACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCAc                                      >  2:969995/1‑139 (MQ=255)
                                                                                               tACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAAc                                                   <  1:415693/125‑1 (MQ=255)
                                                                                               tACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAAc                                                   >  2:415693/1‑125 (MQ=255)
                                                                                                aCTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAAc                                                   >  2:339048/1‑125 (MQ=255)
                                                                                                aCTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAAc                                                   <  1:339048/125‑1 (MQ=255)
                                                                                                      ttATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCg                                                                             <  2:981968/93‑1 (MQ=255)
                                                                                                      ttATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCg                                                                             >  1:981968/1‑93 (MQ=255)
                                                                                                           aaGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGccc                          >  1:776905/1‑139 (MQ=255)
                                                                                                                   cTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGc                  <  2:661733/139‑1 (MQ=255)
                                                                                                                   cTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGc                  <  2:791944/139‑1 (MQ=255)
                                                                                                                   cTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGc                  <  2:812554/139‑1 (MQ=255)
                                                                                                                   cTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTATACATCTTCAATATACAAAAGCCTATTATCAATGCCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGc                  <  1:743218/139‑1 (MQ=255)
                                                                                                                        gTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGCCAgtg             <  1:896077/139‑1 (MQ=255)
                                                                                                                           tCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGCCAGTGTTg          >  1:477572/1‑139 (MQ=255)
                                                                                                                            ccAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGCCAGTGTTGa         >  2:256640/1‑139 (MQ=255)
                                                                                                                            ccAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGCCAGTGTTGa         <  1:520298/139‑1 (MQ=255)
                                                                                                                             cAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGCCAGTGTTGaa        <  1:273308/139‑1 (MQ=255)
                                                                                                                                   tAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCg                                                                  >  2:358243/1‑75 (MQ=255)
                                                                                                                                   tAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGCCAGTGTTGAACAGTAt  >  1:586979/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                       
TAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGCTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGCCAGTGTTGAACAGTAT  >  NC_000913/1212569‑1212838

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