Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,478,324	0	T	G	42.3%	80.4 / 10.9	52	G33G (GGT→GGG)	ynbB	putative CDP‑diglyceride synthase
Reads supporting (aligned to +/- strand):  ref base T (9/21);  new base G (22/0);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 7.80e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACATTTT  >  NC_000913/1478175‑1478474                                                                                                                                                      |                                                                                                                                                        gACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTg                                                                                                                                                        <  2:833430/151‑1 (MQ=255)  aCTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGa                                                                                                                                                       <  2:820891/151‑1 (MQ=255)       gATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGGGAGAAAAAGGGGGGAgggt                                                                                                                                                   >  2:703721/1‑148 (MQ=255)           aaTCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACAtt                                                                                                                                              <  1:822752/151‑1 (MQ=255)           aaTCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGG‑TACAtt                                                                                                                                              <  1:744076/150‑1 (MQ=255)                    cGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCCGAGAGCAAAAGGGGGGAGGGTACATTAcaca                                                                                                                                         >  2:682220/1‑144 (MQ=255)                       aGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAGGGGGGGGGGTTACATTACGTAtt                                                                                                                                       >  1:395155/1‑146 (MQ=255)                            gCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGGGGAAAGGGGGGGAGGGTACATTACGTATTCGTACCtgg                                                                                                                              >  1:665094/1‑150 (MQ=255)                             cttcttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGAGCGTAAAGGGGGTGGGGGTACATTACGTATTCGTATCtggtg                                                                                                                            >  1:560614/1‑151 (MQ=255)                                 ttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATCGGGGAGAAAAGGGGGGGGGGGTATATTTCGTATTCGTACCTGGTgggg                                                                                                                         >  2:363666/1‑150 (MQ=255)                                  tATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTggg                                                                                                                          <  1:770832/148‑1 (MQ=255)                                   aTGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAAGGGGGGAGGGTACATTTCGTATTCGTACCTGGTgggg                                                                                                                         >  1:35639/1‑148 (MQ=255)                                     gAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACGCTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATc                                                                                                                    <  1:661701/150‑1 (MQ=255)                                        cgcTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGAGAGAAAAAGGGGGGGGGGTACATTACGTATTTGTACCTGGTGGGGG‑ATCGTc                                                                                                                 >  2:825171/1‑150 (MQ=255)                                         gcTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCAt                                                                                                               <  2:292791/151‑1 (MQ=255)                                            gAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTg                                                                                                            <  1:235900/151‑1 (MQ=255)                                             aaaTATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGt                                                                                                           <  1:507835/151‑1 (MQ=255)                                              aaTATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATCGGGGAGAAAAAGGGGTGTGGTTTCATTAAGGATTTTTATCTGGGGGGGG‑ATcattcttt                                                                                                             >  1:469132/1‑141 (MQ=255)                                               atatAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGGAAAAGGGGTGGGGGTTCATTGCGGATTTCTAACTCGGGGGGG‑ATCGGCCTTTGttt                                                                                                         >  1:568708/1‑151 (MQ=255)                                                   aaTGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTc                                                                                                      <  2:301602/150‑1 (MQ=255)                                                    aTGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTCaa                                                                                                    <  1:89061/151‑1 (MQ=255)                                                       cTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGGGAGGGTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTTAATgg                                                                                                 >  1:353320/1‑151 (MQ=255)                                                           aaaaaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGGGAGAAAAAGGGGGGGGGGTACATTTCGGATTTGTAACTGGGGGGGG‑ATCGGCATTTTTTTTTTTATgg                                                                                                 >  1:265097/1‑147 (MQ=255)                                                            aaaaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTACGGGGAGAAAAAGGGGGGGGGGGACATTTCGTATTTCTACCTGGGGGGGG‑GTCGGCCTTTTTTTTTTAATGGGGAtt                                                                                            >  2:786683/1‑151 (MQ=255)                                                             aaaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAGGGGGGGGGGGTACATTACGTATTCGTACCTGGTGGGGG‑GTCGTCCTTTGTTTTTTTATTGTGGttt                                                                                           >  2:808080/1‑151 (MQ=255)                                                               aTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTACGGGGAGAAAAGGGGGGGGGGGTTCATTATGTATTCGTACCTGGGGGGGG‑ATCGGCATTTGTTTTTTTATTGTGAttt                                                                                           >  2:117732/1‑149 (MQ=255)                                                                    tGGCAGCGCTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAGCGGTTTTCTGCTGTGGCGGTTGGCGGGGAGAAAAGGGGGTGGGGTTGCATTGCGTGTTCGTGCCTGGTGGGGG‑GTCGTGATTTGTTTTTGAATGGTGATGTcg                                                                                         <  2:250697/146‑2 (MQ=255)                                                                       cAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAGGGGGGGAGGGTACATTACCCATTTGTAGCTGGTGGGGG‑ATCGTCCTTTGTTTTTTAATGGGGATTTTAGGt                                                                                      >  2:716113/1‑146 (MQ=255)                                                                           acTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGGGGGAAAAAGGGGGGAGGGTACATTTCGTTTTTTGACCTTGGGGGGG‑ACCGTCATTTTTTTTTTAATGGTGGTTTTTGGTCCACGCTgg                                                                             >  2:538820/1‑151 (MQ=255)                                                                           acTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAGGGGGGGGGGGTACATTACGTATTCGTAACTGGTGGGGG‑ATCGTCATTTGTTTTTTTATGGGGGTTTCAGGTTCACGc                                                                                >  2:194055/1‑148 (MQ=255)                                                                             ttttCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGGAAAAGGGGTGGGGGTTCTATTCGGATTCGGACCTGGTGGGGC‑AGCGTCATTTTTTTTTCTATGGTGGTTTTAGGTCCACGCTggg                                                                            >  2:157359/1‑149 (MQ=255)                                                                                 cgcgTTGTTAATTTTAGCAACGCTGATATACCGTTTTCTGCTATGGCGGTTAACGGGGAGAAAAGGGGGGGGGGGTACATTACGTATTTGTTCCTGGGGGGGG‑ATCGTCATTTGTTTTTTTATGGTGATTTTAGGTCCACGCTGGGtgacg                                                                       >  2:282110/1‑151 (MQ=255)                                                                                  gcgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGGGAGGGTACAATACGTATTCGTACCTGGGGGGGG‑GTCGTCCTTTGTTTTTTAATCGGGATTTTAGGTCCACGCTGGGtgacgt                                                                      >  2:344098/1‑151 (MQ=255)                                                                                   cgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAGGGGGGGAGGGTACATTTCGTATTCGTACCCGGGGGGGG‑ATCGTCATTTTTTTTTTAATGGGGGTTTTAGGGCCACGCTGGGtgac                                                                        >  2:642886/1‑148 (MQ=255)                                                                                          aaTTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGttt                                                              <  1:208692/151‑1 (MQ=255)                                                                                                   aaCGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTg                                                     <  2:494551/151‑1 (MQ=255)                                                                                                         gATAAACCGTTTTCTGCTATGGCGGTTACCGGGGGGAAAAAGGGGGGGGGGTTCATTACGAAATCGTATCCGGTGGGGCGACGGTTTTTTGTTTTT‑AATGGGGGTTTTATGCTCCCGCCGGGTGAGGGTGGGGTTTTTTTCGCGGAtta                                                 >  2:761588/1‑147 (MQ=255)                                                                                                           tAAACCGTTTTCTGCTATGGCGGTTACCGGGGGGAAAAGGGGGGGGGGTTACATTACGTATTCGTACCTGGGGGGGC‑ATCGTCATTTTTTTTTTAATGGGGGTTTTAGGTCCACGCTGGATGACGGTGACGTTTTTTTCGCTGATAAGttt                                             >  1:579222/1‑151 (MQ=255)                                                                                                           tAAACCGTTTTCTGCTATGGCGGTTAACGGGGAGAAAAGGGGGGGAGGTTACATTACGTATTCTTACCTGGTGGGGG‑ATCGTCATTTGTTTTTTAATGGGGATTTTAGGGCCACGCTGGGTGAAGTTGACGTTTTTTGTGCTGATAAGttt                                             >  1:31634/1‑151 (MQ=255)                                                                                                            aaaCCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGttt                                             <  2:780601/150‑1 (MQ=255)                                                                                                            aaaCCGTTTTCTGCTATGGCGGTTAACCGGGAGAAAAGGGGGGGGGGGTTCATTACGTATTCGTACCTGGTGGGGG‑AACGTCCTTTTTTTTTTTATGGGGGTTTTAGGGCCCCGCCGGGTGACGTTTTCGGTTTTTTTGGTGCTAAAtttt                                            >  1:397508/1‑151 (MQ=255)                                                                                                               ccGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTg                                         <  1:529894/151‑1 (MQ=255)                                                                                                                cGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTgg                                        <  1:4980/151‑1 (MQ=255)                                                                                                                             ggCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAat                            <  2:151248/150‑1 (MQ=255)                                                                                                                                 gTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGt                       <  2:202633/151‑1 (MQ=255)                                                                                                                                      cGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCt                  <  2:175781/151‑1 (MQ=255)                                                                                                                                      cGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCt                  >  1:520820/1‑151 (MQ=255)                                                                                                                                       ggAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCtt                 <  2:52979/151‑1 (MQ=255)                                                                                                                                       ggAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCtt                 >  1:514830/1‑151 (MQ=255)                                                                                                                                             aaaaGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCt           >  2:746174/1‑151 (MQ=255)                                                                                                                                                   gggTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACAt     >  2:56119/1‑151 (MQ=255)                                                                                                                                                    ggTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACAtt    >  2:791580/1‑151 (MQ=255)                                                                                                                                                      tGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACAtttt  >  2:263304/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                        GACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGC‑ATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACATTTT  >  NC_000913/1478175‑1478474

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