Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,636,175 C→T L57L (CTG→TTG)  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,1750CT98.4% 197.9 / ‑4.5 65L57L (CTG→TTG) creBDNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (32/31);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

AAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGG‑TACCCGACGTCATGATT‑CTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATT  >  NC_000913/4635930‑4636408
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         
aaGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCTAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGcct                                                                                                                                                                                                                                      <  2:136470/251‑2 (MQ=255)
aaGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGCTGTTGGTTTGcct                                                                                                                                                                                                                                      >  2:512392/1‑250 (MQ=255)
 aGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGACGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGATCTTTGAGCGAGGCTTGCCGGGGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCAGGATT‑CTCGATGTTGGTTTGcctg                                                                                                                                                                                                                                     >  2:480694/1‑249 (MQ=255)
         aaGGTTCGCCGAAAAACGCGATTAGCGCGGTTCATGTCATGCCGTGGCGGCAATAACAGAGGCGATTTAGGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGAGT‑CTCGATGTTGGTGTGccttatatcagc                                                                                                                                                                                                                             <  2:123948/251‑11 (MQ=255)
           ggTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCgg                                                                                                                                                                                                                           >  1:376213/1‑251 (MQ=255)
            gTTCGCCGAAAAACGCGATTAGAGCGGTGTCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTAGGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGTTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGc                                                                                                                                                                                                                          <  1:381401/251‑1 (MQ=255)
                  cGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGaa                                                                                                                                                                                                                    <  1:313242/251‑1 (MQ=255)
                   gAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAAt                                                                                                                                                                                                                   >  1:7311/1‑251 (MQ=255)
                       aactcgATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAAttatg                                                                                                                                                                                                               <  1:363465/247‑4 (MQ=255)
                       aaCGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAAttatg                                                                                                                                                                                                               >  1:347406/1‑248 (MQ=255)
                        aCGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGCGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATgc                                                                                                                                                                                                              >  2:190960/1‑251 (MQ=255)
                         cgctatcatcgcgGTTACTGTCATGCCGTGGCGGCAAGAACAGAGGCGATTTATGCAACGGGAAACGTTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTGGCAGCAGGAAGGTTTTGCAGGCGATGTCGTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATG‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATgcc                                                                                                                                                                                                             <  2:151184/242‑1 (MQ=255)
                         cgcgATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTGGAATTATgcc                                                                                                                                                                                                             <  1:190960/251‑1 (MQ=255)
                          gcgATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGGTTGGTTAGTGGAAGAGGAGAAAGGGATAGCCGACACGCTGATCTAAATGTTGCAGAAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTTCCGGAGCTGGATAAAGCTCGCCAGCAGGCTACC‑GACGTCA‑GATTTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGagtaatgccg                                                                                                                                                                                                            >  1:534682/1‑242 (MQ=255)
                             aTTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTTCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCa                                                                                                                                                                                                         <  1:354792/251‑1 (MQ=255)
                                 gcgcGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCCGCAGG‑CTCCCGACGTCATGATT‑CTCGAAGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGtt                                                                                                                                                                                                      >  2:197830/1‑250 (MQ=255)
                                  cgcgGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTAc                                                                                                                                                                                                    <  1:592637/251‑1 (MQ=255)
                                          cTGTCATGACGTGGCGGCAATAACAGAGACGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCCAGGGAAAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGAAGTCATGATT‑CTCGATGTTGGCTTGCCTGATATCAGAGGCTTTGAATTATGCCGCCCGTTACTGGCGCTc                                                                                                                                                                                            >  2:116370/1‑251 (MQ=255)
                                            gTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTccc                                                                                                                                                                                          >  2:392681/1‑250 (MQ=255)
                                              cATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTTATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATc                                                                                                                                                                                        <  1:254670/251‑1 (MQ=255)
                                                           cAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTTGCGCTCCATCCGGCGTTACCTGt                                                                                                                                                                           <  1:119765/251‑1 (MQ=255)
                                                           cAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGt                                                                                                                                                                           >  1:390622/1‑251 (MQ=255)
                                                            aaTAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGATTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGt                                                                                                                                                                           >  1:267930/1‑250 (MQ=255)
                                                               aaCAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCAGCCGGCGTTACCTGTACTg                                                                                                                                                                       <  1:508113/251‑1 (MQ=255)
                                                                    agGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCt                                                                                                                                                                  <  1:198970/251‑1 (MQ=255)
                                                                         aTTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACgg                                                                                                                                                             <  1:294104/251‑1 (MQ=255)
                                                                                        aaCGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAgg                                                                                                                                              >  2:589089/1‑251 (MQ=255)
                                                                                                ggTTAGTGGAAGATGAGGAAGGGATAGCCGACACGCTGGTCTTCATGTTGCAGCAGGAAGGTTGTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGTTTTGCCTGATATCAGCGTCTTTGAATTATGCCGCCAGGTACTGGCGCTCCATCCTGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGc                                                                                                                                      <  1:627172/251‑1 (MQ=255)
                                                                                                ggTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGc                                                                                                                                      >  1:122534/1‑251 (MQ=255)
                                                                                                        gAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTgg                                                                                                                              <  1:310030/251‑1 (MQ=255)
                                                                                                         aaGATGAGCAAGGGATAGCTGACACGCAGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCTCGACGGCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTGCCGGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTggg                                                                                                                             <  1:599189/251‑1 (MQ=255)
                                                                                                         aaGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTggg                                                                                                                             <  2:244335/251‑1 (MQ=255)
                                                                                                           gATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCt                                                                                                                           >  1:381017/1‑251 (MQ=255)
                                                                                                                gCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACATGTACTGTTCCTGACGGCCCGAAGGGAAGAGGTCGATCGCCTGCTTGGGCTggga                                                                                                                       >  2:417476/1‑248 (MQ=255)
                                                                                                                 cAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAAt                                                                                                                     >  2:323248/1‑251 (MQ=255)
                                                                                                                  aaGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAAtt                                                                                                                    >  2:363218/1‑251 (MQ=255)
                                                                                                                         aGCCGACACGCTTGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCt                                                                                                              <  2:337972/250‑1 (MQ=255)
                                                                                                                            cGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGGGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGACTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTgacg                                                                                                          >  2:397570/1‑251 (MQ=255)
                                                                                                                               cacGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACt                                                                                                       <  1:132115/251‑1 (MQ=255)
                                                                                                                                     gttCTACATGTTGCAGCAGGAAGGTTTGGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CGCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGTGCTGGAAATTTGTGCTGACGACTATGTgg                                                                                                 <  2:566822/249‑1 (MQ=255)
                                                                                                                                            aTGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGCCGACTATGTGGCTAAAcc                                                                                          >  2:135494/1‑251 (MQ=255)
                                                                                                                                              gTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGt                                                                                        <  1:293437/251‑1 (MQ=255)
                                                                                                                                                                  ccGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGAAGG‑CTCCCGCCGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGAACTGTTCCTGACGGCCCGAAGGGAAGAGGGCGATCGCCTGCTTGGGCTGGAAATTTGTGCTGACGACTATGTGGCTAAACGGCTTTCACCCCGCGCAGTGTgc                                                                    >  2:338264/1‑251 (MQ=255)
                                                                                                                                                                    gTCGATGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGGGAAGTGTgcgc                                                                  <  2:528296/251‑1 (MQ=255)
                                                                                                                                                                       gAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGcc                                                                 >  1:339769/1‑249 (MQ=255)
                                                                                                                                                                              ttGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTgcg                                                         <  2:441699/250‑1 (MQ=255)
                                                                                                                                                                                gAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCAc                                                      <  2:122534/251‑1 (MQ=255)
                                                                                                                                                                                   cGAGGCTTGCCGGTGCAGGATAAAGCTCGCCAGCAGG‑CGCCCGACGTCATGATT‑ATCGATGTTTGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCAACtt                                                   <  2:547477/251‑1 (MQ=255)
                                                                                                                                                                                    gAGGCTTGCCGGTGCTGGCTAACGCTCGCCAGCAGG‑CTCCCGCCGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCCCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGCAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGCGCGCACCtt                                                   >  2:117809/1‑250 (MQ=255)
                                                                                                                                                                                        cTTGCCGGTGCTGGATAAAGCTCGCAAGCAGG‑TACCCGACGTCATGATT‑CTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGTTGCTGACGACTACGTGGCTAAACCGGTTTCAACCAGCGAAGTGTGGGCAAGGGAGAGCACCTTACTGc                                              >  1:385512/1‑251 (MQ=255)
                                                                                                                                                                                         ttGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCCGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCATGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTACACCGTTTTCACCCCGCGAAGTGTGCGCCAg                                                               >  2:227076/1‑233 (MQ=255)
                                                                                                                                                                                         ttGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAg                                                               <  1:227076/233‑1 (MQ=255)
                                                                                                                                                                                             cGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTACACCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGCCTATGTGGCTAAACCGTTTTCACCCAGCGAAGTCGGCGCCAGGGGGCGCACCTTACTGCGTCgg                                         >  1:537035/1‑251 (MQ=255)
                                                                                                                                                                                                  cTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTga                                     >  1:37283/1‑250 (MQ=255)
                                                                                                                                                                                                  cTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTAATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTga                                     >  1:508178/1‑250 (MQ=255)
                                                                                                                                                                                                                   agcaggctcccgccgtactGATT‑CTAGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCg                    >  2:615555/19‑251 (MQ=255)
                                                                                                                                                                                                                   agcaggctcccGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCg                    <  1:469524/243‑1 (MQ=255)
                                                                                                                                                                                                                     caggctcccGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGtc                  <  1:392681/245‑1 (MQ=255)
                                                                                                                                                                                                                             cGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTccc              >  2:334560/1‑247 (MQ=255)
                                                                                                                                                                                                                             cGACGTCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTccc              <  1:334560/247‑1 (MQ=255)
                                                                                                                                                                                                                                  tCATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCgg     >  2:150659/1‑250 (MQ=255)
                                                                                                                                                                                                                                   cATGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGCAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCCGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGt     >  2:577571/1‑250 (MQ=255)
                                                                                                                                                                                                                                     tGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTc                                           <  2:64471/210‑1 (MQ=255)
                                                                                                                                                                                                                                     tGATT‑CTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTc                                           >  1:64471/1‑210 (MQ=255)
                                                                                                                                                                                                                                      gATT‑CTCGATGTTGGTTTGCCTGATAGCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCTCCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCAGCCGTAtt  <  2:37283/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         
AAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGG‑TACCCGACGTCATGATT‑CTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATT  >  NC_000913/4635930‑4636408

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