Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,636,466 G→A A154T (GCG→ACG)  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,4660GA96.5% 167.1 / NA 59A154T (GCG→ACG) creBDNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (38/18);  minor base C (1/1);  total (39/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCA‑CCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCAT  >  NC_000913/4636219‑4636709
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
gcTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGCATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAacgc                                                                                                                                                                                                                                                   >  1:58268/1‑251 (MQ=255)
gcTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCCTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAacgc                                                                                                                                                                                                                                                   >  2:119306/1‑251 (MQ=255)
     aTCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGGGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGAcg                                                                                                                                                                                                                                              >  2:131973/1‑250 (MQ=255)
       ccGGCGTTACCTGTACTGTTACTGACGGCCCGGAGTGAAGAGGTCGCTCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGCGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAGCCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAAAGCTGACTCg                                                                                                                                                                                                                                            >  2:56928/1‑251 (MQ=255)
         ggCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCggg                                                                                                                                                                                                                                          >  2:350513/1‑250 (MQ=255)
          gCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTa                                                                                                                                                                                                                                         <  1:512392/251‑1 (MQ=255)
                 cggTACTGTTCTTGACGGCCCGAATTGAAGAGGTCGATCTCTTGCTTGGGCTGGAAATTGGTGCTGACGACTAGGTGGCTAAACCGTTTTCACCCCGCGAAGGGTGCTCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCTTCTCCCTTCATCCTTATTGGTCATTTTGAATTGAATGAACCCGTGGCGCAGATCAACGGG‑TGTGACACGCCATTAACGCTGACGCGGTATGAGttt                                                                                                                                                                                                                                  <  2:394870/249‑1 (MQ=255)
                     aCTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTAt                                                                                                                                                                                                                              <  1:323248/251‑1 (MQ=255)
                          tgcTGACGGCCCGAAGTGAATAGTTCGATCGCCTGCTTGGGCTGTAGATTGGTGCTGACGACGATGTGGCTAAACCGTTTTTCCCCCGCGAGGTGTGCGCCTGGGTGCGCCCTTTATTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCGGTCATCCCTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTAttgct                                                                                                                                                                                                                          <  2:439896/248‑3 (MQ=255)
                                    ccGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTACTGACGACTATGGGGCTAAACCGTTTTCACCCCGCGAAGTGAGCGCCAGGGTGCGCACCTTACTGCGTCGGGCGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGCATGGGTTTATATTGCTGAGGACGTTAc                                                                                                                                                                                                               >  2:492588/1‑251 (MQ=255)
                                           gAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGt                                                                                                                                                                                                         <  1:577571/250‑1 (MQ=255)
                                              gagGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTcgcc                                                                                                                                                                                                     >  1:12238/1‑248 (MQ=255)
                                                ggTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTcgccag                                                                                                                                                                                                   <  2:7311/251‑6 (MQ=255)
                                                 gTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTCAACCGTTTTCACACCGCGAAGTGTGCGCCAGGGGGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCCTCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCCTTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTcgcccgg                                                                                                                                                                                                  >  2:627592/1‑245 (MQ=255)
                                                   cgagcGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGcc                                                                                                                                                                                               <  1:350398/247‑1 (MQ=255)
                                                    gATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCcg                                                                                                                                                                                              >  2:421100/1‑251 (MQ=255)
                                                        cccTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTc                                                                                                                                                                                          <  1:511383/250‑1 (MQ=255)
                                                         ccTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTc                                                                                                                                                                                          >  2:511383/1‑250 (MQ=255)
                                                          cTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTg                                                                                                                                                                                        >  1:204090/1‑251 (MQ=255)
                                                              ttGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTc                                                                                                                                                                                     <  2:117364/250‑1 (MQ=255)
                                                               tGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCCGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTccc                                                                                                                                                                                   >  2:303407/1‑251 (MQ=255)
                                                                gggCTGGAAATTGGTGCTTACGACTATGTGGCTAACCCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATGGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTccc                                                                                                                                                                                   <  1:615555/250‑1 (MQ=255)
                                                                       aaaTTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCa                                                                                                                                                                            >  1:312264/1‑250 (MQ=255)
                                                                                       ctatgTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGtgt                                                                                                                                                           >  2:149188/5‑249 (MQ=255)
                                                                                           gTGGCTAAACCGTTTTCCCCCCGCTAAGTGGGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCGCGTCGCCGTCTCCCGACTTCCGTAGTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGAcgcg                                                                                                                                                                                                                                            <  1:589709/167‑4 (MQ=255)
                                                                                           gTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTCCTGCGTAGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCCTCCGTATTGGTCATTTTGACTTGGATGAACCCGCGGCGCAAATCAGCTGG‑TTTGACACGCCCTTAACGCTGAcgcg                                                                                                                                                                                                                                            >  2:589709/1‑164 (MQ=255)
                                                                                                  aaCCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGCCGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATgcg                                                                                                                                                >  1:329279/1‑251 (MQ=255)
                                                                                                         ttCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTCATGGATAGCGTGTGGGAAGATGCGCAGGac                                                                                                                                          >  1:328963/1‑250 (MQ=255)
                                                                                                              cccGCGAAGTGTGCGCCAGTGTGCCCCCTTTCCTGCGTCGGGTAAAGAAGTTCTCGTCGCCGTCTCCCTTCATCCGTATTGGTCATTTTGAATGGAATGACCCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTGCTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTAc                                                                                                                                    <  1:270189/251‑1 (MQ=255)
                                                                                                                    aagtggcccccaggtggcgcCCCTTACTGGGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGCCACCCGATTTGGTATTTTTAAATTGAATGACCCCGGGGCCCAGATCAGTTGG‑TTTGACACGCCATTAACGCTGTCGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGc                                                                                                                              <  1:460151/235‑1 (MQ=255)
                                                                                                                     aGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCCTTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACg                                                                                                                                   >  2:90094/1‑245 (MQ=255)
                                                                                                                          gcgcCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTCATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGtt                                                                                                                        >  2:283526/1‑250 (MQ=255)
                                                                                                                                    gcgcACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCcac                                                                                                               >  2:24130/1‑250 (MQ=255)
                                                                                                                                     cgcACCTTACTGCGTCGGGTGAAGACGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGACTGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGcc                                                                                                                                                                                >  2:604425/1‑184 (MQ=255)
                                                                                                                                              cTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACgtt                                                                                                    >  2:166902/1‑249 (MQ=255)
                                                                                                                                               tGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACgttg                                                                                                   >  1:95485/1‑248 (MQ=255)
                                                                                                                                                 cgtgggttgaagaatgtcgcgtcGCCGTCCCCCGTCATCCGTATTGGTTATTTTGATTTAATTGAACCCGCGGCGCAAATCAGTTGT‑TTTGCCACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGTTTCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGc                                                                                                  <  1:417476/228‑1 (MQ=255)
                                                                                                                                                  gTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATAAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGt                                                                                                >  2:391575/1‑251 (MQ=255)
                                                                                                                                                        tGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATgg                                                                                                                                                                    >  1:587863/1‑177 (MQ=255)
                                                                                                                                                            gaatttCTCGTCGCCGTCTCCCTTCATCCTTATTGGTCATTTTAATTTGAATGAACCCGCGGCGCAGATCAGCGGG‑TTTAACACGCCATTACCGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTgc                                                                                      <  1:397570/247‑1 (MQ=255)
                                                                                                                                                            gaagtgccctcccccgcctcCCTTCATCCTTATCGGTTATTTTGAATGGAATGACCCCGCTGCGCAGAACA‑CTGGGTTTCACCCGCCATTAACGCTGACGCGGTATGAGTTCTTATTGCTGAAGACGTTACCCAAGCCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATAACCACATTAAAACGTTGCGTGCCAAGCTgc                                                                                      <  1:338264/234‑1 (MQ=255)
                                                                                                                                                                        gCCGTCGCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCCCAGATCAGCTGG‑TTTGACACGCCATTACCGCTGCCGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCAACATTAAAACGTTGCGTGCCAAGCTGCGCGCAATCAAcc                                                                          >  1:257142/1‑251 (MQ=255)
                                                                                                                                                                            tctcCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCg                                                                       >  1:89430/1‑250 (MQ=255)
                                                                                                                                                                               cccGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGctct                                                                   >  1:249341/1‑251 (MQ=255)
                                                                                                                                                                               cccGTCATCCGTATGGGTCATTTTGAATTGAATGACCCCGCGGCGAAGATCAGTTGG‑TTTGACACGCCATTAACGCTGACGCGGTAGGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGGTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCt                                                                   <  1:480694/251‑1 (MQ=255)
                                                                                                                                                                                    cATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCAc                                                              >  1:92914/1‑251 (MQ=255)
                                                                                                                                                                                           aTTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGTCACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGCTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTa                                                        >  2:143138/1‑250 (MQ=255)
                                                                                                                                                                                             tggtcATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCGGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAAt                                                      >  1:431702/5‑250 (MQ=255)
                                                                                                                                                                                                        aatttaatgaccctgcgtgccAGATCACCGGG‑TTTGACACGCCATTAACGCTGACCCGGTTGGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCGGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCTCAGGACACCTACGATCGCTCCTTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCTCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGc                                          <  1:30458/231‑1 (MQ=255)
                                                                                                                                                                                                        aattgaatGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAg                                                                                                                                                     >  2:453737/1‑144 (MQ=255)
                                                                                                                                                                                                            gaatgaaCCCGCGGCGAAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCAGATCTTTCACCGATTAATACTCATCGCGGCATgg                                      >  1:79234/1‑251 (MQ=255)
                                                                                                                                                                                                                    ccGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGc                              >  1:130873/1‑251 (MQ=255)
                                                                                                                                                                                                                       cggcgCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCCTGTGATATAGCCTg                           >  2:346002/1‑251 (MQ=255)
                                                                                                                                                                                                                            tagacacggtGG‑TTTGACCAGCCATTAGCGCTGCGGCGGTATGAGTTTGTATTGCGGAAGAGGTTATTCAA‑TCGGCCAGGCCGCTATTGGGACCGCCATCAACTAAGGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGAGCGCACCGTTGATACCCGCATTAAAACGGTGTGTGCCAAGCTGCGTGCCATCAGCCCCGAGCTTTCGCCGATTAATACTCATCGCGGCATGGGATATAGCCTGAggg                       <  1:270683/241‑1 (MQ=255)
                                                                                                                                                                                                                                 cAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGt                 >  2:154153/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAAt                                                      >  2:43706/1‑212 (MQ=255)
                                                                                                                                                                                                                                   gCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTaa               <  1:119306/251‑1 (MQ=255)
                                                                                                                                                                                                                                     tGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGAGGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGAGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATAGCGGCATGGGATATAGCCTGAGGGGCCTGTAATg             >  1:304237/1‑251 (MQ=255)
                                                                                                                                                                                                                                       g‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTgcgcg                                                                                   >  1:119331/1‑179 (MQ=255)
                                                                                                                                                                                                                                               acGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCCGCAACTCATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTACTGCGTATCGGc    >  1:560012/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 gCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCG‑CCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCAt  >  2:502871/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
GCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCA‑CCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCAT  >  NC_000913/4636219‑4636709

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