Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,509,5560TC50.0% ‑1.2 / 72.4 62pseudogene (78/315 nt)insOIS911B regulator fragment
Reads supporting (aligned to +/- strand):  ref base T (13/17);  new base C (12/18);  total (25/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATATCTGGACGGGTAAGCGCTGGGCGTACCTCGCCGTTGTTCTCGACCTGTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAACCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAAC  >  NC_000913/4509313‑4509800
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
atatCTGGACGGGTAATCGCTGGTCGTACCTCGCCGTTGTTCGCTACCTGTTCCCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGGCGGACAGCAGGCTTACCATGAAAGCACTGCAAATGACATGGGAAACCCGTGCTAAGCCCGTCTGGGTGATCTTCCA‑ATCGATCAATCCAGTCAGTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTg                                                                                                                                                                                                                                                <  1:237576/251‑1 (MQ=255)
  atCTGGACGGGTAAGCGCTGGGCGTACCTCGCCGTTGTTCTCGACCTGTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTCACCATGAAAGCGCTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCgcg                                                                                                                                                                                                                                                <  1:171885/250‑3 (MQ=255)
        aCGGGTAAGCGCTGGGCGTACCTCGCCGTTGTTCTCGACCTGTGCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTCACCATGAAAGCGCTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTTTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAActg                                                                                                                                                                                                                                         <  1:89148/251‑1 (MQ=255)
          gggTAAGCGCTGGGCGTACCTCGCCGTTGTTCTCGACCTGTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTCACCATGAAAGCGCTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAActgct                                                                                                                                                                                                                                       >  1:32630/1‑251 (MQ=255)
                gcgcTGGGCGTACCTCGCCGTTGTTCTCGACCTGTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTCACCATGAAAGCGCTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATGTTCCACAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTATTGTGGCGATACCAGATCAGACAGAGTATGAGCCGGCGCGGAAACTGCTGGGATa                                                                                                                                                                                                                                 <  1:27843/251‑1 (MQ=255)
                               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)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
ATATCTGGACGGGTAAGCGCTGGGCGTACCTCGCCGTTGTTCTCGACCTGTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCA‑AGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAAT‑GGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAACCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAAC  >  NC_000913/4509313‑4509800

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