Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,926,096 G→A *219* (TAG→TAA yobB → C‑N hydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,926,0960GA98.0% 152.9 / NA 49*219* (TAG→TAAyobBC‑N hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (28/20);  minor base C (1/0);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTATGGGATGAATATGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGT  >  NC_000913/1925984‑1926229
                                                                                                                |                                                                                                                                     
cACTATGGGATGAATATGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATg                                                                                                             <  1:127898/139‑1 (MQ=255)
   tATGGGATGAATATGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTg                                                                                                          >  1:589151/1‑139 (MQ=255)
       ggATGAATATGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCa                                                                                                      <  1:345322/139‑1 (MQ=255)
          tGAATATGGTCGTCTCATAGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTa                                                                                                   >  1:507578/1‑139 (MQ=255)
           gAATATGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTAt                                                                                                  <  2:423847/139‑1 (MQ=255)
             atatGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCg                                                                                                <  2:626411/139‑1 (MQ=255)
             atatGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCg                                                                                                >  2:126266/1‑139 (MQ=255)
                 ggTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATg                                                                                                                  >  2:395973/1‑117 (MQ=255)
                 ggTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATg                                                                                                                  <  1:395973/117‑1 (MQ=255)
                 ggTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATAc                                                                                            >  1:15926/1‑139 (MQ=255)
                 ggTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATAc                                                                                            >  1:441929/1‑139 (MQ=255)
                      tctcATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGATCATGCCATGTTGCGCATTATCGATACAGaaa                                                                                       >  2:11825/1‑139 (MQ=255)
                      tctcATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGt                                                                                                            >  1:487375/1‑118 (MQ=255)
                      tctcATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGt                                                                                                            <  2:487375/118‑1 (MQ=255)
                      tctcATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGaaa                                                                                       >  2:150489/1‑139 (MQ=255)
                           tCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGc                                                                                  <  2:199667/139‑1 (MQ=255)
                               tCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGtt                                                                              <  2:648624/139‑1 (MQ=255)
                               tCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGtt                                                                              >  1:648624/1‑139 (MQ=255)
                                cgcgCCGAACGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTACGCTTTTTCAGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGttt                                                                             >  2:377201/1‑139 (MQ=255)
                                 gcgcCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTg                                                                            <  2:589151/139‑1 (MQ=255)
                                               tCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGt                                                              >  2:556142/1‑139 (MQ=255)
                                                 attattGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGcc                                                                                                                           >  1:629545/1‑76 (MQ=255)
                                                 attattGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGcc                                                                                                                           <  2:629545/76‑1 (MQ=255)
                                                 attattGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTg                                                            >  1:433114/1‑139 (MQ=255)
                                                   tattGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCttttt                                                                                                                                >  1:672233/1‑69 (MQ=255)
                                                   tattGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCttttt                                                                                                                                <  2:672233/69‑1 (MQ=255)
                                                       gTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGtt                                                                                                           >  1:289900/1‑86 (MQ=255)
                                                       gTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGtt                                                                                                           <  2:289900/86‑1 (MQ=255)
                                                       gTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTAt                                                                                                  >  1:189126/1‑95 (MQ=255)
                                                       gTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTAt                                                                                                  <  2:189126/95‑1 (MQ=255)
                                                        ttAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGc                                                     >  2:332914/1‑139 (MQ=255)
                                                          aGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCct                                                   >  1:307654/1‑139 (MQ=255)
                                                           gTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCctc                                                  <  1:556142/139‑1 (MQ=255)
                                                               ttCAGGGTCCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCATCTTCGGTCTGGAGCATGCCAAGTTGCCCATTATGGGTACAGCAACCTGAGGTTTGCAGCGAGGGATCGTCATGGTTGCATCTGtt                                              >  2:499325/2‑139 (MQ=255)
                                                               gTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGtt                                              >  1:712237/1‑139 (MQ=255)
                                                                           cacaGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGTAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGAc                                  <  1:315702/139‑1 (MQ=255)
                                                                               gggTTGGGAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTga                                                                                                                        <  1:473638/49‑2 (MQ=255)
                                                                               gggTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTgg                                                                                                                        >  2:473638/1‑49 (MQ=255)
                                                                                     gCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCg                        >  2:163088/1‑139 (MQ=255)
                                                                                        aGGCGATATCATTCCATTACGCTCAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGTAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCa                     <  2:507578/139‑1 (MQ=255)
                                                                                               atCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCAt              >  1:642772/1‑139 (MQ=255)
                                                                                                 cATTCCATTACGCTAAGCTTTTTCGGCCTGGATCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGa            <  1:11825/139‑1 (MQ=255)
                                                                                                    tCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGcc         <  1:150489/139‑1 (MQ=255)
                                                                                                     ccATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCa        >  2:239609/1‑139 (MQ=255)
                                                                                                        ttACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCt     >  2:735582/1‑139 (MQ=255)
                                                                                                           cGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGt  >  2:473243/1‑139 (MQ=255)
                                                                                                           cGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATAAGCCACCTGGt  >  1:732234/1‑139 (MQ=255)
                                                                                                             cTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGc                                                     >  2:339923/1‑86 (MQ=255)
                                                                                                             cTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGc                                                     <  1:339923/86‑1 (MQ=255)
                                                                                                                |                                                                                                                                     
CACTATGGGATGAATATGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGT  >  NC_000913/1925984‑1926229

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