Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,193	(C)6→7	coding (16/1353 nt)	creD →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,187	1	.	C	93.8%	207.9 / 8.4	64	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base C (32/28);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638093‑4638281                                                                                                |                                                                                                   gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑c‑c‑cc                                                                                            <  2:1114157/100‑1 (MQ=255)   ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑c‑c‑cccc                                                                                          >  2:1202711/1‑99 (MQ=255)   ttCGACTTCACCATCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑c‑c‑cccc                                                                                          >  1:1024991/1‑99 (MQ=255)    tagACGTCACCGTCACTTCACATAGCTCCAAATTCGGCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑c‑c‑cc                                                                                            <  1:560383/95‑1 (MQ=255)     cGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGATGAAAT‑c‑c‑c‑cccctg                                                                                        >  2:1795310/1‑97 (MQ=255)        cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTc                                                                                     >  2:668806/1‑100 (MQ=255)         ttCACCGTCACTTCACATAGCTTCAAATTCATCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑CGC‑CCC‑TGTTCt                                                                                    >  1:1193882/1‑100 (MQ=255)          tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTg                                                                                   >  1:878715/1‑100 (MQ=255)           cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTgg                                                                                  >  1:309482/1‑100 (MQ=255)           cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTg                                                                                   >  1:358020/1‑99 (MQ=255)            aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGa                                                                                 <  2:646721/100‑1 (MQ=255)            aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGa                                                                                 <  1:335604/100‑1 (MQ=255)              tgTCACTTCACATAGCTTCAAATTCTGCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGaaa                                                                               <  2:1184636/99‑1 (MQ=255)               gTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGaa                                                                                <  1:1575071/98‑1 (MQ=255)                tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAAt                                                                             <  1:97382/100‑1 (MQ=255)                tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAAt                                                                             >  1:580643/1‑100 (MQ=255)                     tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTa                                                                        >  2:1530031/1‑100 (MQ=255)                      cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAg                                                                       >  2:1728421/1‑100 (MQ=255)                        caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGcc                                                                     >  1:1020767/1‑100 (MQ=255)                          tAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTg                                                                   >  1:32008/1‑100 (MQ=255)                           aGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGt                                                                  >  1:1387605/1‑100 (MQ=255)                            gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTg                                                                   <  1:1100179/98‑1 (MQ=255)                            gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGtt                                                                 >  2:99743/1‑100 (MQ=255)                             cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGttt                                                                <  1:1425903/100‑1 (MQ=255)                              ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  2:1312223/1‑100 (MQ=255)                               tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTgg                                                              >  2:1549278/1‑100 (MQ=255)                               tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTgg                                                              <  1:1450598/100‑1 (MQ=255)                               tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTgg                                                              >  1:1492811/1‑100 (MQ=255)                                cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  1:674018/1‑100 (MQ=255)                                cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  1:1724532/1‑100 (MQ=255)                                   aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCa                                                          <  1:1443811/100‑1 (MQ=255)                                      cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGt                                                        >  1:355399/1‑99 (MQ=255)                                      cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                       >  2:775080/1‑100 (MQ=255)                                      cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAAAAACTAGACTGATTGGTGCAGTa                                                       >  2:732645/1‑100 (MQ=255)                                        tcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAtt                                                     <  1:234685/100‑1 (MQ=255)                                           cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  <  2:990407/100‑1 (MQ=255)                                            acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 <  1:1017555/100‑1 (MQ=255)                                            acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 <  2:106786/100‑1 (MQ=255)                                            acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 >  1:1768558/1‑100 (MQ=255)                                                 ttcttctTTTCCTTCTGCCATTGTAAATGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  2:784668/99‑1 (MQ=255)                                                 gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  2:1025001/100‑1 (MQ=255)                                                 gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  1:1463670/100‑1 (MQ=255)                                                  tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  2:1785730/100‑1 (MQ=255)                                                    ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         <  1:893119/100‑1 (MQ=255)                                                     tCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAtt                                        >  2:1572295/1‑100 (MQ=255)                                                       gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  1:41537/100‑1 (MQ=255)                                                        tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     <  2:1318324/100‑1 (MQ=255)                                                        tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     >  2:1742779/1‑100 (MQ=255)                                                        tATCCTGCTGACATTGCAAAGTAGAAGACTATGTTGACAT‑C‑C‑CAAAA‑AAAAAAAAAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  1:1196558/99‑1 (MQ=255)                                                          tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  2:1652133/100‑1 (MQ=255)                                                             tgctgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                                >  1:1702741/1‑100 (MQ=255)                                                               ctgcCATTGCTAAGGAGAAGATTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                              <  2:1013184/100‑1 (MQ=255)                                                               ctgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                              <  2:358022/100‑1 (MQ=255)                                                               ctgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                              <  1:1763877/100‑1 (MQ=255)                                                                 gcCATTGCAAAGGAGAAGACTATGTTGAAAG‑CGC‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  1:842058/99‑1 (MQ=255)                                                                    aTTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TTTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                         <  2:1057294/100‑1 (MQ=255)                                                                     ttGCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                        >  1:657621/1‑100 (MQ=255)                                                                       gCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTCTTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  1:1783730/1‑100 (MQ=255)                                                                       gCAAAGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  1:1026096/1‑100 (MQ=255)                                                                          aaGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                   >  1:1716951/1‑100 (MQ=255)                                                                           aGGAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAgg                  <  1:762437/100‑1 (MQ=255)                                                                            ggAGAAGACTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAgtt                 <  2:1233485/100‑3 (MQ=255)                                                                               gaagaCTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              >  2:843502/1‑100 (MQ=255)                                                                                 agaCTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  1:394331/1‑100 (MQ=255)                                                                                    cTATGTTGAAATCC‑C‑C‑CCC‑TGTTCTTGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         <  2:32008/100‑1 (MQ=255)                                                                                       tGTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  2:1197432/1‑100 (MQ=255)                                                                                        gTTGAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:1119972/1‑100 (MQ=255)                                                                                         ttGAAAT‑C‑C‑C‑CCCGTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGCTGTTGATTCCGATAATGCTGATTCGACAGGTGATTGTCGaa       >  2:1220746/1‑97 (MQ=255)                                                                                           gAAATCC‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  >  1:595407/1‑100 (MQ=255)                                                                                                |                                                                                                   GCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑C‑C‑CCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638093‑4638281

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