Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 854,765 G→A *373* (TAG→TAA)  ybiR → putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913854,7650GA98.0% 157.0 / ‑4.3 51*373* (TAG→TAA) ybiRputative transporter
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (28/22);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

GGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTG  >  NC_000913/854657‑854898
                                                                                                            |                                                                                                                                     
ggCGCACGATCGACGCATCTGGCGGCGTTTCCCTATCTATTCAATACCGATGCTGTTGTGGGCGGCGTCGGTGGGATATGTTTTGTTAGTTATCCCCCCGGCCAACTAAGGTATGGCAAATAAAGAAGGCGGATTATgt                                                                                                         >  1:117183/1‑138 (MQ=255)
     aCGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCa                                                                                                                              >  1:137926/1‑113 (MQ=255)
     aCGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCa                                                                                                                              <  2:137926/113‑1 (MQ=255)
     aCGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCg                                                                                                    >  1:926209/1‑139 (MQ=255)
       gATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCCATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGTCGGATTATGAGTCCGcc                                                                                                  >  2:990953/1‑139 (MQ=255)
                          gTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCa                                                                                                                              >  2:338667/1‑92 (MQ=255)
                          gTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCa                                                                                                                              <  1:338667/92‑1 (MQ=255)
                               cATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTg                                                                                 <  1:405557/132‑1 (MQ=255)
                               cATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTg                                                                                 >  2:405557/1‑132 (MQ=255)
                                               cGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCg                                                                                                                 <  2:1344659/84‑1 (MQ=255)
                                               cGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCg                                                                                                                 >  1:1344659/1‑84 (MQ=255)
                                               cGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGc                                                          >  1:1056103/1‑139 (MQ=255)
                                                gATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCa                                                         <  2:926209/139‑1 (MQ=255)
                                                 aTGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAgg                                                                                                                   <  1:952655/80‑1 (MQ=255)
                                                 aTGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAgg                                                                                                                   >  2:952655/1‑80 (MQ=255)
                                                 aTGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTg                                                                                 <  2:595929/114‑1 (MQ=255)
                                                 aTGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTg                                                                                 >  1:595929/1‑114 (MQ=255)
                                                 aTGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACg                                                                 <  2:182055/130‑1 (MQ=255)
                                                 aTGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACg                                                                 >  1:182055/1‑130 (MQ=255)
                                                 aTGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCAt                                                        >  2:1121343/1‑139 (MQ=255)
                                                 aTGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCAt                                                        >  1:549432/1‑139 (MQ=255)
                                                 aTGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCAt                                                        >  1:877082/1‑139 (MQ=255)
                                                        tgtgGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAg                                                 <  1:1264845/139‑1 (MQ=255)
                                                        tgtgGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAg                                                 <  2:1009751/139‑1 (MQ=255)
                                                        tgtgGGCGGCGCTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAg                                                 <  2:1056103/139‑1 (MQ=255)
                                                          tgGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAg                                                 >  1:819247/1‑137 (MQ=255)
                                                          tgGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAg                                                 <  2:819247/137‑1 (MQ=255)
                                                          tgGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCg                                               >  2:277264/1‑139 (MQ=255)
                                                            ggcggcGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAG‑TTTGAAGATTAATTCAGACGAACCGGCa                                                         <  1:417466/126‑1 (MQ=255)
                                                            ggcggcGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAG‑TTTGAAGATTAATTCAGACGAACCGGCa                                                         >  2:417466/1‑126 (MQ=255)
                                                             gcggcgTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGtt                                            >  2:1378578/1‑139 (MQ=255)
                                                                 cgTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAg                                                 <  1:1325789/130‑1 (MQ=255)
                                                                 cgTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAg                                                 >  2:1325789/1‑130 (MQ=255)
                                                                 cgTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTa                                        >  2:346248/1‑139 (MQ=255)
                                                                   ttGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAAt                                      <  1:1473542/139‑1 (MQ=255)
                                                                   ttGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAAt                                      <  1:346248/139‑1 (MQ=255)
                                                                                ttttGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCtt                                                                                                >  2:1467323/1‑68 (MQ=255)
                                                                                ttttGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCtt                                                                                                <  1:1467323/68‑1 (MQ=255)
                                                                                        gatataCTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTg                 >  2:1445928/3‑139 (MQ=255)
                                                                                                gtataagagacaGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACAt         >  1:794041/12‑139 (MQ=255)
                                                                                                cccGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACAt         >  1:199495/1‑139 (MQ=255)
                                                                                                  cGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCt                                  <  2:1129892/112‑1 (MQ=255)
                                                                                                  cGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCt                                  >  1:1129892/1‑112 (MQ=255)
                                                                                                   gggccaCTAAGCTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATCAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTg      >  1:936025/6‑139 (MQ=255)
                                                                                                   ggCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACg                                                                 >  2:1100348/1‑80 (MQ=255)
                                                                                                   ggCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACg                                                                 <  1:1100348/80‑1 (MQ=255)
                                                                                                   ggCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCt                                  >  1:1089636/1‑111 (MQ=255)
                                                                                                   ggCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAACTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCt                                  <  2:1089636/111‑1 (MQ=255)
                                                                                                      cAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTtcatc                             <  2:1572175/113‑1 (MQ=255)
                                                                                                      cAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTtcatc                             >  1:1572175/1‑113 (MQ=255)
                                                                                                       aaCTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTg  <  1:1121343/139‑1 (MQ=255)
                                                                                                            |                                                                                                                                     
GGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTG  >  NC_000913/854657‑854898

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: ‑