Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,080,837	0	G	T	39.4%	105.6 / 32.6	71	intergenic (+24/‑519)	putP/efeU	proline:sodium symporter/ferrous iron permease (pseudogene)
Reads supporting (aligned to +/- strand):  ref base G (9/34);  new base T (28/0);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 2.72e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCA  >  NC_000913/1080693‑1080983                                                                                                                                                 |                                                                                                                                                      tATTGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGtttttttttg                                                                                                                                                  >  1:534244/1‑148 (MQ=255)  aTTGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTtttgttt                                                                                                                                                    >  2:959546/1‑147 (MQ=255)     gggATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTca                                                                                                                                             <  1:916701/151‑1 (MQ=255)           gTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGGTTCagcagc                                                                                                                                        >  2:668607/1‑150 (MQ=255)              gtgtTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGtttttttttggttttgcaaccgg                                                                                                                                     >  1:409910/1‑135 (MQ=255)              gtgtTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGt                                                                                                                                    <  1:182414/151‑1 (MQ=255)                gtTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGtt                                                                                                                                   <  1:118955/150‑1 (MQ=255)                gtTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGtt                                                                                                                                   <  1:727080/150‑1 (MQ=255)                gtTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGtt                                                                                                                                   <  2:549800/150‑1 (MQ=255)                gtTTAGTTTGCTGGGTAAAGCGCCGGCAGCGGTGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGt                                                                                                                                    <  2:28458/149‑1 (MQ=255)                 tttAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCAGCggggt                                                                                                                                   >  2:638014/1‑147 (MQ=255)                   tAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTgcg                                                                                                                                <  1:696305/150‑1 (MQ=255)                     gTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑c                                                                                                                           <  1:230144/151‑1 (MQ=255)                      tttGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGGTTTGTTTTGCGTCAATAGCAGGGTGAG‑CT‑cc                                                                                                                          >  2:92486/1‑151 (MQ=255)                      tttGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGGTTTGGTTTGCTTTATCAGACGGGTGGG‑CT‑cc                                                                                                                          >  2:95835/1‑151 (MQ=255)                         gCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCt                                                                                                                        <  1:38559/150‑1 (MQ=255)                          cTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCAGCGGGGTGCG‑Cc‑ccc                                                                                                                         >  1:730253/1‑148 (MQ=255)                            gggTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCAGCGGGGTGCG‑CT‑CCCCTaa                                                                                                                     >  2:822879/1‑150 (MQ=255)                            gggTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTaa                                                                                                                     <  1:309956/150‑1 (MQ=255)                            gggTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAAt                                                                                                                    <  2:349044/151‑1 (MQ=255)                            gggTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAAt                                                                                                                    <  2:678726/151‑1 (MQ=255)                            gggTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAAt                                                                                                                    <  1:984243/151‑1 (MQ=255)                             ggTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGGTGCG‑CC‑CCCttta                                                                                                                     >  2:1176273/1‑147 (MQ=255)                              gTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAAtgt                                                                                                                  <  2:1202418/151‑1 (MQ=255)                              gTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAAtgt                                                                                                                  <  1:882750/151‑1 (MQ=255)                                aaaGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAAtgtg                                                                                                                 <  2:121336/150‑1 (MQ=255)                                aaaGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑c                                                                                                               <  1:491604/151‑1 (MQ=255)                                 aaGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCCCCAGCGGGGTGCG‑CT‑CCCTTTAtgtg                                                                                                                 >  2:1226554/1‑149 (MQ=255)                                  aGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGCTTCAGCAGCGGGTTGCG‑CT‑CCCCTAATGTG‑cc                                                                                                              >  2:863166/1‑150 (MQ=255)                                    cgcCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCAGCGGGTTGCG‑CT‑CCCCTTATGTG‑CCTc                                                                                                            >  2:68416/1‑150 (MQ=255)                                    cgcCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGGTTTTTTTGGGTTCAGCAGCGGGGTGCG‑CC‑CCCCTAATGTG‑CCTc                                                                                                            >  1:733657/1‑150 (MQ=255)                                    cgcCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTTGTTGGCTTCTTCAACGGGGTGCG‑GT‑CCCCTAATATG‑cccc                                                                                                            >  1:187107/1‑148 (MQ=255)                                     gcCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTTCGGAACGGGTTTCG‑CC‑CCCCTAATGTG‑cc                                                                                                              >  1:546107/1‑147 (MQ=255)                                        gTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCACCAGCGGGTTGCG‑CT‑CCCTTAATTTG‑CCTCGccct                                                                                                       >  2:235765/1‑149 (MQ=255)                                          cAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCatat                                                                                                     <  2:459915/151‑1 (MQ=255)                                          cAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCata                                                                                                      <  2:616795/150‑1 (MQ=255)                                            gcggcgATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCCCCAGCGGGGTGGG‑CT‑CCCTTCATGTG‑CCTCGCCCtata                                                                                                    >  2:129727/1‑150 (MQ=255)                                                cgATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTg                                                                                               <  1:801099/151‑1 (MQ=255)                                                 gATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGa                                                                                              <  1:152629/151‑1 (MQ=255)                                                  aTGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGaa                                                                                             <  1:835915/151‑1 (MQ=255)                                                             cTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCAGCGGGTTGGG‑CT‑CCCTTAATTTG‑CCTCCCCCATTAAATTTAATGGGGCAggg                                                                                   >  2:1189138/1‑150 (MQ=255)                                                               ttGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGAATGGTGCAGGGa                                                                                  <  2:447496/149‑1 (MQ=255)                                                                tGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTCAGCAGCGGGGTGCG‑CT‑CCCCTAATGTG‑CCTCCCCCTATAAAATGAATGGTGCAggg                                                                                   >  2:457279/1‑147 (MQ=255)                                                                tGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGCTTTTTTTGTCCTCACGAGCGGGGTTCG‑CC‑CCCTTTATATG‑CCTCCCCCCATAAATTGAATGGGGGAGGGAacg                                                                               >  1:465421/1‑151 (MQ=255)                                                                 gCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTTAGCAGCGGGGTGCG‑CT‑CCCCTTATGTG‑GCTCGCCCTATAAATTGAATGGTGCAggg                                                                                   >  2:1190695/1‑146 (MQ=255)                                                                  ccGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCCCCAGCCGGTTGTG‑CC‑CCCCTTATGTGCCCTCCCCTTTTAAATT‑AATgggggcggg                                                                                   >  1:6745/1‑138 (MQ=255)                                                                   cGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGTTTAGCAGCGGGGTTCG‑CT‑CCCCTTATGTG‑CCTCGCCCTATAAATTTAATGGGGGAGGGGgggcg                                                                             >  2:519321/1‑150 (MQ=255)                                                                       gCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCAGCGGGGTGCC‑CC‑CCCCTAATGTG‑CCTCCGCCTATAAAATTAATGGGGCGGGGAGCGCCCAggg                                                                        >  1:145755/1‑151 (MQ=255)                                                                       gCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAggg                                                                        <  2:1151937/151‑1 (MQ=255)                                                                       gCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAggg                                                                        <  1:448750/151‑1 (MQ=255)                                                                          gATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGggc                                                                     <  2:8950/151‑1 (MQ=255)                                                                           aTGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCAGCGGGTTGCG‑CC‑CCCTTAATGGG‑CCTCGCCATATAAAATGAATGGGGCAGGGAGCGCGCAGGGggcg                                                                    >  2:904608/1‑151 (MQ=255)                                                                           aTGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGggcg                                                                    <  1:208635/151‑1 (MQ=255)                                                                            tGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAACAACGGGGTTCG‑CC‑CCCCTTATTTG‑CCCCGCCCTTTAAATTGAAAGGTGGAGGGAGCGCGCAGGGggcgg                                                                   >  2:879106/1‑151 (MQ=255)                                                                                 aCTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGCTTCAGCAGCGGGGTGCG‑CT‑CCCCTAATGTG‑CCTCGCCCTATAAATTGAATGGTGCAGGGGGCGCGCAgggggggggc                                                                 >  1:1008080/1‑146 (MQ=255)                                                                                    aTCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATcgc                                                           <  1:1103645/151‑1 (MQ=255)                                                                                      cATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATcgcc                                                          <  2:1215893/150‑1 (MQ=255)                                                                                       aTTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCCCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGGGGCCAATcgccgc                                                        >  2:759402/1‑151 (MQ=255)                                                                                          cGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATcgccgtcg                                                      <  1:254793/150‑1 (MQ=255)                                                                                             cTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTTATGTG‑CCTCCCCCTATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCcccccc                                                  >  2:285364/1‑151 (MQ=255)                                                                                               ccgccgTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCct                                                <  2:100262/151‑1 (MQ=255)                                                                                                cgccgTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGGTTTTTTTCGACTCAGCAGCGGGGTGTGCCT‑CCCTT‑ATTCG‑CCTGGCCCGATAAAAAGGAGAGTGCCAGGAACGGCCCGGGGGCGGCCAAGCGGCTCCGcccc                                                  >  1:1236721/1‑148 (MQ=255)                                                                                                            tGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGcc                                    <  2:312190/150‑1 (MQ=255)                                                                                                               aGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGGTTCAGCAACGGGGTTCG‑CT‑CCCCTTATTGG‑CCTCGCCCTATAAAATTAATGGTGCAGGGAGCGCGCAGGGGGCGGGCCAACCCCCCCCCCCCCCGCCGTTCCCGCCCTc                                 >  1:580039/1‑150 (MQ=255)                                                                                                                ggAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTGTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑GCTTTCCATAAAAATTTAATGGGGCAGGGAGAGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCCGCTGTCCCGGCCTTCgg                               >  1:548348/1‑151 (MQ=255)                                                                                                                   aaGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTTGCTTCAGCAGCGGGGTGCG‑CT‑CCCCTTATTTG‑CCTCGCCATATAAATTGAATGGGGCAGGGGGCGCGCCAGGGGCGGCCAAACGCCCCCGCCCCCCGCTGTCCCCGGCCTCggggg                            >  2:450405/1‑150 (MQ=255)                                                                                                                   aaGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCCTTAGCAACGCGGTGCG‑CTCCCCTTAATGGG‑CC‑CGCCCTATAAATTTAAAAGGGCAGGGGGCGCGCCCGGGGCGGGCAAAAGCCCCCGCCCCCTGCCGCCCCGGGCCTCtgggg                            >  1:138692/1‑150 (MQ=255)                                                                                                                       tAAGGGACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCAGCGGGGTGCG‑CT‑CCCCTAATGTG‑CCTCGCCCTATAAATTGAATGGTGGAGGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACg                         >  2:243414/1‑150 (MQ=39)                                                                                                                        aaGGGACTTAGCCTGCGGCGGCTTTTTTTTGGTTCTGCAGCGGGGTGCG‑CT‑TCCCTTACGTG‑GCTCCCCATATAAATTGAATGGGGGAGGGGGCGCGCACGGGGGGGCCCATCTCCGGCCCCCCCCTCCGTGCCGGCCTTTgggg                             >  1:424497/1‑145 (MQ=255)                                                                                                                           ggACTTAGCCTGCGGCGGTTTTTTTTGGCTTCAGCAGCGGGGTGCG‑CT‑CCCTTAATGTG‑CCTCGCCCTATAAATTGAATGGGGCAGGGAGCGCGCAGGGGGGGGCCCATCGCCGCCCCCCCCCGCTGTCCCCGCCCTCGGGGGACGCtt                      >  2:671369/1‑149 (MQ=255)                                                                                                                            gACTTAGCCTGCGGCGGTTTTGTTTTGCTTCAGCAACGGGGTGCG‑CT‑CCCCTTATTTG‑CCTCGCCCTATAAAATTAATGGTGCAGGGGGCGCGCAGGGGGCGGCCAAACGCCCCCGCCCCCCGCTGTCCCCGCCTTCGGGGGACGCTTc                     >  2:876159/1‑149 (MQ=255)                                                                                                                                  gCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGAttt             <  1:959494/151‑1 (MQ=255)                                                                                                                                        ggcggTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGc       >  2:193599/1‑151 (MQ=255)                                                                                                                                            gTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGccacc   <  1:806790/151‑1 (MQ=255)                                                                                                                                              tttgtttgGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGTGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGccacca  >  1:21948/1‑150 (MQ=255)                                                                                                                                              tttgtttgGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGccacca  >  2:367443/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                      TATTGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTGTTTGGCTTCAGCAGCGGGTTGCG‑CT‑CCCTTAATGTG‑CCTCGCCATATAAATTGAATGGTGCAGGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCA  >  NC_000913/1080693‑1080983

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