Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	362,730	G→A	G150G (GGC→GGT)	lacY ←	lactose permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	362,730	0	G	A	97.1%	107.4 / NA	35	G150G (GGC→GGT)	lacY	lactose permease
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (17/17);  total (17/18)
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.73e-01

ACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTGAACAGTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGG  >  NC_000913/362485‑362955                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  aCGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTGAACAGTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCcaaccaac                                                                                                                                                                                                                              <  1:359919/251‑1 (MQ=255)                   gACAAAAACCACAGTGTTGGCTGTCTGAACAGGTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCcgcgc                                                                                                                                                                                                            <  1:35076/251‑1 (MQ=255)                         aaCCACAGTTTTGGCTGTCTGAACAGTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGAcc                                                                                                                                                                                                      <  2:156078/251‑1 (MQ=255)                              caGTTGTGGCTGTCTGAACAGTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAAtt                                                                                                                                                                                                 <  1:30020/251‑1 (MQ=255)                                        tGTCTGAACAGTTCCAGTGCCAGCTTAAGGCTACATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACt                                                                                                                                                                                       >  2:109790/1‑251 (MQ=255)                                         ggcTGAACAGTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATTGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACATCCAAACAGCCTCGCGCTACCAAATTCGAAATTACTg                                                                                                                                                                                      <  2:483470/249‑1 (MQ=255)                                         gTCTGAACAGTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACt                                                                                                                                                                                       <  1:320000/250‑1 (MQ=255)                                                cAGTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGc                                                                                                                                                                               >  2:4002/1‑251 (MQ=255)                                                  gTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTg                                                                                                                                                                             >  1:259299/1‑251 (MQ=255)                                                                        aaaTGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCctct                                                                                                                                                        <  1:580876/250‑1 (MQ=255)                                                                         aaTGCCGAATGGTTGGAACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCCACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCatctcc                                                                                                                                                      >  2:33510/1‑249 (MQ=255)                                                                                aaTGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACGGCGACGGCTGACTTGCGCAATAAATGCCTCGACTGCTGgc                                                                                                                                               <  2:197308/251‑1 (MQ=255)                                                                                 aTGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGgcg                                                                                                                                              <  1:184650/251‑1 (MQ=255)                                                                                  tGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGgcgc                                                                                                                                             >  1:249958/1‑251 (MQ=255)                                                                                   ggTTGGCACCTACCTCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCTACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGgcgc                                                                                                                                             <  1:504071/250‑1 (MQ=255)                                                                                         cACCTACCGCATTGGCAACCGTGGCAGAAGAAGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATTCCTCTACTGCTGGCGCACCGGCg                                                                                                                                      >  1:336060/1‑251 (MQ=255)                                                                                              aCCGCAGTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGATTAAAACGGCGATGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGAGGCCGACAATCGAGGCACACAGCGTCCAACCAACACAGCCAAACATCCGCGCGCTACCAAATTCGAAAGTGCTGCGTCGGCTGACATTCGCAATAAATGCCTCTACTGGTGGCGGACCCGCGTTaa                                                                                                                                  <  2:37783/250‑1 (MQ=255)                                                                                                  cATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACa                                                                                                                              >  1:85327/1‑250 (MQ=255)                                                                                                   aTTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACaaa                                                                                                                            <  1:22275/251‑1 (MQ=255)                                                                                                               ggCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGCGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGCTTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCCACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATaa                                                                                                                 >  1:92795/1‑250 (MQ=255)                                                                                                                             cATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATc                                                                                                  <  1:333858/251‑1 (MQ=255)                                                                                                                                ccGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGAt                                                                                               <  1:378269/251‑1 (MQ=255)                                                                                                                                                   tAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGttgtatt                                                                            >  2:470043/1‑251 (MQ=255)                                                                                                                                                     aaaCGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGttgtattgt                                                                          <  1:517151/251‑1 (MQ=255)                                                                                                                                                           cGAGGATGAGTGCACAGCCAGAGCCAAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCACATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATACATACCACCAACAATCGATCCTCCTAAAATGTTGTATTGTAACAGt                                                                    >  1:330692/1‑251 (MQ=255)                                                                                                                                                             aGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGCCTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTgg                                                                  >  1:105792/1‑251 (MQ=255)                                                                                                                                                                          aGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATGCCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATaaa                                                      >  1:343610/1‑250 (MQ=255)                                                                                                                                                                                              acaaacTGATTATAGATGGTGAACATGATGCCGACACTCGAGGAACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACAt                                 >  2:153520/1‑251 (MQ=255)                                                                                                                                                                                                        tattGATGGTGAACATGAGGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAGACATCACTAACAt                        <  2:92795/250‑1 (MQ=255)                                                                                                                                                                                                               ggTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGc                                                                                                                                                            >  1:91301/1‑111 (MQ=255)                                                                                                                                                                                                               ggTGAACAGGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCTACCAAAGGCGAAATGACTGCGACGGCTGACTTTCGCAATAAATGc                                                                                                                                                            <  2:91301/111‑1 (MQ=255)                                                                                                                                                                                                                gTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGtaat               >  2:122719/1‑251 (MQ=255)                                                                                                                                                                                                                      atgatgCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGAGAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCca         <  1:33510/251‑1 (MQ=255)                                                                                                                                                                                                                       tgatgCCGACAATCGAGGCACACAGCGCCCCACCAACACAGACAACCATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCcac        >  1:456746/1‑251 (MQ=255)                                                                                                                                                                                                                        gatgCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCcac        >  1:5859/1‑250 (MQ=255)                                                                                                                                                                                                                             cGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGAACAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAgg  <  2:423411/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  ACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTGAACAGTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGG  >  NC_000913/362485‑362955

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