Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,509,586	0	A	C	57.1%	32.8 / 61.7	63	pseudogene (108/315 nt)	insO	IS911B regulator fragment
Reads supporting (aligned to +/- strand):  ref base A (11/16);  new base C (17/19);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCCGTTGTTCTCGACCTGTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAACCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCC  >  NC_000913/4509344‑4509828                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    cgtcGTTGTTCTCGACCTGTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACCGCCCCATGGAgcg                                                                                                                                                                                                                                             <  2:21224/248‑1 (MQ=255)             cGACCTGTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTCACCATGAAAGCGCTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGGGATGTTCCACAGCGTTCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTCTGAGCCGGCGCGGAACCCGCTGGGATAACTGCCCAGTGGAACCCTTTTTCAgtg                                                                                                                                                                                                                                  >  2:68485/1‑249 (MQ=255)              gACCTGTTCGCAAGAAAACCAGTGGGCTGGGCAATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGGGCGCTTCTTCAGGAGt                                                                                                                                                                                                                                >  2:88312/1‑251 (MQ=255)                ccTGTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCt                                                                                                                                                                                                                              >  1:4968/1‑251 (MQ=255)                   gTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTCACCATGAAAGCGCTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTTGGATAACAGCCCAATGTAACGCTTCTTCAGGAGTCt                                                                                                                                                                                                                              <  2:32630/249‑1 (MQ=255)                    ttCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGAACAGCAGGCTTACCATGAAAGCAATGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTgagg                                                                                                                                                                                                                          >  2:232155/1‑249 (MQ=255)                      cGCAAGAAAACCAGGGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTCACCATGAAAGCGCTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGATCCGGCGCGGAAACTGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTgaag                                                                                                                                                                                                                          <  2:103123/250‑1 (MQ=255)                         aaGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGCCAGCAGGCTCACCATGAAAGCGCTGGAACTGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGCTGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAAATGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACg                                                                                                                                                                                                                      >  2:239222/1‑251 (MQ=255)                           gAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACAATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAt                                                                                                                                                                                                                   >  1:244420/1‑251 (MQ=255)                           gAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTCACCATGAAAGCGCTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGaa                                                                                                                                                                                                                    <  2:85359/251‑1 (MQ=255)                            aaaaCCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑g                                                                                                                                                                                                                 <  1:221013/251‑1 (MQ=255)                                ccAGTGGGCTGGGCCATGTCGTGCTCGCCGGACAGCAGGCTCACCATGAAAGCGCTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATTTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGAGAACAGCCCAATGGAACGCGTCTTCAGGAGTCTGAAGAACGAAT‑ggat                                                                                                                                                                                                              <  2:13354/251‑3 (MQ=255)                                       gctgggcCATGTCGTTCTCGCCGGACAGCAGGCTCACCATGAAAGCGCTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGccggtg                                                                                                                                                                                                       <  1:171423/251‑5 (MQ=255)                                          gggcCAGGTCGTTCGCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCA‑AGCGAGCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGGCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCGTCAGGAGGCTGAAGATCGAAT‑GGGTGCCAGCGAc                                                                                                                                                                                                     <  2:176865/249‑1 (MQ=255)                                                    gTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAg                                                                                                                                                                                         <  2:10557/251‑1 (MQ=255)                                                     ttCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGc                                                                                                                                                                                        >  2:231977/1‑251 (MQ=255)                                                           ccGGACAGCAGGCTCACCATGAAAGCGCTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCCGGTGGTGGGTTACGTAAGCTTCAg                                                                                                                                                                                   <  2:206607/251‑1 (MQ=255)                                                            cGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCg                                                                                                                                                                                 >  2:70770/1‑251 (MQ=255)                                                             ggACAGCAGGCTCACCATGAAAGCGCTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCCGGTGGTGGGTTACGTAAGCTTCAGCg                                                                                                                                                                                 <  2:103841/251‑1 (MQ=255)                                                                       cTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCAc                                                                                                                                                                      >  2:138525/1‑251 (MQ=255)                                                                        ttACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACg                                                                                                                                                                     >  1:203262/1‑251 (MQ=255)                                                                         caCCATGAAAGCGCTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCCGGTGGTGGGTTACGTAAGCTTCAGCGAGGCAGCTCACg                                                                                                                                                                     <  2:106220/250‑1 (MQ=255)                                                                           ccATGAAAGCGCTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCCGGTGGTGGGTTACGTAAGCTTCAGCGAGGCAGCTCACGcc                                                                                                                                                                   <  2:265456/251‑2 (MQ=255)                                                                             aTGAAAGCTCTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCCGGTGGTGGGTTACGTAAGCTTCAGCGAGGCAGCTCACGccat                                                                                                                                                                 <  1:165432/251‑4 (MQ=255)                                                                               gAAAGCTCTGTAAGTGGCATGGTAAACCCGTGGTAAGCCCTTCGGGGTGATGTGCCA‑AGCGAGTAAGGCAGTCATTATACTAGCAGGCAGTTCCGGCAGTTACTGTGGCGATTCCGGATCAGGCAGAGGATTAGTCGGCGTGGAAACGGCTGGGATAACAGCCCCATTGAGCGCTTCTTCAGTGGTCTGAAGAACCAATTGGGT‑CCAGCGACGGGCTATGTAAGCTGCAGCGATTCTGCTTACTCAATAAc                                                                                                                                                              <  2:32337/251‑1 (MQ=255)                                                                                aaaGCGCTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCCGGTGGTGGGTTACGTAAGCTTCAGCGAGGCAGCTCACGCCATAAc                                                                                                                                                              >  1:209623/1‑251 (MQ=255)                                                                                  agcgcTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCCGGTGGTGGGTTACGTAAGCTTCAGCGAGGCAGCTCACGCCATAACgg                                                                                                                                                            >  1:26608/5‑251 (MQ=255)                                                                                       tGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCCGGTGGTGGGTTACGTAAGCTTCAGCGAGGCAGCTCACGCCATAACGGACtct                                                                                                                                                       >  2:4386/1‑249 (MQ=255)                                                                                         gAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGATAACAGCCCAATGGAACGCTTCTGCAGGAGGCTTAAGAACGAAT‑GGATGCCGGTGGTGGGTTACGTAAGCTTCAGCGAGGCAGCTCACGCCATAACGGACTatat                                                                                                                                                     <  1:239222/251‑1 (MQ=255)                                                                                              ggCATGGGAAACCCGGTGTAAGCCCGTCGGGGTGATGTTCCA‑AGAGATCAAGGCAGTCATGATGCGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGGTGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCGTCTTCAGGAGTCTGAAGAACGAAG‑GGGGGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAAATCACTATATCGtggg                                                                                                                                               <  2:39026/251‑4 (MQ=255)                                                                                                    ggAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATCACAGCCCCATGGCGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATAtt                                                                                                                                          >  1:141962/1‑250 (MQ=255)                                                                                                      aaaCCCGTGGTAAGCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCCGGTGGTGGGTTACGTAAGCTTCAGCGAGGCAGCTCATGCCATAACGGACTATATCGTTGGATATTAc                                                                                                                                        <  1:46698/251‑1 (MQ=255)                                                                                                        aaccgaggTAAGCCCGTCTGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATGATACGAGCAGTCAGTTCCTGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTAGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCATGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGTCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTAGATCTTTGGATATTACAgc                                                                                                                                     <  2:177056/245‑1 (MQ=255)                                                                                                              ggTAAGCCCGTCTGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTTGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATTTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTa                                                                                                                               <  1:47511/251‑1 (MQ=255)                                                                                                                 aaGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGa                                                                                                                            <  1:252830/251‑1 (MQ=255)                                                                                                                  aGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGAc                                                                                                                           <  1:225923/251‑1 (MQ=255)                                                                                                                  aGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGAc                                                                                                                           <  2:4968/251‑1 (MQ=255)                                                                                                                   gCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCCGGTGGTGGGTTACGTAAGCTTCAGCGAGGCAGCTCACGCCATAACGGACTATATCGTTGGATATTACAGCGCACTAAGAc                                                                                                                           <  1:243870/251‑1 (MQ=255)                                                                                                                       gTCGGGGTGATGTTCCA‑AGCGATCCAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCTCACTAAGCCCGCAc                                                                                                                      <  1:67350/251‑1 (MQ=255)                                                                                                                             gTGATGTTCCA‑AGCGATCAAGGCAGTCATTAGACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAatat                                                                                                                <  1:259939/251‑1 (MQ=255)                                                                                                                                 tGTTCCACAGCGAACAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCCAATCCGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGCTCAAAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCCGGTGGTGGGCTACGTAAGCTTCAGCGAGGCAGCTAACGCCATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAatatcac                                                                                                             >  2:21155/1‑248 (MQ=255)                                                                                                                                      cA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTggg                                                                                                       >  1:123756/1‑251 (MQ=255)                                                                                                                                         aGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCCGGTGGTGGGTTACGTAAGCTTCAGCGAGGCAGCTCACGCCATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAACGGTGGGtt                                                                                                     >  2:239449/1‑251 (MQ=255)                                                                                                                                               cagtgCAGTCATTATACGAGCAGGCAGTTTCGGCAGTTATTGTGGCGATCCCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGATAACCGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCCGGTGGTGGGTTACGTAAGCTTCAGCGAGGCAGCTCACGCCATAACGGACTATATCGTTGGATATTACAGCGCACTAAGAACGCACGAATATAACGGTGGGTTaccccc                                                                                               >  2:96026/5‑248 (MQ=255)                                                                                                                                               cAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTaccacc                                                                                               <  1:138525/251‑1 (MQ=255)                                                                                                                                                aaGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGCTCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGCTAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTcccacc                                                                                               >  1:153231/1‑250 (MQ=255)                                                                                                                                                 aGGCAGTCATTATACGAGCAGGCAGTGCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCGTCAGCGATGCAGCTCACGCAATAACGGACTATATCGGTGGATATTACAGCGCACTAAGACCGCACGAGTATAATGGTGGGTTACCACCaa                                                                                             <  2:244420/251‑1 (MQ=255)                                                                                                                                                        cATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACGGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCCGGTGGTGGGTTACGTAAGCTTCAGCGAGGCAGCTCACGCCATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAACGGTGGGTTACCCCCAAACGAATc                                                                                      <  2:26608/251‑1 (MQ=255)                                                                                                                                                                gAGCAGGCAGTTACGGCAGTTATTGTGGCGATACCAAATCAGCCAGCGGCTGAGCCGCCGAGGAAACTGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCCGGTGGGGGGTTACGTAAGCTTCAGCGAGGACGCTCACGCCATAACGGACTATATCGTTGGATCTTACAGCGCACTAAGACCGCACGCATATAACGGGGGGCTAACCCCAAACGAATCGGaaaatc                                                                              >  2:53889/1‑249 (MQ=255)                                                                                                                                                                     ggcagTTCAGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCAGGAGGTGGGTTACGTAAGCTTCAGCGAGGCAGCTCACGCCATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAACGGTGGGTTACCCCCAAACGAATCGGaaaatagaaac                                                                         >  2:202802/1‑244 (MQ=255)                                                                                                                                                                                                   ggATCAGGCAGAGTATGAGTCGGCGGGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAACCGATACTGGAAAAACTCTAACGCGGAGGCCAGtttt                                           <  2:257816/251‑1 (MQ=255)                                                                                                                                                                                                     aTCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGATAACAGCCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCCGGTGGTGGGTTACGTAAGCTTCAGCGAGGCAGCTCACGCCATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAACGGTGGGTTACCCCCAAACGAATCGGAAAATCGATACTGGAAAAACTCTAACGCGGTGGCCAGtttttg                                         >  2:115186/1‑249 (MQ=255)                                                                                                                                                                                                        aGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTCCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAACCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTg                                      <  1:74411/251‑1 (MQ=255)                                                                                                                                                                                                          ccAGAGTATGAGCCGGCGCGGAAACTGCTGGGATAACAGCCCAATGGACCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCCGGTGGTGGGTTACGTAAGCTTCAGCGAGGCAGCTCACGCCATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAACGGTGGGTTACCCCCAAACGAATCGGAAAATCGATACTGGAAAAACTCTAACGCGGTGGCCAGTTTTTGTTGAc                                    <  2:105451/250‑1 (MQ=255)                                                                                                                                                                                                            agagTATGAGTCGGCGCGGAAACTGCTGGGATAACAGTCCAATGGAACGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGATGCCGGTGGTGGGTTACGTAAGCTTCAGCGAGGCAGCTCACGCCATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAACGGTGGGTTACCCCCAAACGAATCGGAAAATCGATACTGGAAAAACTCTAACGCGGTGGCCAGTTTTTGTTGACca                                  <  1:212853/251‑1 (MQ=255)                                                                                                                                                                                                              agTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGACT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAACCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACcac                                 >  2:59443/1‑250 (MQ=255)                                                                                                                                                                                                                  tGAGTCGGCGTGGAAACTGCTGAGATAACAGCCCCCTGGAGCGCTTCTTCAGGAGTCGGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTAAGCGATGCAGCTCACGCAATAAAGGAATATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAACCGATACTGGAAAAACTCTAACGCGGATGCCAGTTTTAGTTGACCacaac                              >  1:12382/1‑249 (MQ=255)                                                                                                                                                                                                                    atccGGCGTGGAATCTGCTGTGAGATCCGCCCCATGGAGCGCTGCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTTAGGCCGCACTCATATAATGGTTGGTTACCACCAAACGAATCAGAAATCCTATACTTGAAACACTCTAGCGCGGACGCCAGTTTTAGTTCATCACCACAGAc                          <  2:45710/248‑1 (MQ=255)                                                                                                                                                                                                                        ggCGTGGAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCTATGCAGCTCACGCAATAACGGACTATATCGTTGGATAGTACAGCGCACTAAGACCGCACGAATATAATGGTTGGTTACCACCAAACGAATCAGAAAACCGATTTTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTAcc                      <  1:263841/251‑1 (MQ=255)                                                                                                                                                                                                                              gAAACTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAACCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAgg                >  1:247187/1‑251 (MQ=255)                                                                                                                                                                                                                               aaaCTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAACCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAggg               >  1:169332/1‑251 (MQ=255)                                                                                                                                                                                                                                 aCTGCTGGGATAACAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAACCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAggg               >  1:61818/1‑249 (MQ=255)                                                                                                                                                                                                                                            aaCAGCCCCATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAACCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGcc  >  1:143239/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    CGCCGTTGTTCTCGACCTGTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAACCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCC  >  NC_000913/4509344‑4509828

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