Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,381,915	G→T	8.0%	intergenic (+13/‑32)	ycjV → / → ompG	pseudogene; putative ATP‑binding component of a transport system/outer membrane porin G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,381,915	0	G	T	8.0%	166.3 / 5.4	50	intergenic (+13/‑32)	ycjV/ompG	pseudogene; putative ATP‑binding component of a transport system/outer membrane porin G
Reads supporting (aligned to +/- strand):  ref base G (23/23);  new base T (2/2);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-01

CCACGGTTGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAA  >  NC_000913/1381767‑1382062                                                                                                                                                     |                                                                                                                                                    ccacGGTTGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAg                                                                                                                                                     >  4:345586/1‑149 (MQ=255)  cacGGTTGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAgg                                                                                                                                                    >  4:394986/1‑149 (MQ=255)  cacGGTTGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAgg                                                                                                                                                    >  2:50453/1‑149 (MQ=255)        tGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTcc                                                                                                                                              <  2:429865/149‑1 (MQ=255)         gggggCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTccc                                                                                                                                             <  2:177762/149‑1 (MQ=255)          ggggCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGATAATATCACTATTCATTTTGATATGACGAAATGACATTTCTTTGATTCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAATGCCTTcccc                                                                                                                                            <  4:430026/149‑1 (MQ=255)           gggCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTccccc                                                                                                                                           <  3:549234/149‑1 (MQ=255)           gggCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTccccc                                                                                                                                           >  3:288900/1‑149 (MQ=255)            ggCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAATGCCTTcccc                                                                                                                                            <  2:557328/147‑1 (MQ=255)            ggCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAATGCCTTcccc                                                                                                                                            >  1:557328/1‑147 (MQ=255)                 gAGTTAGTGGTCCGTGCTGGTGCGTTACATGATTATAATGCAGGAGAAAATATCAATATTCATGTTGATATGACGAAATGTCATGTCTTTGATGCAGAAACTGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGata                                                                                                                                     >  2:159559/1‑149 (MQ=255)                    ttAGTTGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGataata                                                                                                                                  <  1:50453/149‑1 (MQ=255)                                  tGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACaataat                                                                                                                    <  2:273843/149‑1 (MQ=255)                                      gCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGaaa                                                                                                                >  4:133695/1‑149 (MQ=255)                                        gTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATAGAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGaaaaa                                                                                                              <  1:385550/149‑1 (MQ=255)                                        gTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACaa                                                                                                                               <  1:113538/132‑1 (MQ=255)                                        gTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACaa                                                                                                                               >  2:113538/1‑132 (MQ=255)                                        gTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGaaaaa                                                                                                              <  2:272492/149‑1 (MQ=255)                                          tAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGt                                                                                                            <  3:582492/149‑1 (MQ=255)                                             aTGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGttat                                                                                                         >  4:333428/1‑149 (MQ=255)                                             aTGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGttat                                                                                                         >  2:471856/1‑149 (MQ=255)                                             aTGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGttat                                                                                                         >  1:5079/1‑149 (MQ=255)                                                aTTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTAc                                                                                                      <  4:215798/149‑1 (MQ=255)                                                  tATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTAccc                                                                                                    >  1:357562/1‑149 (MQ=255)                                                   aTCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCt                                                                                                   <  2:357562/149‑1 (MQ=255)                                                    tCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTg                                                                                                  >  4:481922/1‑149 (MQ=255)                                                      aTGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTa                                                                                                <  3:394986/149‑1 (MQ=255)                                                                aaaTATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGt                                                                                      <  1:112158/149‑1 (MQ=255)                                                                  atatCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGa                                                                                    <  3:198057/149‑1 (MQ=255)                                                                   tatCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAgg                                                                                                                             <  2:498900/107‑1 (MQ=255)                                                                   tatCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAgg                                                                                                                             >  1:498900/1‑107 (MQ=255)                                                                         tATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGc                                                                             <  1:552532/149‑1 (MQ=255)                                                                         tATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGc                                                                             >  3:331262/1‑149 (MQ=255)                                                                           ttCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCgg                                                                           <  1:393791/149‑1 (MQ=255)                                                                              aTTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGaa                                                                        >  2:233367/1‑149 (MQ=255)                                                                               ttttGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTgcgcaggcca                                                                       >  2:40432/1‑143 (MQ=255)                                                                               ttttGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAAt                                                                       >  4:332068/1‑149 (MQ=255)                                                                                                gTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGCGAATGGCCTGCGCACAGGCCg                                                      >  2:256737/1‑149 (MQ=255)                                                                                                 tCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGa                                                     >  3:81850/1‑149 (MQ=255)                                                                                                    tttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGagga                                                  >  4:215667/1‑149 (MQ=255)                                                                                                           gATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAAc                                           >  4:328373/1‑149 (MQ=255)                                                                                                           gATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAAc                                           >  1:461681/1‑149 (MQ=255)                                                                                                            aTGCAGAAACGGAAATAGCACTTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTAACGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACg                                          >  1:393297/1‑149 (MQ=255)                                                                                                              gCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGAc                                        <  3:133695/149‑1 (MQ=255)                                                                                                                 gAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTgg                                     <  2:5079/149‑1 (MQ=255)                                                                                                                           tAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTaa                             <  2:198529/147‑1 (MQ=255)                                                                                                                           tAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTaa                             >  1:198529/1‑147 (MQ=255)                                                                                                                             gCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTTGCACTTTAATATc                         <  3:247146/149‑1 (MQ=255)                                                                                                                                      tAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATg                >  2:557268/1‑149 (MQ=255)                                                                                                                                             aGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGaaa         <  4:295352/149‑1 (MQ=255)                                                                                                                                              gggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAAt        >  1:492483/1‑149 (MQ=255)                                                                                                                                                  gAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGaa    >  2:10233/1‑149 (MQ=255)                                                                                                                                                    aGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTATGAAATAGaaaa  <  3:333428/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CCACGGTTGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAA  >  NC_000913/1381767‑1382062

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