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

GGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACATTT  >  NC_000913/1478174‑1478473                                                                                                                                                        |                                                                                                                                                      ggACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGt                                                                                                                                                       <  2:919338/151‑1 (MQ=255)           tAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGA‑GGGAAAAGGGGGGGAgggttaat                                                                                                                                             >  2:479558/1‑145 (MQ=255)            aaTCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATCGGG‑GAGAAAAGGGGGGGAGGGTACAtt                                                                                                                                            >  2:277850/1‑151 (MQ=255)             aTCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGG‑GCGAAAAGGGGGGGAGGTTACa                                                                                                                                              >  1:314464/1‑148 (MQ=255)               cGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTAc                                                                                                                                          <  2:945231/150‑1 (MQ=255)                  tGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGA‑GAGGAAAGGGGGGGAGGGTTCATTACGTa                                                                                                                                       >  1:684126/1‑150 (MQ=255)                    tCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCTGAGGGGAAAAGGGGGGGGGGTTTAATT‑CGTAtt                                                                                                                                     >  2:258754/1‑150 (MQ=255)                     cGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGG‑GGGGAAAGGGGGGGGGGGTACATTACGTATTc                                                                                                                                    >  1:117723/1‑150 (MQ=255)                     cGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTc                                                                                                                                    <  1:932666/150‑1 (MQ=255)                        aGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTAc                                                                                                                                <  1:488518/151‑1 (MQ=255)                             gCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGA‑GAGAAAAGGGGGGGGGGGTACATTACGTATTCGTACCtgg                                                                                                                            >  1:192072/1‑150 (MQ=255)                             gCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCtgg                                                                                                                            <  1:177429/150‑1 (MQ=255)                               ttcttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATCGGA‑GGGAAAAGGGGGGGAGGGTACATTACGTATTCGTAACttgtgg                                                                                                                         >  2:551055/1‑151 (MQ=255)                                tcttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCtggtgg                                                                                                                         <  2:423585/150‑1 (MQ=255)                                 cttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGG‑GAGGAAAGGGGGGGAGGGTACATTACGTCTTCGTTCCCTGTgggg                                                                                                                       >  2:736446/1‑151 (MQ=255)                                  ttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCCGG‑GAGAAAAAGGGGGGAGGTTACATTACGGATTTGTACCTGGTgggg                                                                                                                       >  1:251457/1‑150 (MQ=255)                                     tGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGG‑GAGAAAAAGGGGGGGGGGTACATTACGTTTTTGTACCTGGGGGGGGAt                                                                                                                    >  1:373803/1‑150 (MQ=255)                                        gcgcTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGG‑GAGGAAAAGGGGGGAGGGTACATTACGTATTCGTAACTGGGGGGGGGTCGt                                                                                                                 >  2:373584/1‑150 (MQ=255)                                         cgcTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGG‑GAGAAAAGGGGGTGAGGGTACATTACGTATTCGTACCTGGTGGGGGATCGTc                                                                                                                >  2:484377/1‑150 (MQ=255)                                             gAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGG‑GAGAAAAGGGGGGGGGGGTACATTTCGTATTTTTACCTGGTGGGGGATCGTCAttt                                                                                                            >  1:547177/1‑150 (MQ=255)                                              aaaTATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGt                                                                                                          <  2:444991/151‑1 (MQ=255)                                                  ataATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGttttt                                                                                                      <  1:384125/151‑1 (MQ=255)                                                      tGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCCGG‑GAGAAAAGGGGGGGAGGGTTAATTACGTATTTGTAACTGGTGGGGGGTCGTCATTTGttttttta                                                                                                   >  2:789055/1‑147 (MQ=255)                                                        cTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGG‑GAGGAAAAGGGGGGGGGGTAAATTACGTATTTGTACCTGGGGGGGGATCGTCATTTTTTTTTTAATgg                                                                                                >  2:532720/1‑151 (MQ=255)                                                        cTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGG‑GAGAAAAAGGGGGGGGGGTTCATTTCGGATTCGTAACCGGGGGGGCCTCGTCCTTTTTTTTTTAATgg                                                                                                >  1:518670/1‑151 (MQ=255)                                                             aaaaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGC‑GGGAAAAGGGGGGGGGGGTACATTAAGTATTCGTACCTGGTGGGGGATGGGCATTTGTTTTTTAAt                                                                                                  >  1:830563/1‑144 (MQ=255)                                                               aaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATCGGG‑GAGGAAAAGGGGGGGGGGTAAATTATGGATTCGTACCCGGTGGGGGGTCGTCATTTTTTTTTTTATGGGGAtttt                                                                                         >  2:70402/1‑150 (MQ=255)                                                               aaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTc                                                                                         <  1:421930/151‑1 (MQ=255)                                                               aaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTc                                                                                         <  2:256653/151‑1 (MQ=255)                                                                             cTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGA‑GGGAAAAGGGGGGGAGGGTACATTACGTATTCGTACCTGGTGGGGGATCGTCTTTTGTTTTTTTATGGGGGttttcggtacccg                                                                                >  2:93351/1‑136 (MQ=255)                                                                              ttttCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATCGGC‑AAGAAAAGGGGGGGATGGTACATTTTGTATTTGTACCCGGTGGGGGGTCGTAATTTGTTTTTTTATTGGGAGTTCGGGTCCACGCTGGAt                                                                          >  1:97325/1‑151 (MQ=255)                                                                               tttCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCCGA‑GAGAAAAAGGGGGGGGGGTACATTACGTATTCGTACCTGGGGGGGCGACGTCATTTGTTTTTTCATGGGGATTTTTGGTCCACGc                                                                               >  2:443489/1‑145 (MQ=255)                                                                                    cgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacgtt                                                                    <  1:840469/151‑1 (MQ=255)                                                                                              tttAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGG‑GAGAAAAGGGGGGGAGGGTACATTTCGGATTTGTACCTGGGGGGGGATCGTCATTTGTTTTTTTATGGGGGTTTTAGGGCCACCCCGGGTGGAGTTGGCGgttttt                                                          >  1:332550/1‑151 (MQ=255)                                                                                               ttAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTACGGA‑GGGAAAAAGGGGGGAGGGTACATTAAGTATTTTTACCTGGGGGGGGATTGGCATTTTTTTTTTAATGGTGATTTTAGGGCCCCGCTGGATGGCGTTGGTGtttttt                                                          >  1:786704/1‑150 (MQ=255)                                                                                               ttAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGG‑GAGAAAAGGGGGGGAGGGTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTTAATGGTGGTTTCAGGGCCCCGCTGGATGACGTTGACGtttttt                                                          >  1:115258/1‑150 (MQ=255)                                                                                                       gCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATa                                                 <  1:361274/151‑1 (MQ=255)                                                                                                       gCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATa                                                 <  1:445074/151‑1 (MQ=255)                                                                                                          gATAAACCGTTTTCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATaa                                                <  1:467848/149‑1 (MQ=255)                                                                                                          gATAAACCGTTTTCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAg                                               <  1:759306/150‑1 (MQ=255)                                                                                                              aaCCGTTTTCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTc                                          <  2:764051/151‑1 (MQ=255)                                                                                                              aaCCGTTTTCTGCTATGGCGGCTACCGGA‑GAGAAAAAGGGGGGAGGGTACATTACGTATTCGTATCTGGGGGGGGATCGGCATTTGTTTTTCAATGGGGATTTTAGGTCCACGCCGGGTGACGCTGAGGTTTTTTTCGCTGATAAGtttt                                           >  2:50018/1‑150 (MQ=255)                                                                                                                      tCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCAtt                                   <  1:277854/150‑1 (MQ=255)                                                                                                                              ggCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAat                           <  2:185643/150‑1 (MQ=255)                                                                                                                                 ggTTACCTTA‑GAGAAAAGGGTTTAAGGTTACATTACGTATTCTTACCTGGTGTGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACTTTGACTTTTTTTGCGCTTATCATTTTTCTTGCATTgcttgcatatt                        <  2:701508/150‑11 (MQ=255)                                                                                                                                 ggTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTg                       <  2:788012/151‑1 (MQ=255)                                                                                                                                     aCCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACg                   >  2:466655/1‑151 (MQ=255)                                                                                                                                       cGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCt                 >  2:107805/1‑151 (MQ=255)                                                                                                                                          a‑gagaAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTtat              <  1:789072/151‑1 (MQ=255)                                                                                                                                                aaaGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTg         <  2:251453/151‑1 (MQ=255)                                                                                                                                                   gggggTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTAc      <  1:660815/151‑1 (MQ=255)                                                                                                                                                   gggggTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTAc      >  1:477269/1‑151 (MQ=255)                                                                                                                                                       gTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACAttt  <  2:103301/151‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                      GGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGA‑GAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACATTT  >  NC_000913/1478174‑1478473

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