Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	986,184	C→A	29.7%	intergenic (+50/‑200)	U712_04810 → / → U712_04815	putative oxidoreductase yhdF/putative amino acid permease YhdG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	986,184	0	C	A	29.7%	89.9 / 10.4	56	intergenic (+50/‑200)	U712_04810/U712_04815	putative oxidoreductase yhdF/putative amino acid permease YhdG
Reads supporting (aligned to +/- strand):  ref base C (7/31);  new base A (10/5);  total (20/36)
Fisher's exact test for biased strand distribution p-value = 2.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCTATGTTCTGCTGGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGACCATTG  >  CP006881/986047‑986318                                                                                                                                             |                                                                                                                                       agcCTATGTTCTGCTGTCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCa                                                                                                                                       <  5:200754/138‑1 (MQ=255)              ctgGCGTCTGACGAATCTTCCTATATTAAAGGGCACACCAATCAAGTGAATTGCGGCCGAGTTATTTAAACCTAAACAAACAAAAAAAACCGAAACAGCCTGTTTTTTTTG‑TGTTGA‑‑AAAATTTTAAAAAAAAAAAtaa                                                                                                                          >  6:200754/1‑137 (MQ=255)              ctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTTAAAGTAATAAAACAAAAAACACCGCAAAAGCTGGGTTTTTTTT‑TTGTGT‑‑AAAATTTTAAAATAAAAAAatt                                                                                                                          >  2:233895/1‑139 (MQ=255)              ctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTAAACCTAATAAACCAAAAAACACCACAAGAGCTGGGTTTTTGTG‑TGGTGA‑‑AAAAATTAAAAATAAAAAAttt                                                                                                                          >  2:84631/1‑139 (MQ=255)              ctgGCGTCTGACGAATATTCCTATATGACAGGGCAGACCAATCATGTGAATGGCGGCCGGGTTTTTTTAAACTATAAAAACAAAAAAAAACACAACAGCGGGGTTTTTGTG‑TGGGGAATAAAATATAAAAAAAAAAAAt                                                                                                                            >  6:370408/1‑139 (MQ=255)                   gTCTGACGAATCTTCCTATATTACAGGGCAGCCCATTCATGTGAAAGGCGGCCGGTTTTTTTGAAACGAATAAAAGAAAAAACACCAAAAGAACGGGTTTTTTTTGTTGTTAA‑‑AAAAATTAAAAAAAAAAAAttttaat                                                                                                                      >  8:8356/1‑135 (MQ=255)                         cGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTc                                                                                                               <  2:27980/139‑1 (MQ=255)                                          cAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTAAAACTAATAAAACAAAAAACAACACAAGAGCTGGGTTTTTTTG‑TTGTGT‑‑AAAAATTAAAAAAAAAAAAATTTTATATTTTTTTAAAAAATtttttt                                                                                                 >  4:39647/1‑136 (MQ=255)                                             ggCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTAAACGTAATAAACCAAAAACCAACGCAATAGGTGGGTTTTTTTG‑TGGTTA‑‑AAAATTTAAAAAAAAAAAAATTTTGGAGTTTTATTAAAAATTT‑TTTTTAtttt                                                                                          >  8:91291/1‑139 (MQ=255)                                             ggCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAAAGTTATAAACGAAAAAACATCTCCAGATCTGGTGTTTTGTG‑TGGTTT‑‑AAAATTTCAAAATAAAAAATTTTTCGATTTTTCGAACAATTTTTTTTTGAttt                                                                                           >  8:345264/1‑139 (MQ=255)                                             ggCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTAAACCTAATAAAAGAAAAACCAGCAAAAGAGCGGGTTTTTTTTG‑TGGTGA‑‑AAAATTTAAAAATAAAAAATTTTTGGATTTTTATAACAAAttttttttgtttt                                                                                           >  5:380543/1‑134 (MQ=255)                                              gCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTATACCTAATAAAAGAAAAAAAAACACAAGAGCGGGGTTTTTTTG‑TGGTGA‑‑AAAAATTAAAAAAAAAAAATTTTTAGATTTTTTTTAAAATTTTTTTTTTAt                                                                                             >  4:394265/1‑136 (MQ=255)                                                  aCCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCa                                                                                      <  1:420021/139‑1 (MQ=255)                                                  aCCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCa                                                                                      <  5:182805/139‑1 (MQ=255)                                                  aCCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCa                                                                                      <  2:192052/139‑1 (MQ=255)                                                  aCCATTCATGTGAATGGAGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGc                                                                                       >  5:236751/1‑138 (MQ=255)                                                  aCCATTCATGAGAATGGCGGCCGTTTTATTTAAACAGAATGAACGAAAAAACACCAAAACAACGGGTTTTTTTTT‑TTGTGT‑‑AAAAATAAAAAATTAAAAAATTTTATATTTTTTAAAAAATTTTTTTTTTAatttttc                                                                                       >  6:325305/1‑136 (MQ=255)                                                     aTTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtt                                                                                   <  7:157526/139‑1 (MQ=255)                                                     aTTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtt                                                                                   <  7:147132/139‑1 (MQ=255)                                                       tCATGTGAATGGCGGCCGTTTTTTTTAAAACTAATAAAAGAAAAACAAGCTGAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑TAAATTTTAAAATAAAAAAATTGTCGATTTTTATAACAATTTTGTTTTTATTTTTGGAAtttt                                                                                 >  7:136511/1‑139 (MQ=255)                                                       tCATGTGAATGGCGGCCGTTTTATTTAAAAGTAATCAACCAAAAACCAACAAAAGCGCTGGGGTTTTGTT‑GGGGGA‑‑AAAATTTCAAAATAAAAAATTTGTCGTTTTTTGTAAAAATTTTTTTTTTAATTTTGTAAttat                                                                                 >  8:406126/1‑137 (MQ=255)                                                        cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                 <  1:84631/138‑1 (MQ=255)                                                        cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                 <  7:345264/138‑1 (MQ=255)                                                        cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                 <  5:325305/138‑1 (MQ=255)                                                        cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAt                                                                                    <  4:260219/135‑1 (MQ=255)                                                         aTGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtat                                                                               <  7:91291/139‑1 (MQ=255)                                                              aaTGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtata                                                                              <  4:352653/135‑1 (MQ=255)                                                              aaTGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCGGGGTTTTTTTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTATAACAATTTTGTTTTCATTTTTTCAATTTtatatat                                                                           >  7:127451/1‑138 (MQ=255)                                                                 ggcggcCGTGTTTTTTAAACGTAATCAACGAAAAAACAGCGCAAGCGCGGGGTTTTTTTG‑TGGTGT‑‑AAAATTTAAAAATAAAAAATTTTTAGATTTTTGTAACAATTTTTTTTTTATTTTTTGAATTTtttatat                                                                           >  3:352653/1‑135 (MQ=255)                                                                  gcgggCGTTTTATATCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatatata                                                                          <  3:79287/135‑1 (MQ=255)                                                                  gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                       <  7:247366/138‑1 (MQ=255)                                                                  gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                       <  6:380543/138‑1 (MQ=255)                                                                      ccGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATAtt                                                                     <  4:83616/136‑1 (MQ=255)                                                                      ccGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg                                                                  <  7:193591/139‑1 (MQ=255)                                                                        gTTTTATTTAAAAGTAATAAAAAAAAAAACACCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑TAAAATTAAAAATAAAAAATTTTTCGATTTTTTTAAAAAATTTTTTTTAATTTTTTGAATTTTATATATAATAttta                                                                   >  4:378085/1‑134 (MQ=255)                                                                         ttttATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAg                                                               <  6:236751/139‑1 (MQ=255)                                                                                      tAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCAttt                                                                                           <  4:418483/98‑1 (MQ=255)                                                                                      tAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCAttt                                                                                           >  3:418483/1‑98 (MQ=255)                                                                                              gAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                            <  8:355587/105‑1 (MQ=255)                                                                                              gAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                            >  7:355587/1‑105 (MQ=255)                                                                                               aaaaaCCTGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTcgcg                                            <  3:378085/136‑1 (MQ=255)                                                                                               aaaaaCCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTTA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGAAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTcgcg                                            <  3:39647/136‑1 (MQ=255)                                                                                               aaaaaCCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTcgcg                                            <  3:394265/136‑1 (MQ=255)                                                                                               aaaaaCCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                         <  5:113887/139‑1 (MQ=255)                                                                                               aaaaaCCAGCTAAATTGCTGGTTTTTTGTG‑TGGTTA‑‑AAAATTTAAAAATAAAAAATTTGTCGATTTTTAAAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTTAGTATGAGAAAGGAATGATTCGCGAcc                                         <  7:8356/139‑1 (MQ=255)                                                                                                aaaaCCAGCTCATGAGCTGGTTTTTTGTG‑TGGTTA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGATAAAGCAATGATTcgcg                                            <  4:185831/135‑1 (MQ=255)                                                                                                aaaaCCAGCTAAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                         <  8:197957/138‑1 (MQ=255)                                                                                                 aaaCCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGCAGAAAGGAATGATTCGCGACCAg                                       <  1:233895/139‑1 (MQ=255)                                                                                                 aaaCCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTAAAAATAAAAATTTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTTGGAGAACGGAATGATTCGCGACCAg                                       <  5:370408/139‑1 (MQ=255)                                                                                                  aaCCAGGTCAAGTCGTGGTTTTTTTTG‑TGGTTA‑‑AAAATTTAAAAATAAAAAATTTTTCGTTTTTTCTAAAAATTTTGTTTTCATTTTTGCAATTTTATATATAATTTTTATTAGGATAAAGGAATGATTCGCGACCAGt                                      <  2:18060/139‑1 (MQ=255)                                                                                                    ccAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGAAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTAc                                    <  5:348873/139‑1 (MQ=255)                                                                                                                  ggTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATAAAGTAGaa                       >  2:381604/1‑138 (MQ=255)                                                                                                                   gTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAAAAAAAAATTTGTCGATTTTTCAAACAATTTTGTTTTAATTTTTGAAATTTAATATATAAAATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTTCAAATATAGTAGaaaa                     <  2:5696/139‑1 (MQ=255)                                                                                                                                   a‑‑aaaaTTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATTATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGAc       <  8:38895/138‑1 (MQ=255)                                                                                                                                   a‑‑aaaaTTTAAAAATAAAAATTTTGTGGTTTTTTCTAACAATTTTGTTTTCATTTTTTCATTTTTAATTATATAATTGTTTTGGCGAAAGCATTGTTTCGCGACCAGTACAAATATAGTAGAAAAATAAGAATAATGAc       <  6:261365/138‑1 (MQ=255)                                                                                                                                         aTTTAAAAATAAAAAATTTGTCTATTTTTCTAAAAATTTTGTTTTCATTTTTGAAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGACCATTg  <  8:127451/139‑1 (MQ=255)                                                                                                                                             |                                                                                                                                       CGCCTATGTTCTGCTGGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTG‑TGGTGA‑‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGACCATTG  >  CP006881/986047‑986318

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