Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	548,374	G→A	A6A (GCG→GCA)	ylbE →	NO‑induced DUF1116 protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	548,374	0	G	A	94.5%	224.7 / 3.4	73	A6A (GCG→GCA)	ylbE	NO‑induced DUF1116 protein
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base A (38/31);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

GCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAG  >  NC_000913/548278‑548470                                                                                                 |                                                                                                 gCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACaa                                                                                               >  2:1533375/1‑100 (MQ=255) gCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACaa                                                                                               >  1:497393/1‑100 (MQ=255)   aTTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGc                                                                                             <  1:1185014/100‑1 (MQ=255)   aTTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGCTTGCAATAAGGGGTACCCCTGTCTACATCCGTCGCACAAGc                                                                                             >  2:1489805/1‑100 (MQ=255)     tACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCa                                                                                           >  2:1309563/1‑100 (MQ=255)      aacAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCaa                                                                                          >  1:1065913/3‑100 (MQ=255)       ccAATGGTCGCCAGTCGCCGGTGGCAATAAAAACCTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCGAAt                                                                                         >  2:1265121/1‑100 (MQ=255)         aaTGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAAtgc                                                                                       >  1:102201/1‑100 (MQ=255)           tGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAAtgctg                                                                                     >  2:1381275/1‑100 (MQ=255)             gTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCg                                                                                   >  1:654757/1‑100 (MQ=255)               cGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGt                                                                                 <  2:568813/100‑1 (MQ=255)               cGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGt                                                                                 >  2:440300/1‑100 (MQ=255)                   aGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATc                                                                             >  2:1111086/1‑100 (MQ=255)                      caccGGTGGCAATAAAAAACTGTCTCTTTTATTAGAACGTTTGCAATAGGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCg                                                                            <  2:1043088/96‑1 (MQ=255)                          ggTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACaaa                                                                      >  2:1361637/1‑100 (MQ=255)                            tGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAAtt                                                                    <  1:1348580/100‑1 (MQ=255)                              gCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATtcg                                                                  <  1:1369912/100‑1 (MQ=255)                              gCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGAACAAGCCAATGCTGCGGTTATCGAACAAATtcg                                                                  <  2:1023669/100‑1 (MQ=255)                               cAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATtcg                                                                  >  2:1407508/1‑99 (MQ=255)                                naTAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATtcgtc                                                                >  1:1197225/2‑100 (MQ=255)                                aaTAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATtcgtc                                                                >  2:433411/1‑100 (MQ=255)                                 aTAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATtcgtcg                                                               <  1:862355/100‑1 (MQ=255)                                  tAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTcgc                                                              <  2:1396299/100‑1 (MQ=255)                                  tAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGGTATCGAACAAATTCGTcgc                                                              >  2:1416338/1‑100 (MQ=255)                                     aaaaCTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCt                                                           <  1:173946/100‑1 (MQ=255)                                       aaCTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCg                                                         <  1:282512/100‑1 (MQ=255)                                       aaCTGGCTCATTTAGTAGAACGTTCGCAATAATAGATTCCCATTCTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCCCg                                                         <  2:1470994/100‑1 (MQ=255)                                        aCTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGt                                                        <  1:664460/100‑1 (MQ=255)                                          tGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTcc                                                      >  2:798693/1‑100 (MQ=255)                                             cTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCcaca                                                   >  2:1025660/1‑100 (MQ=255)                                             cTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCcaca                                                   >  1:453604/1‑100 (MQ=255)                                             cTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCcaca                                                   >  2:1326665/1‑100 (MQ=255)                                                 tttATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACActgg                                               >  1:21872/1‑100 (MQ=255)                                                  ttattaGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACActggc                                              <  1:1519373/100‑1 (MQ=255)                                                  ttattaGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACActggc                                              <  1:442668/100‑1 (MQ=255)                                                   tattaGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACActggct                                             <  1:945253/100‑1 (MQ=255)                                                    attaGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACActggctn                                            <  2:887753/100‑2 (MQ=255)                                                    attaGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACActggctg                                            >  2:877417/1‑100 (MQ=255)                                                     ttaGAACGTTTGCAATAAGGCGTTCGCATGTTGACAAGAGTGGCACAAGCCAATGCTGCGGTCATAGAACAAATTCGGCGCGCTCGTCCACActggctg                                            <  2:1405649/99‑1 (MQ=255)                                                      taGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGa                                          <  2:8333/100‑1 (MQ=255)                                                       aGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGTACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGAt                                         >  1:711276/1‑100 (MQ=255)                                                       aGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGAt                                         >  1:1124830/1‑100 (MQ=255)                                                       aGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGAt                                         >  2:1267425/1‑100 (MQ=255)                                                       aGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGAt                                         >  2:1209070/1‑100 (MQ=255)                                                          aCGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGAtgtg                                      <  2:654754/100‑1 (MQ=255)                                                             tttGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCaa                                   >  1:528619/1‑100 (MQ=255)                                                             nttGCAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCaa                                   >  1:1354681/2‑100 (MQ=255)                                                                 cAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCgg                               >  1:624187/1‑100 (MQ=255)                                                                 cAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCgg                               <  1:618882/100‑1 (MQ=255)                                                                 cAATAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCgg                               >  1:940389/1‑100 (MQ=255)                                                                  aaTAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGc                              >  2:80024/1‑100 (MQ=255)                                                                   aTAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGc                              >  2:1007262/1‑99 (MQ=255)                                                                    tAAGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGctt                            >  2:1233764/1‑100 (MQ=255)                                                                     aaGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGcttc                           >  2:860947/1‑100 (MQ=255)                                                                     aaGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGcttc                           <  2:102201/100‑1 (MQ=255)                                                                     aaGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGcttc                           <  2:215679/100‑1 (MQ=255)                                                                     aaGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGcttc                           <  2:333925/100‑1 (MQ=255)                                                                      aGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTTGCGCTCGTCCACACTGGCTGGATGTGCAACCGGcttct                          <  1:925945/100‑1 (MQ=255)                                                                      aGGGGTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGcttct                          <  2:1376376/100‑1 (MQ=255)                                                                        gggTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCttcttc                        <  1:1309572/100‑1 (MQ=255)                                                                        gggTTCCCATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCttcttc                        <  2:526711/100‑1 (MQ=255)                                                                         ggTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCa                       >  2:537518/1‑100 (MQ=255)                                                                             cccATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTa                   <  2:1065906/100‑1 (MQ=255)                                                                               naTGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATc                 >  2:1430870/2‑100 (MQ=255)                                                                               cATGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATc                 >  1:262815/1‑100 (MQ=255)                                                                                 tGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAg               <  1:433413/100‑1 (MQ=255)                                                                                 tGTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAg               <  1:1460369/100‑1 (MQ=255)                                                                                  gTTTACATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGc              >  1:1190843/1‑100 (MQ=255)                                                                                      aCATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAAc          <  2:452891/100‑1 (MQ=255)                                                                                       cATCAGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACt         <  2:124206/100‑1 (MQ=255)                                                                                           aGTGGCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAAc     >  1:1173526/1‑100 (MQ=255)                                                                                              ggCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAg  >  1:70329/1‑100 (MQ=255)                                                                                              ggCACAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAg  >  2:821925/1‑100 (MQ=255)                                                                                                 |                                                                                                 GCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAATGCTGCGGTTATCGAACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAG  >  NC_000913/548278‑548470

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑