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

GCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCAT  >  NC_000913/4636219‑4636709
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
gcTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGCATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAacgc                                                                                                                                                                                                                                                   >  1:58268/1‑251 (MQ=255)
gcTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCCTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAacgc                                                                                                                                                                                                                                                   >  2:119306/1‑251 (MQ=255)
     aTCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGGGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGAcg                                                                                                                                                                                                                                              >  2:131973/1‑250 (MQ=255)
       ccGGCGTTACCTGTACTGTTACTGACGGCCCGGAGTGAAGAGGTCGCTCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGCGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAGCCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAAAGCTGACTCg                                                                                                                                                                                                                                            >  2:56928/1‑251 (MQ=255)
         ggCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCggg                                                                                                                                                                                                                                          >  2:350513/1‑250 (MQ=255)
          gCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTa                                                                                                                                                                                                                                         <  1:512392/251‑1 (MQ=255)
                 cggTACTGTTCTTGACGGCCCGAATTGAAGAGGTCGATCTCTTGCTTGGGCTGGAAATTGGTGCTGACGACTAGGTGGCTAAACCGTTTTCACCCCGCGAAGGGTGCTCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCTTCTCCCTTCATCCTTATTGGTCATTTTGAATTGAATGAACCCGTGGCGCAGATCAACGGGTGTGACACGCCATTAACGCTGACGCGGTATGAGttt                                                                                                                                                                                                                                  <  2:394870/249‑1 (MQ=255)
                     aCTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTAt                                                                                                                                                                                                                              <  1:323248/251‑1 (MQ=255)
                          tgcTGACGGCCCGAAGTGAATAGTTCGATCGCCTGCTTGGGCTGTAGATTGGTGCTGACGACGATGTGGCTAAACCGTTTTTCCCCCGCGAGGTGTGCGCCTGGGTGCGCCCTTTATTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCGGTCATCCCTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTAttgct                                                                                                                                                                                                                          <  2:439896/248‑3 (MQ=255)
                                    ccGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTACTGACGACTATGGGGCTAAACCGTTTTCACCCCGCGAAGTGAGCGCCAGGGTGCGCACCTTACTGCGTCGGGCGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGCATGGGTTTATATTGCTGAGGACGTTAc                                                                                                                                                                                                               >  2:492588/1‑251 (MQ=255)
                                           gAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGt                                                                                                                                                                                                         <  1:577571/250‑1 (MQ=255)
                                              gagGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTcgcc                                                                                                                                                                                                     >  1:12238/1‑248 (MQ=255)
                                                ggTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTcgccag                                                                                                                                                                                                   <  2:7311/251‑6 (MQ=255)
                                                 gTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTCAACCGTTTTCACACCGCGAAGTGTGCGCCAGGGGGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCCTCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCCTTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTcgcccgg                                                                                                                                                                                                  >  2:627592/1‑245 (MQ=255)
                                                   cgagcGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGcc                                                                                                                                                                                               <  1:350398/247‑1 (MQ=255)
                                                    gATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCcg                                                                                                                                                                                              >  2:421100/1‑251 (MQ=255)
                                                        cccTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTc                                                                                                                                                                                          <  1:511383/250‑1 (MQ=255)
                                                         ccTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTc                                                                                                                                                                                          >  2:511383/1‑250 (MQ=255)
                                                          cTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTg                                                                                                                                                                                        >  1:204090/1‑251 (MQ=255)
                                                              ttGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTc                                                                                                                                                                                     <  2:117364/250‑1 (MQ=255)
                                                               tGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCCGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTccc                                                                                                                                                                                   >  2:303407/1‑251 (MQ=255)
                                                                gggCTGGAAATTGGTGCTTACGACTATGTGGCTAACCCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATGGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTccc                                                                                                                                                                                   <  1:615555/250‑1 (MQ=255)
                                                                       aaaTTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCa                                                                                                                                                                            >  1:312264/1‑250 (MQ=255)
                                                                                       ctatgTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGtgt                                                                                                                                                           >  2:149188/5‑249 (MQ=255)
                                                                                           gTGGCTAAACCGTTTTCCCCCCGCTAAGTGGGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCGCGTCGCCGTCTCCCGACTTCCGTAGTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGAcgcg                                                                                                                                                                                                                                            <  1:589709/167‑4 (MQ=255)
                                                                                           gTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTCCTGCGTAGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCCTCCGTATTGGTCATTTTGACTTGGATGAACCCGCGGCGCAAATCAGCTGGTTTGACACGCCCTTAACGCTGAcgcg                                                                                                                                                                                                                                            >  2:589709/1‑164 (MQ=255)
                                                                                                  aaCCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGCCGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATgcg                                                                                                                                                >  1:329279/1‑251 (MQ=255)
                                                                                                         ttCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTCATGGATAGCGTGTGGGAAGATGCGCAGGac                                                                                                                                          >  1:328963/1‑250 (MQ=255)
                                                                                                              cccGCGAAGTGTGCGCCAGTGTGCCCCCTTTCCTGCGTCGGGTAAAGAAGTTCTCGTCGCCGTCTCCCTTCATCCGTATTGGTCATTTTGAATGGAATGACCCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTGCTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTAc                                                                                                                                    <  1:270189/251‑1 (MQ=255)
                                                                                                                    aagtggcccccaggtggcgcCCCTTACTGGGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGCCACCCGATTTGGTATTTTTAAATTGAATGACCCCGGGGCCCAGATCAGTTGGTTTGACACGCCATTAACGCTGTCGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGc                                                                                                                              <  1:460151/235‑1 (MQ=255)
                                                                                                                     aGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCCTTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACg                                                                                                                                   >  2:90094/1‑245 (MQ=255)
                                                                                                                          gcgcCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTCATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGtt                                                                                                                        >  2:283526/1‑250 (MQ=255)
                                                                                                                                    gcgcACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCcac                                                                                                               >  2:24130/1‑250 (MQ=255)
                                                                                                                                     cgcACCTTACTGCGTCGGGTGAAGACGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGACTGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGcc                                                                                                                                                                                >  2:604425/1‑184 (MQ=255)
                                                                                                                                              cTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACgtt                                                                                                    >  2:166902/1‑249 (MQ=255)
                                                                                                                                               tGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACgttg                                                                                                   >  1:95485/1‑248 (MQ=255)
                                                                                                                                                 cgtgggttgaagaatgtcgcgtcGCCGTCCCCCGTCATCCGTATTGGTTATTTTGATTTAATTGAACCCGCGGCGCAAATCAGTTGTTTTGCCACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGTTTCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGc                                                                                                  <  1:417476/228‑1 (MQ=255)
                                                                                                                                                  gTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATAAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGt                                                                                                >  2:391575/1‑251 (MQ=255)
                                                                                                                                                        tGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATgg                                                                                                                                                                    >  1:587863/1‑177 (MQ=255)
                                                                                                                                                            gaatttCTCGTCGCCGTCTCCCTTCATCCTTATTGGTCATTTTAATTTGAATGAACCCGCGGCGCAGATCAGCGGGTTTAACACGCCATTACCGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTgc                                                                                      <  1:397570/247‑1 (MQ=255)
                                                                                                                                                            gaagtgccctcccccgcctcCCTTCATCCTTATCGGTTATTTTGAATGGAATGACCCCGCTGCGCAGAACACTGGGTTTCACCCGCCATTAACGCTGACGCGGTATGAGTTCTTATTGCTGAAGACGTTACCCAAGCCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATAACCACATTAAAACGTTGCGTGCCAAGCTgc                                                                                      <  1:338264/234‑1 (MQ=255)
                                                                                                                                                                        gCCGTCGCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCCCAGATCAGCTGGTTTGACACGCCATTACCGCTGCCGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCAACATTAAAACGTTGCGTGCCAAGCTGCGCGCAATCAAcc                                                                          >  1:257142/1‑251 (MQ=255)
                                                                                                                                                                            tctcCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCg                                                                       >  1:89430/1‑250 (MQ=255)
                                                                                                                                                                               cccGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGctct                                                                   >  1:249341/1‑251 (MQ=255)
                                                                                                                                                                               cccGTCATCCGTATGGGTCATTTTGAATTGAATGACCCCGCGGCGAAGATCAGTTGGTTTGACACGCCATTAACGCTGACGCGGTAGGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGGTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCt                                                                   <  1:480694/251‑1 (MQ=255)
                                                                                                                                                                                    cATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCAc                                                              >  1:92914/1‑251 (MQ=255)
                                                                                                                                                                                           aTTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGTCACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGCTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTa                                                        >  2:143138/1‑250 (MQ=255)
                                                                                                                                                                                             tggtcATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCGGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAAt                                                      >  1:431702/5‑250 (MQ=255)
                                                                                                                                                                                                        aatttaatgaccctgcgtgccAGATCACCGGGTTTGACACGCCATTAACGCTGACCCGGTTGGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCGGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCTCAGGACACCTACGATCGCTCCTTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCTCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGc                                          <  1:30458/231‑1 (MQ=255)
                                                                                                                                                                                                        aattgaatGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAg                                                                                                                                                     >  2:453737/1‑144 (MQ=255)
                                                                                                                                                                                                            gaatgaaCCCGCGGCGAAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCAGATCTTTCACCGATTAATACTCATCGCGGCATgg                                      >  1:79234/1‑251 (MQ=255)
                                                                                                                                                                                                                    ccGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGc                              >  1:130873/1‑251 (MQ=255)
                                                                                                                                                                                                                       cggcgCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCCTGTGATATAGCCTg                           >  2:346002/1‑251 (MQ=255)
                                                                                                                                                                                                                            tagacacggtGGTTTGACCAGCCATTAGCGCTGCGGCGGTATGAGTTTGTATTGCGGAAGAGGTTATTCAATCGGCCAGGCCGCTATTGGGACCGCCATCAACTAAGGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGAGCGCACCGTTGATACCCGCATTAAAACGGTGTGTGCCAAGCTGCGTGCCATCAGCCCCGAGCTTTCGCCGATTAATACTCATCGCGGCATGGGATATAGCCTGAggg                       <  1:270683/241‑1 (MQ=255)
                                                                                                                                                                                                                                 cAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGt                 >  2:154153/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAAt                                                      >  2:43706/1‑212 (MQ=255)
                                                                                                                                                                                                                                   gCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTaa               <  1:119306/251‑1 (MQ=255)
                                                                                                                                                                                                                                     tGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGAGGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGAGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATAGCGGCATGGGATATAGCCTGAGGGGCCTGTAATg             >  1:304237/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTgcgcg                                                                                   >  1:119331/1‑179 (MQ=255)
                                                                                                                                                                                                                                               acGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCCGCAACTCATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTACTGCGTATCGGc    >  1:560012/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 gCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCAt  >  2:502871/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
GCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCAT  >  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: