Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,194,572	(T)6→5	intergenic (+47/+149)	yqiK → / ← sibD	protein YqiK/small RNA SibD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,194,567	0	T	.	100.0%	121.4 / NA	30	intergenic (+42/+154)	yqiK/sibD	protein YqiK/small RNA SibD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (14/16);  total (14/16)

TAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAG  >  NC_000913/3194488‑3194642                                                                                |                                                                            tAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtt                                                                            >  2:593453/1‑81 (MQ=255)    ccTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAttttt                                                                         <  2:593411/81‑1 (MQ=255)    ccTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAttttt                                                                         <  1:593410/81‑1 (MQ=255)          caacTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCt                                                                  <  1:593412/81‑1 (MQ=255)            acTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGc                                                                >  1:593461/1‑81 (MQ=255)                ccGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGt                                                              <  1:593413/79‑1 (MQ=255)                    cGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGc                                                         >  2:593450/1‑80 (MQ=255)                     gaagaaAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCt                                                        >  2:593458/1‑80 (MQ=255)                     gaagaaAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTa                                                       <  1:593414/81‑1 (MQ=255)                      aagaaAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTaa                                                      <  1:593415/81‑1 (MQ=255)                      aagaaAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTa                                                       >  1:593445/1‑80 (MQ=255)                         aaaaaGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCTGGGGA‑TTTTTATATCTGCGTTCCGCTaaa                                                     >  2:593447/1‑79 (MQ=25)                            aaGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTg                                                >  1:593460/1‑81 (MQ=255)                             aGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGc                                               <  2:593416/81‑1 (MQ=255)                              gCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGc                                               >  2:593459/1‑80 (MQ=255)                                agagTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGCaaa                                            <  2:593417/81‑1 (MQ=255)                                 gagTAAACCCTCTTTCAGGAATAATCGCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGCAAAt                                           <  1:593418/81‑1 (MQ=18)                                 gagTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGCAAAt                                           <  2:593419/81‑1 (MQ=255)                                  agTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATg                                          <  2:593420/81‑1 (MQ=255)                                       aCCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAg                                     <  1:593422/81‑1 (MQ=255)                                       aCCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAg                                     <  2:593421/81‑1 (MQ=255)                                        cccTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAgg                                    >  1:593452/1‑81 (MQ=255)                                               ccAGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgttgt                                <  2:593423/78‑1 (MQ=255)                                                 aGGAATAATCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTgcag                           >  1:593473/1‑81 (MQ=255)                                                    aataatCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCAt                        <  1:593424/81‑1 (MQ=255)                                                    aataatCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCAt                        >  2:593469/1‑81 (MQ=255)                                                        atCCCTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCAtttt                     <  1:593425/80‑1 (MQ=255)                                                          cccTGGCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCg                  >  2:593468/1‑81 (MQ=255)                                                               gCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGa               <  2:593427/79‑1 (MQ=255)                                                               gCCCAAAAAGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACg             <  2:593426/81‑1 (MQ=255)                                                                    aaaaGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCg        >  1:593463/1‑81 (MQ=255)                                                                     aaaGCCGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCgg       >  1:593471/1‑81 (MQ=255)                                                                          cGGGGA‑TTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAg  >  1:593470/1‑81 (MQ=255)                                                                                |                                                                            TAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAG  >  NC_000913/3194488‑3194642

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