Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,636,411 A→T G135G (GGA→GGT)  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,4110AT92.3% 138.0 / 1.6 54G135G (GGA→GGT) creBDNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base T (24/24);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAGCGCTGACTCG‑GTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACC  >  NC_000913/4636165‑4636650
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
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‑GAATGAACCCGCGGCGCAGTTCAGCT‑gg                                                                                                                                                                                                           >  2:90798/3‑251 (MQ=255)
                                           cAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGGGAAGAAGTTCTCGTCGCCGGCTCCCGTCATCCGGATTGGTCATTTTGAATT‑GAGTGAACCCGCGGCGCCGCTCAGCT‑GGTTTGacac                                                                                                                                                                                                   >  2:309595/1‑251 (MQ=255)
                                                aCTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTGGACACGCCa                                                                                                                                                                                               <  2:284138/250‑1 (MQ=255)
                                                aCTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCGTACGGCCCGAAGTGAAGAGGTCGATCGCCGGCTTGGGCTGGAAATTGGTG‑CTGACGATTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCAt                                                                                                                                                                                              <  2:264949/251‑1 (MQ=255)
                                                 cTGGCGATCCATCCGGCGTTTCCTGTACTGTTCCTGACGGCCCGAAGTGACGAGGTCGATCGCATGCTTGGGCTGTAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCAACCCGCGAAGTGAGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCCTTTTGAATT‑GAATGAACCCGGGGCGACGATCAGCT‑TGTTTGACACGCCAtt                                                                                                                                                                                             >  2:311463/1‑251 (MQ=255)
                                                   ggCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGCGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTa                                                                                                                                                                                            >  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‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTaacg                                                                                                                                                                                         <  1:336927/251‑4 (MQ=255)
                                                         ccttCCGGCTTTACCTGTACTTTTCCTGACGGCCCGATGTGAAGAGGTCGACCGCCTGCTTGGGCTGGAAATTGGGG‑CGGACGTCTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGTTGCGCCCCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCCCCCGTCATCCGTATTGGTAATTTTGAATT‑GAATGACCCCGCGGCGAAGATCAGCT‑GGTTTGACACGCCATTAACGCTg                                                                                                                                                                                      <  1:265946/247‑1 (MQ=255)
                                                             ccgtcgttactggttcTGTACCTGACGCCCTGATGGGAAGAGGTCGAGCGCCTGCTTGGGCTGGAAAT‑GGTGGCTGACGACTTTGTGGCTAACCCGTTTTCGCCCGGCGAAGTGTGCGCCGGGGTTCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTAGTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGACCAGCG‑GGTTTGACACGCCATTAACGCTGAcgcg                                                                                                                                                                                 <  2:350709/236‑4 (MQ=255)
                                                                     ccttGTACTGTTCCTGACGGCCCGAAGGGAGGAGGTCGATCGCCTGCTTGGGCTGGAAAGTGGTG‑CTCACGCATATGTGGCTAAACCTTTTTCCCCCCGCTAAGTTTGCCCCAGGTGGCGCACTTTACTGCGTCGGTTGAAGAATTTCTCGTCGCCGTCTCCCGTCATCCTTTTTGTTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGt                                                                                                                                                                        <  2:342039/248‑1 (MQ=255)
                                                                      cctgtatggTTCCTGTCGGCCCGAAGGGAAGAGGCCGAGCGCCTGCTTGGGTTGGAAATGGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGAGAAGTGTGCGCCAGGGTGCGCCCCTTACTGCGGCGGGTGAAGAAGTTCTCGTCGCCGTCCCCCTTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGtt                                                                                                                                                                       <  2:429059/243‑1 (MQ=255)
                                                                          tACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCCTTAACGCTGACGCG‑GGATGAGttttt                                                                                                                                                                    >  2:494703/1‑250 (MQ=255)
                                                                              gtgccggACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCCCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTAttg                                                                                                                                                                <  1:103301/244‑1 (MQ=255)
                                                                              gTTCCTGACGGCCAGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGAAGACTATGTGGCTCAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGGAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCCGACCAGAT‑GGTTTGACACGCCCTTAACGCTGAACCG‑GTATGAGTTTTTAttgc                                                                                                                                                               >  2:259027/1‑250 (MQ=255)
                                                                                    aaCGGCCCGAAGGGACGGGGTCGAGTGCCTGCTTGGGTGGGAAATTGGTG‑CTGACGACTATGTGGCTAACCCGTTTCCACCCCGCGAAGTGGGCTCGAGGTTGCGCACCTTACGTCGTCGGGTGAAGAAGTTCTAGTAGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGTGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGa                                                                                                                                                         <  2:239618/250‑1 (MQ=255)
                                                                                         cccGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTa                                                                                                                                                    <  1:159992/251‑1 (MQ=255)
                                                                                          ccGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGCAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTAc                                                                                                                                                   >  1:200455/1‑251 (MQ=255)
                                                                                                 gAAGAGGGCGATCGCTTGCTGGGGCAGTAAATTGGTG‑CTGACGACTATGTGGCTAAACCTTTTTCACCCCGCGAAGTGTGCTCCAGGGTGCGTATCTTACGGCGTCGGGTAAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTTCTGAAGACGGTACTCAAGTc                                                                                                                                            <  1:436268/251‑1 (MQ=255)
                                                                                                  aaGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTcg                                                                                                                                           <  1:224148/251‑2 (MQ=255)
                                                                                                      ggTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑ATGAAGACTATGTGGATAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGtggcc                                                                                                                                         >  1:435995/1‑245 (MQ=255)
                                                                                                        tCGATCGCCTGCTTGGGCTGGAAATTGGGG‑CTGACGACTATGTGGCGACCCCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGACCGAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTAGTCGTTTTGAATT‑GAATGAACCCTCGGCGCAGATACGCT‑GGTTTGACACGCCATTATAGCTGCCGCG‑GTATGAGTTTTTATTGCTGAAGCCGTTACTCacatcgccggcc                                                                                                                                     >  2:170/1‑240 (MQ=255)
                                                                                                             cGCCTGCTTGGGCTGGAAATTGGTG‑CTGCCGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCCTTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTAATCAAGTCGCCAGGCCGCGt                                                                                                                                >  2:492715/1‑251 (MQ=255)
                                                                                                               ccTGCTTGGGCTGGAAATTGGTG‑CTGACGTCTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCt                                                                                                                              <  1:346616/251‑1 (MQ=255)
                                                                                                                      gggCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGGCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTccc                                                                                                                        <  1:492715/250‑1 (MQ=255)
                                                                                                                       ggCTGGAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGc                                                                                                                      <  1:374008/251‑1 (MQ=255)
                                                                                                                           ggAAATTGGTG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCcgcg                                                                                                                                 >  1:138445/1‑236 (MQ=255)
                                                                                                                           ggAAATTGGGG‑CTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCcgcg                                                                                                                                 <  2:138445/236‑1 (MQ=255)
                                                                                                                                                    ggCTAAACCGCTTTCACCCCGCGAAGTGTGCGCAAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGACTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTATGGTCCCGCCAGCAACTAATGGATAGCATGTGGGa                                                                                            >  1:181492/1‑249 (MQ=255)
                                                                                                                                                              ttttCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGACTGGTCCCGCCCGCAACTAATGGATAGCGGGTGGGAAGATgcgccgg                                                                                 >  2:413571/1‑247 (MQ=255)
                                                                                                                                                                  cACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACAcc                                                                            >  1:167717/1‑251 (MQ=255)
                                                                                                                                                                              ggggcgcCAGGGTGCGCACTTTACTGCGTCGGGTGAAGAATTTCTCGTCGCCTTCCCCCGTCATCCGTATTGGTCATTTTGAATT‑AAATGACCCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCAcc                                                                <  1:262215/248‑1 (MQ=255)
                                                                                                                                                                                      aGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGACGTCTCCCGTCATACGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATAc                                                        >  1:197167/1‑251 (MQ=255)
                                                                                                                                                                                          tGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCTTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCACGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGAGCAGGACACCCACGAACGCACCGTTGATCCCcac                                                    >  2:419354/1‑251 (MQ=255)
                                                                                                                                                                                          tGCGCACCTTACTGCGTAGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACCCCTACGAACGCACCGTTGATACCcac                                                    >  2:326222/1‑251 (MQ=255)
                                                                                                                                                                                                    aCTGCGTCGGGTGAAGAAATTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAAcct                                          >  2:167476/1‑249 (MQ=255)
                                                                                                                                                                                                       gCGTGGGGGGAAGAAGTTCTCTTCGCCGCCCCCCTTCCTCCGTTTTGGTCATTT‑GAATTTGAAGGACCCCGCGGCGCAGACCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACgttg                                        <  1:394652/250‑4 (MQ=255)
                                                                                                                                                                                                        cGTCGGGTGAAGAAGTTCTCGACGCAGTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGAT‑GGATTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTAGCCAGTCCGCGTCTGGTCCCGCCAGCAACTAATGGGTACCGGGGGGGAAGATGCGAAGGACCCCTAAGATAGCACAGTTGATACCCACATTAAAACGTTGCg                                      >  2:387618/1‑251 (MQ=255)
                                                                                                                                                                                                         gtggggtaaaaaagttccggaggcggTCTCCCGTCACCCGTGTGGGAAATTTTAAATG‑AAATGACCCCGGGGCGCAGA‑CAACTTGGTTGGACACGCCATTAGCGCTGACTCGTGT‑TGAGTTTTTATTGGGGAAGCCGTTCCTCAAGTCCCCGGGCCGCGTTGGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGt                                     <  2:73296/226‑1 (MQ=255)
                                                                                                                                                                                                               caccccccttttccgatctgtcccccGTCATCCGTTTTGGTAATTTGAAATT‑AAATAAACCCGGGGCGCATATCAGCT‑GTTTTGCCACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAg                               <  2:202321/228‑1 (MQ=255)
                                                                                                                                                                                                                               gCCGTCTCCCGTCCTCCGTATTGGTCATTTTGAATT‑GAATGAACCAGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GAATGAGTTTTTATTGCTGAAGACGTTACTAAAGTCGCCAGGCCGCGTATGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCCGGACACCTACGATCGACCCGTTGATCCCCACATTAAAACGTTGCGTGCCAAGCGGCGCGCCATCAAcc               >  2:448636/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gTCTCCCGTCATCCGTATTGGTCATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAACGCTGACGCG‑GTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAg                               >  1:202321/1‑232 (MQ=255)
                                                                                                                                                                                                                                            aTCCGTATTGGACATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAGCGCTGACTCG‑GTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCAACATTAAAACGCTGCGTGCCAAGCTGCGCGCCCTCAACCCCGCTCTTTCAcc  >  2:450372/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
CGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTG‑CTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATT‑GAATGAACCCGCGGCGCAGATCAGCT‑GGTTTGACACGCCATTAGCGCTGACTCG‑GTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACC  >  NC_000913/4636165‑4636650

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