Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,636,385 A→T T127S (ACG→TCG)  creB → DNA‑binding transcriptional regulator CreB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,636,3850AT94.8% 147.5 / ‑0.6 58T127S (ACG→TCG) creBDNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (28/27);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.18e-01

CAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCG  >  NC_000913/4636138‑4636614
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
ccagcaggctcccGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGCAGTTCTCGTCGc                                                                                                                                                                                                                                     >  1:155920/11‑251 (MQ=255)
   gcaggctcccGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGt                                                                                                                                                                                                                                  >  1:46505/8‑251 (MQ=255)
             gACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGCCTATGTGGCTAAACCGATTTCACCCCGCGAAGTGTGCGCCAGGGGGCGCACCTTACTGCGTCGGGCGCAGAAGGTCTCGTCGCCGTATCCCGACAt                                                                                                                                                                                                                        >  2:155610/1‑251 (MQ=255)
             gACGTCATGATTCTCGAGGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCAt                                                                                                                                                                                                                        <  2:394211/251‑1 (MQ=255)
                gTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCg                                                                                                                                                                                                                     >  1:172694/1‑251 (MQ=255)
                     gatgctcGATGTTGGTTTTCCTGATATCAGCGGCTTTTAATTATGCCGCCAGTTACTGGCGCTCCAGCCGGCGTTACCTGTACGGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTAtt                                                                                                                                                                                                                 <  2:291267/246‑1 (MQ=255)
                     gATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTAtt                                                                                                                                                                                                                 <  1:80121/250‑1 (MQ=255)
                     gATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTg                                                                                                                                                                                                                >  1:360972/1‑251 (MQ=255)
                           ctatgttgtgttgccttatatCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCACCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTggtc                                                                                                                                                                                                             <  2:326079/231‑3 (MQ=255)
                           cGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGAc                                                                                                                                                                                                             <  1:16259/248‑1 (MQ=255)
                            gATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCAttt                                                                                                                                                                                                         >  1:48596/1‑251 (MQ=255)
                            gATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGTTCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACAttt                                                                                                                                                                                                         >  1:73296/1‑251 (MQ=255)
                               gTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTttg                                                                                                                                                                                                       <  1:12032/250‑1 (MQ=255)
                                     ttGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGGAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGGCATTttgaatt‑ga                                                                                                                                                                                                >  2:11497/2‑250 (MQ=255)
                                          tgATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGACGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGACGGGTGAAGAAGTTCTCGTCGCCGTATCCCGTCCTCCGTATTGGTCATTTTGAATT‑GAATGAAc                                                                                                                                                                                          >  2:403328/2‑251 (MQ=255)
                                                        ttGAATTAGGCCGCCAGTTGCTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGa                                                                                                                                                                              <  2:447001/249‑1 (MQ=255)
                                                          gAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCGCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCa                                                                                                                                                                           <  2:479954/250‑1 (MQ=255)
                                                           aaTTATGCCGCCTGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGc                                                                                                                                                                         <  2:158594/251‑1 (MQ=255)
                                                            attatgcggccAGTTACTGGCGCTCCATCATGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGCAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGAGGCGCACCTTACTGCGTCGGGTGTCGAAGTTCTCGTCGCCGCCTCCCGTCCTCCGTATTGGTCATTTTTAATT‑GAATGAACCAGCGgcgccgtacagct                                                                                                                                                                        >  2:354579/9‑242 (MQ=255)
                                                              tatgCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTAATGTTCCTGACGGTCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGTAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGAGCGCACCTTACTGAGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGGCATTTTGAGTT‑GAATGAACCCGCGGCGCAGTTCAGCTgg                                                                                                                                                                      >  2:90798/3‑251 (MQ=255)
                                                                      cAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGGGAAGAAGTTCTCGTCGCCGGCTCCCGTCATCCGGATTGGTCATTTTGAATT‑GAGTGAACCCGCGGCGCCGCTCAGCTGGTTTGacac                                                                                                                                                              >  2:309595/1‑251 (MQ=255)
                                                                           aCTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTGGACACGCCa                                                                                                                                                          <  2:284138/250‑1 (MQ=255)
                                                                           aCTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCGTACGGCCCGAAGTGAAGAGGTCGATCGCCGGCTTGGGCTGGAAATTGGTG‑CTGACGATTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCAt                                                                                                                                                         <  2:264949/251‑1 (MQ=255)
                                                                            cTGGCGATCCATCCGGCGTTTCCTGTACTGTTCCTGACGGCCCGAAGTGACGAGGTCGATCGCATGCTTGGGCTGTAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCAACCCGCGAAGTGAGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCCTTTTGAATT‑GAATGAACCCGGGGCGACGATCAGCTTGTTTGACACGCCAtt                                                                                                                                                        >  2:311463/1‑251 (MQ=255)
                                                                              ggCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGCGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTa                                                                                                                                                       >  2:329550/1‑250 (MQ=255)
                                                                                cgcTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATt‑gaatga                                                                                                                                                                                            >  2:82883/1‑211 (MQ=255)
                                                                                cgcTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATt‑gaatga                                                                                                                                                                                            <  1:82883/211‑1 (MQ=255)
                                                                                cgcTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTaacg                                                                                                                                                    <  1:336927/251‑4 (MQ=255)
                                                                                    ccttCCGGCTTTACCTGTACTTTTCCTGACGGCCCGATGTGAAGAGGTCGACCGCCTGCTTGGGCTGGAAATTGGGG‑CGGACGTCTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGTTGCGCCCCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCCCCCGTCATCCGTATTGGTAATTTTGAATT‑GAATGACCCCGCGGCGAAGATCAGCTGGTTTGACACGCCATTAACGCTg                                                                                                                                                 <  1:265946/247‑1 (MQ=255)
                                                                                        ccgtcgttactggttcTGTACCTGACGCCCTGATGGGAAGAGGTCGAGCGCCTGCTTGGGCTGGAAAT‑GGTGGCTGACGACTTTGTGGCTAACCCGTTTTCGCCCGGCGAAGTGTGCGCCGGGGTTCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTAGTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGACCAGCGGGTTTGACACGCCATTAACGCTGAcgcg                                                                                                                                            <  2:350709/236‑4 (MQ=255)
                                                                                                ccttGTACTGTTCCTGACGGCCCGAAGGGAGGAGGTCGATCGCCTGCTTGGGCTGGAAAGTGGTG‑CTCACGCATATGTGGCTAAACCTTTTTCCCCCCGCTAAGTTTGCCCCAGGTGGCGCACTTTACTGCGTCGGTTGAAGAATTTCTCGTCGCCGTCTCCCGTCATCCTTTTTGTTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGt                                                                                                                                    <  2:342039/248‑1 (MQ=255)
                                                                                                 cctgtatggTTCCTGTCGGCCCGAAGGGAAGAGGCCGAGCGCCTGCTTGGGTTGGAAATGGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGAGAAGTGTGCGCCAGGGTGCGCCCCTTACTGCGGCGGGTGAAGAAGTTCTCGTCGCCGTCCCCCTTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGtt                                                                                                                                   <  2:429059/243‑1 (MQ=255)
                                                                                                     tACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCCTTAACGCTGACGCGGGATGAGttttt                                                                                                                                >  2:494703/1‑250 (MQ=255)
                                                                                                         gtgccggACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCCCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTAttg                                                                                                                            <  1:103301/244‑1 (MQ=255)
                                                                                                         gTTCCTGACGGCCAGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGAAGACTATGTGGCTCAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGGAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCCGACCAGATGGTTTGACACGCCCTTAACGCTGAACCGGTATGAGTTTTTAttgc                                                                                                                           >  2:259027/1‑250 (MQ=255)
                                                                                                               aaCGGCCCGAAGGGACGGGGTCGAGTGCCTGCTTGGGTGGGAAATTGGTG‑CTGACGACTATGTGGCTAACCCGTTTCCACCCCGCGAAGTGGGCTCGAGGTTGCGCACCTTACGTCGTCGGGTGAAGAAGTTCTAGTAGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGTGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGa                                                                                                                     <  2:239618/250‑1 (MQ=255)
                                                                                                                    cccGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTa                                                                                                                <  1:159992/251‑1 (MQ=255)
                                                                                                                     ccGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGCAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTAc                                                                                                               >  1:200455/1‑251 (MQ=255)
                                                                                                                            gAAGAGGGCGATCGCTTGCTGGGGCAGTAAATTGGTG‑CTGACGACTATGTGGCTAAACCTTTTTCACCCCGCGAAGTGTGCTCCAGGGTGCGTATCTTACGGCGTCGGGTAAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTTCTGAAGACGGTACTCAAGTc                                                                                                        <  1:436268/251‑1 (MQ=255)
                                                                                                                             aaGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTcg                                                                                                       <  1:224148/251‑2 (MQ=255)
                                                                                                                                 ggTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑ATGAAGACTATGTGGATAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGtggcc                                                                                                     >  1:435995/1‑245 (MQ=255)
                                                                                                                                   tCGATCGCCTGCTTGGGCTGGAAATTGGGG‑CTGACGACTATGTGGCGACCCCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGACCGAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTAGTCGTTTTGAATT‑GAATGAACCCTCGGCGCAGATACGCTGGTTTGACACGCCATTATAGCTGCCGCGGTATGAGTTTTTATTGCTGAAGCCGTTACTCacatcgccggcc                                                                                                 >  2:170/1‑240 (MQ=255)
                                                                                                                                        cGCCTGCTTGGGCTGGAAATTGGTG‑CTGCCGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCCTTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTAATCAAGTCGCCAGGCCGCGt                                                                                            >  2:492715/1‑251 (MQ=255)
                                                                                                                                          ccTGCTTGGGCTGGAAATTGGTG‑CTGACGTCTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCt                                                                                          <  1:346616/251‑1 (MQ=255)
                                                                                                                                                 gggCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGGCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTccc                                                                                    <  1:492715/250‑1 (MQ=255)
                                                                                                                                                  ggCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGc                                                                                  <  1:374008/251‑1 (MQ=255)
                                                                                                                                                      ggAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCcgcg                                                                                             >  1:138445/1‑236 (MQ=255)
                                                                                                                                                      ggAAATTGGGG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCcgcg                                                                                             <  2:138445/236‑1 (MQ=255)
                                                                                                                                                                               ggCTAAACCGCTTTCACCCCGCGAAGTGTGCGCAAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGACTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTATGGTCCCGCCAGCAACTAATGGATAGCATGTGGGa                                                        >  1:181492/1‑249 (MQ=255)
                                                                                                                                                                                         ttttCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGACTGGTCCCGCCCGCAACTAATGGATAGCGGGTGGGAAGATgcgccgg                                             >  2:413571/1‑247 (MQ=255)
                                                                                                                                                                                             cACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACAcc                                        >  1:167717/1‑251 (MQ=255)
                                                                                                                                                                                                         ggggcgcCAGGGTGCGCACTTTACTGCGTCGGGTGAAGAATTTCTCGTCGCCTTCCCCCGTCATCCGTATTGGTCATTTTGAATT‑AAATGACCCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCAcc                            <  1:262215/248‑1 (MQ=255)
                                                                                                                                                                                                                 aGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGACGTCTCCCGTCATACGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATAc                    >  1:197167/1‑251 (MQ=255)
                                                                                                                                                                                                                     tGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCTTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCACGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGAGCAGGACACCCACGAACGCACCGTTGATCCCcac                >  2:419354/1‑251 (MQ=255)
                                                                                                                                                                                                                     tGCGCACCTTACTGCGTAGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACCCCTACGAACGCACCGTTGATACCcac                >  2:326222/1‑251 (MQ=255)
                                                                                                                                                                                                                               aCTGCGTCGGGTGAAGAAATTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAAcct      >  2:167476/1‑249 (MQ=255)
                                                                                                                                                                                                                                  gCGTGGGGGGAAGAAGTTCTCTTCGCCGCCCCCCTTCCTCCGTTTTGGTCATTT‑GAATTTGAAGGACCCCGCGGCGCAGACCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACgttg    <  1:394652/250‑4 (MQ=255)
                                                                                                                                                                                                                                   cGTCGGGTGAAGAAGTTCTCGACGCAGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGATGGATTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTAGCCAGTCCGCGTCTGGTCCCGCCAGCAACTAATGGGTACCGGGGGGGAAGATGCGAAGGACCCCTAAGATAGCACAGTTGATACCCACATTAAAACGTTGCg  >  2:387618/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
CAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCG  >  NC_000913/4636138‑4636614

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