Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,842,914	0	A	G	38.2%	81.5 / 12.2	57	M182V (ATG→GTG)	gdhA	glutamate dehydrogenase, NADP‑specific
Reads supporting (aligned to +/- strand):  ref base A (2/32);  new base G (20/1);  total (24/33)
Fisher's exact test for biased strand distribution p-value = 9.08e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGTTTTGA  >  NC_000913/1842765‑1843062                                                                                                                                                      |                                                                                                                                                           aTTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGgat                                                                                                                                                           <  2:299533/151‑1 (MQ=255) aTTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGgat                                                                                                                                                           <  2:806430/151‑1 (MQ=255)    tCGATCCGAAAGGAAAAATCGAAGGTGAAGTGTTGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGATGTTCCGGCAGGTGATATCGGGGTTGGTAGTCGTGAAGTCGGCTTTATGGGGGGgatg                                                                                                                                                          <  1:1066424/149‑1 (MQ=255)    tCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGgatga                                                                                                                                                         <  2:11436/150‑1 (MQ=255)    tCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGgatga                                                                                                                                                         <  1:341169/150‑1 (MQ=255)      gctCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGAAAGAACTGTATCGCCAACTGGGCGCGGATCCCGGCGTTCCGGCAGGTGAGATCGGGGTTGGTGGTCTTTACGTCGGCTTTATGGAGGGgatg                                                                                                                                                          <  1:169521/145‑1 (MQ=255)      gATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGgtggtg                                                                                                                                                       >  1:1328051/1‑150 (MQ=255)      gATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGATCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGgatg                                                                                                                                                          >  1:293064/1‑147 (MQ=255)           gAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACGGAACTGTATCGCCACCTGGGCGCGGATACCGACGTGCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGaa                                                                                                                                                     <  1:812231/147‑1 (MQ=255)            aaaGGAAAAGGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGc                                                                                                                                                <  2:743234/151‑1 (MQ=255)            aaaGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGaaaaa                                                                                                                                                  >  2:76478/1‑149 (MQ=255)               ggAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGctct                                                                                                                                             <  1:1333524/151‑1 (MQ=255)                  aaaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGAAAAAACTCTcc                                                                                                                                           >  2:52777/1‑150 (MQ=255)                  aaaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGCTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGCTGAGGataaaagtcccc                                                                                                                                           >  2:11346/1‑143 (MQ=255)                   aaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCcaa                                                                                                                                         <  2:538437/151‑1 (MQ=255)                        gAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGATGAAAAAACTCTCcaacaa                                                                                                                                      >  2:1067121/1‑149 (MQ=255)                        gAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATa                                                                                                                                    <  2:73270/151‑1 (MQ=255)                         aaGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATAc                                                                                                                                   <  2:92128/151‑1 (MQ=255)                              gAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCt                                                                                                                              <  1:996083/151‑1 (MQ=255)                               aaGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTg                                                                                                                             <  2:140343/151‑1 (MQ=255)                                 gTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAGGGAGGAAAAACACTCCCACCCAACCGCCTTC‑g                                                                                                                          >  2:649585/1‑151 (MQ=255)                                  tGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGAAGGAAAAAAACTCccacaacaaccc                                                                                                                                >  2:291652/1‑139 (MQ=255)                                      gCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑c                                                                                                                   <  1:356010/151‑1 (MQ=255)                                          ttttGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑cc                                                                                                              <  2:129416/150‑1 (MQ=255)                                               ccAGGCGCGGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAg                                                                                                       <  2:1304172/151‑1 (MQ=255)                                                 aGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTTGAAAAAATCTCCCCCAAAATCCCCCGC‑GTG‑TT‑C‑T‑GCGGT‑GAggg                                                                                                     >  2:390029/1‑151 (MQ=255)                                                  ggCGCTGATGACTGAACTGTATCGCCCCCTGGGCGCGGATCCCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAgg                                                                                                      <  2:990265/149‑1 (MQ=255)                                                   gcgcTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGAGGGGGTGTAAAAGCACTGCAACAAAACACCCTGC‑GTC‑GC‑C‑A‑CCGGG‑ACGGGcg                                                                                                   >  1:572579/1‑150 (MQ=255)                                                    cgcTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTGAAAAAAAACTCCCACAAAAACGCCCCC‑GTC‑ttttcccggaaagg                                                                                                         >  2:1392560/1‑136 (MQ=255)                                                    cgcTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCt                                                                                                  <  2:601251/151‑1 (MQ=255)                                                       tgatgaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGATGAAAAAAATACCCAACAAAAACACCCGC‑GTC‑CT‑T‑A‑CCGGC‑AAGGGCCctt                                                                                                >  2:1046051/1‑150 (MQ=255)                                                          tgaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTGAAAAAAAACTCCCACAAAAACCCCCCC‑GTCTTT‑C‑C‑CCGGA‑AAGGCCCtttttt                                                                                             >  2:1200857/1‑148 (MQ=255)                                                            aCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGAGGTGAAAAAAATCTCCCACAAAACCGCCCGC‑GGC‑TT‑T‑A‑CCCGG‑AAGGGGCTTTTAtttt                                                                                          >  2:299570/1‑150 (MQ=255)                                                               gAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTggc                                                                                        <  2:1208101/150‑1 (MQ=255)                                                                 aCTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTGAAAAAAATCTCCCACAAAACCCCCTGC‑GTC‑TT‑C‑C‑CCCGT‑AAGGGGCTTTCAttttgggg                                                                                      >  2:15033/1‑145 (MQ=255)                                                                    gTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAg                                                                                   <  2:1263619/150‑1 (MQ=255)                                                                         gCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGACGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTAt                                                                             >  1:189943/1‑151 (MQ=255)                                                                            aCCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGATGAAAAAACACTCCCACAATAACCCCCGC‑GGC‑TT‑T‑C‑CCCGT‑AAGGGCCTTTTATTTTGCGGCAGGCTtttt                                                                            >  2:810375/1‑149 (MQ=255)                                                                            aCCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCg                                                                          <  2:1034158/151‑1 (MQ=255)                                                                                gggCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGccc                                                                       <  2:201693/150‑1 (MQ=255)                                                                                 ggCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCgg                                                                     <  1:702004/151‑1 (MQ=255)                                                                                     cgGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTGAAAAAAATCTCCCACAAAACCCCCCGC‑GGC‑TT‑C‑C‑CCGGG‑AAGGGGCCTTTATTTTGCGGCCGTCTTTTTCTccccg                                                                     >  2:285847/1‑145 (MQ=255)                                                                                          aCCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGATGAAAAAAATCTCCCACAATAACCCCTGC‑GCC‑TT‑C‑A‑CCGGG‑AAGGGGCCTTTTTTTTGCGGCGGTCTTTTTTTCCCCGGAGCTAcc                                                             >  2:1273990/1‑150 (MQ=255)                                                                                            cGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGc                                                          <  2:1154616/151‑1 (MQ=255)                                                                                              aCGTTCCGGCGGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCt                                                         <  2:632549/150‑1 (MQ=255)                                                                                                gTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACg                                                      <  2:1206475/151‑1 (MQ=255)                                                                                                   ccGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGATGAAAAAAAACCCCAACAAAAACCCCCGC‑GTC‑TT‑C‑T‑CCCGTAAAGGGCCTTT‑ATTTTGCGGCAGACTTTTTTTCCCGGCAGCGAcagcgtcagggg                                                   >  1:358504/1‑140 (MQ=255)                                                                                                      gCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTgg                                                <  2:705965/151‑1 (MQ=255)                                                                                                        aGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTGAAAAAAGTCGCCCACAATAACGCCCGC‑GTC‑TT‑C‑A‑CCCGT‑AAGGGGCTTTTATTTGGCGGCAGTCTTATTTGCCCCGGAGGCACCGGcggccggccg                                                 >  1:512370/1‑139 (MQ=255)                                                                                                            gATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGAAAAAAATCTCCCACCAAAACGCCCGC‑GCC‑TT‑T‑ACCCGGG‑AAGGCCCTTT‑ATTTGTGGGCGAGCCTATTTTCCCGGGAGGTACCCGCCACGGTCTGGGTtttt                                          >  1:1108713/1‑151 (MQ=255)                                                                                                             atatCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGATGGAAAAAATCTCCAACAATAACCCCCCC‑GGC‑TTCC‑C‑CCGGA‑AAGGGCC‑TTCATTTTTCGGCAGTCTTATTTTCCCCGGAGATACCAGCTACGGGCTGGTtttttt                                         >  2:1253970/1‑151 (MQ=255)                                                                                                                  ggggTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTca                                       <  1:1309253/148‑1 (MQ=255)                                                                                                                  ggggTTGGTGGGCGGGAGGTCGGCGTTATGGCGGGGGTGATGAAAAAGCTCTCCAACAATGCCGCCGGC‑GTC‑TT‑C‑G‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCGCGGAAGCTGCCGGCTGCGGTCTGGTTTATTTCACAg                                    <  1:643350/151‑1 (MQ=255)                                                                                                                    ggTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGTGGATGAAAAAGCCCTCCAAAAAAAACGCCCGCGGCC‑TT‑C‑C‑CCGG‑‑AATGGCCTTTCATTTGGGGGGGGTGTTTTTCCCCCCGAAGGTACCGGCCACGGTCTGGGtttttt                                         >  2:458349/1‑144 (MQ=255)                                                                                                                       tggtgtTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGc                                <  1:431877/150‑1 (MQ=255)                                                                                                                        ggttgtCGTTATGTCCGCTTAACGGCGGGGATAATGAAAAAGCTCTCCAACCGATCCGCCTGC‑GTC‑CT‑A‑A‑CGGGT‑AAGGGCCTTTCATTTGGAGGCAGTCTTATTCGCCCGGAAGCTACCGGCTGCGGTCTAGTTTAGTTCACAGAAGCaa                              <  1:194464/151‑1 (MQ=255)                                                                                                                        ggtggtCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCa                               <  1:589837/150‑1 (MQ=255)                                                                                                                         gtggtCGTAAAGTGGGCTTTATGGCGGGGATGATGAAAAAGTTCTCCAACAATCCCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAAt                             <  2:1050431/151‑1 (MQ=255)                                                                                                                                tGAAGTCGGCTTTATGGCGGGGGTGATTAAAAAAATCTCCAAAAATAACCCCCGC‑GGC‑TT‑CAA‑CCGGG‑AGGGGCTTTTC‑TTTTGCGGCAGGCTTTTTTTCCCGGGAGCTAACCGCTACGGGCTGGGTTTTTTTACAGAAGCAAAGCTaaa                       >  2:662576/1‑150 (MQ=255)                                                                                                                                tGAAGTCGGCTTTATGGCGGGGGTGATGAAAAAAATCTCCCACAATAACCCCCGC‑GTC‑TT‑C‑A‑CCGGG‑AAGGGGCTTTCATTTTGGGGCAGGCTTATTTTCCCGGAAGGTACCCGCTACGGGCTGGGTTTTTTTACAGAAGCAAAGCTaaa                       >  2:432881/1‑150 (MQ=255)                                                                                                                                  aaGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACg                    <  2:1087893/151‑1 (MQ=255)                                                                                                                                         cTTTATGGCGGGGATGATGAAAAAGCTTTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTGCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACTCCa                 <  1:1175549/147‑1 (MQ=255)                                                                                                                                                gCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGtt      <  1:76465/151‑1 (MQ=255)                                                                                                                                                gCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGtt      <  1:1317656/151‑1 (MQ=255)                                                                                                                                                  ggggATGATAAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGtttt    <  1:646179/151‑1 (MQ=255)                                                                                                                                                      atgatgAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGTTTTGa  >  2:1351674/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                           ATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TT‑C‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGTTTTGA  >  NC_000913/1842765‑1843062

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