Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,112,578	0	T	G	44.0%	132.2 / 17.5	76	G88G (GGT→GGG)	yiiR	DUF805 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (3/39);  new base G (33/0);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 3.51e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTCGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAG  >  NC_000913/4112436‑4112725                                                                                                                                               |                                                                                                                                                         tACTCGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACCgggg                                                                                                                                                  >  1:515828/1‑150 (MQ=255)  aCTCGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTgg                                                                                                                                                    <  2:1178265/148‑1 (MQ=255)  aCTCGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGc                                                                                                                                                  <  2:658721/150‑1 (MQ=255)  aCTCGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCg                                                                                                                                                 <  2:1295787/151‑1 (MQ=255)  aCTCGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGTACTGGGCg                                                                                                                                                 >  2:726231/1‑151 (MQ=255)    tCGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACCggggg                                                                                                                                                 >  1:1090950/1‑147 (MQ=255)    tCGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAATggg                                                                                                                                                   >  2:1047033/1‑147 (MQ=255)     cGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑tt                                                                                                                                            <  2:199544/151‑1 (MQ=255)      gATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGCAACTGGGGG‑A‑ttt                                                                                                                                           >  1:1154067/1‑151 (MQ=255)       atatTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGGTGCTGCTGGCGGGGAAATGGGGG‑G‑tttt                                                                                                                                          >  1:109718/1‑151 (MQ=255)       atatTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGGG‑G‑tttt                                                                                                                                          >  1:412809/1‑151 (MQ=255)         atTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTa                                                                                                                                         <  1:45882/150‑1 (MQ=255)          ttCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGGG‑G‑TTTTTcc                                                                                                                                       >  2:371510/1‑151 (MQ=255)           tCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGGG‑G‑ttttttc                                                                                                                                       >  2:890663/1‑147 (MQ=255)            cAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGGTGCTGCTGGCGGGGAACTggggcgttttttc                                                                                                                                         >  1:68582/1‑139 (MQ=255)            cAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGCG‑C‑tttttacgg                                                                                                                                     >  1:250219/1‑146 (MQ=255)             aGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAACGGGCgggttttactgg                                                                                                                                      >  1:897545/1‑140 (MQ=255)                   cggcgTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGGG‑G‑TTTTACC‑GGgggt                                                                                                                                >  2:173277/1‑148 (MQ=255)                         tttGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAAt                                                                                                                      <  2:386443/150‑1 (MQ=255)                                 gtgtgCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTg                                                                                                              <  1:721307/150‑1 (MQ=255)                                 gtgtgCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTg                                                                                                              <  2:845129/150‑1 (MQ=255)                                   gtgCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGGGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGGCGGGGGGGAAGTGGGGG‑ACTTTTACC‑GGGGGTG‑GGC‑GAATgggggggggg                                                                                                            >  1:902772/1‑151 (MQ=255)                                   gtgCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTggg                                                                                                            <  2:440695/150‑1 (MQ=255)                                     gcttgctCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGGG‑G‑TTTTTCC‑CGGGGTT‑TGG‑GAATGgggggggggggg                                                                                                         >  1:1230070/1‑149 (MQ=255)                                     gcttgctCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAATTGGGG‑G‑TTTTTCC‑GGGGGTC‑TGG‑GAATGgggggggggggg                                                                                                         >  2:930440/1‑149 (MQ=255)                                        tgctCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAATTGGCG‑G‑TTTTTTC‑GGGGGGC‑TGG‑CAAAGgggggggggggggtt                                                                                                      >  2:155275/1‑146 (MQ=255)                                           tctGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGGTGCTGCTGGCGGGGAAATGTGGG‑C‑TTTTTCC‑GGGGGTT‑TGG‑CGAAGGGTGGGGGGGGGgttttt                                                                                                   >  2:1092554/1‑149 (MQ=255)                                            ctGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTg                                                                                                    <  1:922656/149‑1 (MQ=255)                                              ggCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTc                                                                                                  <  1:540442/149‑1 (MQ=255)                                              ggCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTc                                                                                                  <  2:1019818/149‑1 (MQ=255)                                                    cAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGtt                                                                                          <  1:1070936/151‑1 (MQ=255)                                                     aGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGGG‑A‑TTTTTTC‑CGGGGTCTTGG‑CAATGGCGGGGGGGGGGtttttccc                                                                                                 >  1:1114810/1‑140 (MQ=255)                                                       cGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAATGGGCG‑G‑TTTTTCC‑GGGGGGT‑TGG‑GAATGGGGGGGGGGGGGGTTTGTTCCGGCggttgt                                                                                       >  2:760209/1‑146 (MQ=255)                                                                  gTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAATTGGCG‑A‑TTTTTCC‑GGGGGTC‑TGG‑GAATGGGCGGGGGGGGGGTTTGTCCCCACGTTGGTTTTgg                                                                                  >  2:279115/1‑145 (MQ=255)                                                                    tAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGCG‑G‑TTTTATC‑GGGGGGC‑TGG‑CAATGGGGGGGGGGGGGGTTTGTTCCCGCGTTTGTTTTGGGggtgat                                                                           >  1:167098/1‑150 (MQ=255)                                                                      aGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACAGGGGG‑G‑TTTTTAC‑GGGGGTT‑TGG‑GCATGGGGGGGGGGGGGGTTTGTCCCGACGTTGATTCTtgtg                                                                                >  1:283947/1‑143 (MQ=255)                                                                       gcgcTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAATGGGGG‑G‑TTTTTCC‑GGGGGTC‑TGG‑CAATGGGCGGGGGGGGGGTTTTTCCCGACGGTGATTTTGGGGGTGATGCtt                                                                       >  2:1125984/1‑151 (MQ=255)                                                                        cgcTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTa                                                                      <  1:100805/151‑1 (MQ=255)                                                                          cTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAATGGGCG‑G‑TTTTTTC‑CGGTGTT‑TGG‑GAATGGGGGGTGGGGGGGTTTTTCCCGACGTTTATTTTGGGGATGATGCtt                                                                       >  2:741633/1‑148 (MQ=255)                                                                              cATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATc                                                                <  2:520841/151‑1 (MQ=255)                                                                                tGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGCGCT‑TTTTACC‑GGGGTGC‑TGG‑CATAGGGC‑GTGGGGCGTTTTTTCCCCACCTTTGTTTTGGTGGTTATTCTTTTCggtttt                                                              >  2:1152293/1‑146 (MQ=255)                                                                                    cgcgGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACATGGCG‑G‑TTTTTCC‑GGGGGGC‑TTG‑CGATGGGGGGGGGGGGGGTTTTTCCCCACGTTGATTCTGGGGGGGATGCTTTTTGATCTTGGt                                                           >  1:149271/1‑150 (MQ=255)                                                                                        ggCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGtt                                                      <  1:1228692/151‑1 (MQ=255)                                                                                        ggCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGtt                                                      <  1:558217/151‑1 (MQ=255)                                                                                            ttCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTAGGAG‑C‑TTTTACC‑CGGTGGA‑TTG‑CAAAAGGCGGTGGGGGGTTTTGGTCCGAAGCTGGATCTGGGGAAGGAGCTTATCTATCCTGGTGCGCttttt                                                  >  2:514132/1‑151 (MQ=255)                                                                                             tCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGCG‑G‑TTTTTCC‑GGGGGTC‑TGGCCAATGGGGGGGGGGGGGTTTT‑TTCCGACGGTGATTTTTGGGATGATGCTTATCGGTCTTGGTGCGTttgtt                                                  >  2:209781/1‑150 (MQ=255)                                                                                             tCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGCG‑G‑TTTTTCC‑CGGTGTT‑TGG‑GCAAGGGCGGGGGGGCGGTTTTTCCCGCCGTTGATTTTGGGGGTGGTGCTTATTGATTTTgggggggtt                                                     >  1:753427/1‑140 (MQ=255)                                                                                               cGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAACGGGGG‑G‑TTTTTTC‑GGGGGTC‑TGG‑CGATGGGCGGGGGGGCGGGTTTTCCCGCCGTTTATTTTGGGGGGGATGTTTTTCGGTCTTGGTGCGTttgtt                                                  >  1:154755/1‑148 (MQ=255)                                                                                                   gcATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACCGGGCG‑G‑TTTTTCCGGGGGGTT‑GGG‑CAAGGGGGGGGGGGGGGTTTT‑TTCCGACGGTGATTTTTGGGGTGATGGTTTTTGATCTTTGTGGGGTTGTTTGgaca                                            >  2:217347/1‑150 (MQ=255)                                                                                                     aTGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGcacaca                                          <  2:869919/150‑1 (MQ=255)                                                                                                       gggCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAgg                                       <  1:1087182/151‑1 (MQ=255)                                                                                                       gggCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAgg                                       <  2:897586/151‑1 (MQ=255)                                                                                                       gggCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAg                                        <  1:190029/150‑1 (MQ=255)                                                                                                        ggCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGc                                      <  1:472746/151‑1 (MQ=255)                                                                                                            tttCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGCG‑G‑TTTTTCC‑GGGGGTC‑TGG‑CAATGGGCGGGGGGGGGTTTTTTCCCCACGTTGATTTTGGGGGTGATGCTTTTCGATCTTTGTGCGGTTTGTGGGACACAAgggg                                     >  2:661102/1‑146 (MQ=255)                                                                                                             ttCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGaaa                                  <  1:964390/150‑1 (MQ=255)                                                                                                                 gatgatTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAaaatg                               <  2:1185995/149‑2 (MQ=255)                                                                                                                 gatgatTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAaaataaa                             <  1:401399/151‑1 (MQ=255)                                                                                                                     atTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAatat                          <  1:1131591/150‑1 (MQ=255)                                                                                                                     atTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATg                         <  1:1043313/151‑1 (MQ=255)                                                                                                                     atTGTGGCGTGGATGCTGCTGGCGGGGTACCGGGGG‑G‑TTTTTCC‑CGGTGTT‑TGG‑GCATGGGGGGGGGGGGGGTTTGTTCCGACGGTGGTTTTTGGGATGATTGTTTTTGATCTTGGTGGGTTTGTTTGGACAACAGGCGGaaataaa                             >  1:902113/1‑147 (MQ=255)                                                                                                                     atTGTGCCGTGGATGGTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGTTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCAAACAAGGCTAAaataaata                           <  1:761255/149‑1 (MQ=255)                                                                                                                       tgtgGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATgg                        <  1:350715/150‑1 (MQ=255)                                                                                                                          ggCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTaa                     <  1:719410/150‑1 (MQ=255)                                                                                                                             gTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGAt                 <  1:371484/151‑1 (MQ=255)                                                                                                                               ggATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATAc               <  2:222062/151‑1 (MQ=255)                                                                                                                                 aagctgctGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGGGGGGCGTTTTGTCCCGACGTTGGTTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACt              >  1:1042375/3‑150 (MQ=255)                                                                                                                                 aTGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGGGCGTTTGTTGGCACACAAGGCGAACATAAATATGGTAAAGAt                 <  1:550337/147‑1 (MQ=255)                                                                                                                                      gctgGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGa         <  1:813710/150‑1 (MQ=255)                                                                                                                                      gctgGCGGGGAACTGGGCG‑A‑TTTTACC‑CGGTGTC‑TGG‑CAATGGGGGGTGGGGGGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGGGTTTGTTGGCACACAAGGGGAAAATAAATATGGTAAAGATACTCAgg          >  1:372945/1‑149 (MQ=255)                                                                                                                                        tgGCGGGTAACTGGGCC‑C‑TGTTACC‑CGGCGGC‑CGG‑GAATGGGCGGTGGGGCGTTTTGTCGCGAGGTTCATTCTGGTGGAGCAGCTTATCGATCTTGATGAGTTTGTTGGCACACAACGCGAACATAAGTATGGTAAAGATACTCAGGAtg       >  1:666283/1‑150 (MQ=255)                                                                                                                                         ggCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGAtgt      <  1:706884/150‑1 (MQ=255)                                                                                                                                         ggCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGAtgt      <  2:131233/150‑1 (MQ=255)                                                                                                                                         ggCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGAtgt      <  2:258478/150‑1 (MQ=255)                                                                                                                                           cGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGaa   <  2:1280860/151‑1 (MQ=255)                                                                                                                                            gggTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAg  <  1:281121/151‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                         TACTCGATATTCAGACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCG‑A‑TTTTACC‑GGGTGTC‑TGG‑CAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAG  >  NC_000913/4112436‑4112725

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