Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,636,146 2 bp→CT coding (140‑141/690 nt) 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,1460TC90.3% 160.2 / 4.6 64V47A (GTA→GCT) creBDNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base T (4/2);  new base C (31/25);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01
*NC_0009134,636,1470AT90.5% 173.0 / 8.3 64V47A (GTA→GCT) creBDNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base T (33/24);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

AACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATT‑GTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTT  >  NC_000913/4635899‑4636380
                                                                                                                                                                                                                                                       ||                                                                                                                                                                                                                                          
aacgcGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGGGCTGGATAAAGCTCGCAAGCAGGTAcc                                                                                                                                                                                                                                          >  2:196431/1‑251 (MQ=255)
      ctcgctTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGt                                                                                                                                                                                                                                    >  1:360979/1‑251 (MQ=255)
                     tACTCCGACAAAGTTGAAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGAAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGAGGTCAGGATTCTCGATgg                                                                                                                                                                                                                     >  2:103301/1‑250 (MQ=255)
                            aCAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGtttg                                                                                                                                                                                                              >  1:447001/1‑248 (MQ=255)
                            aCAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGATCCTGCCATGCCGTGGCGGCAATAACAGAGGCGATTGATGCAACGGGAAACGGTCTGGTTATTGGAAGATGAGCAAGTGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGATGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGtttg                                                                                                                                                                                                              <  2:94436/251‑4 (MQ=255)
                               aaGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGGCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCGGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGcct                                                                                                                                                                                                           <  1:32489/251‑2 (MQ=255)
                                                    aaaaCGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAAtt                                                                                                                                                                                       >  1:291267/1‑250 (MQ=255)
                                                          ggATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTATATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATgccgc                                                                                                                                                                               <  2:79231/250‑1 (MQ=255)
                                                             ttAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATgccgcccg                                                                                                                                                                            >  2:480165/1‑249 (MQ=255)
                                                               aGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCg                                                                                                                                                                                                 <  2:428645/229‑1 (MQ=255)
                                                               aGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCg                                                                                                                                                                                                 >  1:428645/1‑229 (MQ=255)
                                                                        ccTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTAGGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCt                                                                                                                                                                 <  2:360979/251‑1 (MQ=255)
                                                                         cTGTCCTGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTCGAATT‑ATGCCGCCAGTTACTGGCGCTc                                                                                                                                                                >  2:346616/1‑251 (MQ=255)
                                                                         cTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGGTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATTCCGCCAGTTAATGGCGCTc                                                                                                                                                                >  2:263977/1‑251 (MQ=255)
                                                                          tGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTcc                                                                                                                                                               >  1:306133/1‑251 (MQ=255)
                                                                            tCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGACTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCAt                                                                                                                                                             >  2:403827/1‑251 (MQ=255)
                                                                             cATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGTCTTTGAATT‑GTGCCGCCAGTTACTGGCGCTCCATc                                                                                                                                                            >  1:74350/1‑251 (MQ=255)
                                                                                gCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCgg                                                                                                                                                         <  2:470257/251‑1 (MQ=255)
                                                                                     ggcggcAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTa                                                                                                                                                    <  2:423365/251‑1 (MQ=255)
                                                                                          cAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGt                                                                                                                                               >  1:347478/1‑251 (MQ=255)
                                                                                           aaTAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGATCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGt                                                                                                                                               >  1:255824/1‑250 (MQ=255)
                                                                                                cAGAGGCGATTTATGCAACGGGAAACGGTCTGGTGAGTGGCAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTTCAGAAGGAAGGTTTTGACGGCGCCGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGATCCCGACGACATGATTCTCGCTGTTGGTTTGCCTGCAATCAGCGGCTTTGCATTAGT‑AAGCCAGTTATTGGCGCTCCATCCgtggttaacgtaacggtt                                                                                                                                         >  2:409656/1‑234 (MQ=255)
                                                                                                    ggCGATATATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTg                                                                                                                                     <  1:29865/251‑1 (MQ=255)
                                                                                                              gCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCCTCCGGCGTTACCTGTACTGTTCCTGACGGCCCggg                                                                                                                           >  2:115016/1‑249 (MQ=255)
                                                                                                              gCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGaa                                                                                                                           <  1:263977/251‑1 (MQ=255)
                                                                                                               cAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGACTACATGTTGAAGAAGGAAGGTTTTGCCGTCGAGGACTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCCGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTCGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGGTACCTGTACTGTTCCTGACGGCCCGa                                                                                                                            >  2:464420/1‑249 (MQ=255)
                                                                                                                       aaCGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATT‑GTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAgg                                                                                                                  >  1:451005/1‑251 (MQ=255)
                                                                                                                           gTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCg                                                                                                               >  1:394211/1‑250 (MQ=255)
                                                                                                                                   aGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATTGTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTg                                                                                                       <  2:95761/250‑1 (MQ=255)
                                                                                                                                    gTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCt                                                                                                     <  1:243195/251‑1 (MQ=255)
                                                                                                                                            tGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGACTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTTGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGGCGGCCCGAAGTGAAGAGGTCGACCGCCTTCTTGGGCTgg                                                                                             >  2:15205/1‑251 (MQ=255)
                                                                                                                                                  aGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCATGCTCCCGACGTCATGATTCTCGATGTTGGTTTGTCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCa                                                                                                                                                              <  1:236647/180‑1 (MQ=255)
                                                                                                                                                        aGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATgcc                                                                                                                                                                                 >  2:418140/1‑155 (MQ=255)
                                                                                                                                                        aGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATgcc                                                                                                                                                                                 <  1:418140/155‑1 (MQ=255)
                                                                                                                                                          ccGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTgac                                                                               <  1:403827/251‑1 (MQ=255)
                                                                                                                                                           cGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTgacg                                                                              >  1:377718/1‑251 (MQ=255)
                                                                                                                                                                   tGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTAtgtg                                                                      >  1:264949/1‑251 (MQ=255)
                                                                                                                                                                        tACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTa                                                                  >  1:284138/1‑250 (MQ=255)
                                                                                                                                                                                     agGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCGAAACCGTTTTCAccc                                                    >  2:394652/1‑251 (MQ=255)
                                                                                                                                                                                     agGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATT‑GTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCAccc                                                    >  2:16259/1‑251 (MQ=255)
                                                                                                                                                                                      gtaaGGTTTTGCCGTCGATGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCAcccc                                                   <  1:127509/249‑1 (MQ=255)
                                                                                                                                                                                       gAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAACGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGACGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGCGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTCAACCGTTTTCAccccc                                                  >  2:476576/1‑250 (MQ=255)
                                                                                                                                                                                               tGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAgtgt                                          >  2:252816/1‑251 (MQ=255)
                                                                                                                                                                                                  cGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTgcg                                       <  1:480165/251‑1 (MQ=255)
                                                                                                                                                                                                  cGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGGCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCGCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTgcg                                       <  2:53615/251‑1 (MQ=255)
                                                                                                                                                                                                   gTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTgcgc                                      >  1:326079/1‑251 (MQ=255)
                                                                                                                                                                                                      gAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGcc                                     >  1:342039/1‑249 (MQ=255)
                                                                                                                                                                                                         gTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGGGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGGGCGCCAgggg                                >  2:265946/1‑250 (MQ=255)
                                                                                                                                                                                                          tCTTTGAGCGAGGCTTGCCGGTGCTGGATGAAGCTAGCCAGCAGGGTCCGGACGTCATGATTGGCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGTTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTCAAGAGGTCGCTCGCCTGCTTGGGCTGGAAATTGGTGCTGTCGACTATATGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAggggg                               >  1:203169/1‑249 (MQ=255)
                                                                                                                                                                                                             ttGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTgcg                             <  1:250559/250‑1 (MQ=255)
                                                                                                                                                                                                                 gCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCTATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCt                        <  2:209685/251‑1 (MQ=255)
                                                                                                                                                                                                                 gCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCCCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTATGAATT‑ATGCCGCCAGTTTCTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGTCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCt                        <  1:488818/251‑1 (MQ=255)
                                                                                                                                                                                                                  cGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGTCGTCATGATTCTCGAGGTTGGTCTGCCGGATATTAGCGGCTTTGAATT‑GTGCCGCCAGGTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCtt                       <  2:451005/251‑1 (MQ=255)
                                                                                                                                                                                                                    aGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTAc                     >  1:362677/1‑251 (MQ=255)
                                                                                                                                                                                                                      gCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTg                   <  1:252816/251‑1 (MQ=255)
                                                                                                                                                                                                                       cgtGCCGGTGCTGGATAAAGCCCGCCAGCAGGCTCCCGACGTCAGGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTGTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGACCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTACACCCCGCGAAGTGTGCGCCAGTGTGCGCACCTTACTGc                  <  2:255824/249‑1 (MQ=255)
                                                                                                                                                                                                                       cTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGc                  >  1:479954/1‑251 (MQ=255)
                                                                                                                                                                                                                          gCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGCCTGCCGGATATTAGCGGCTTTGAATT‑GTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTc               <  1:196431/251‑1 (MQ=255)
                                                                                                                                                                                                                             ggTGCTGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAGCTTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGCAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCggg            >  2:150218/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTgac                                                                               <  2:127312/179‑1 (MQ=255)
                                                                                                                                                                                                                                  tGGATAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTgac                                                                               >  1:127312/1‑179 (MQ=255)
                                                                                                                                                                                                                                     aTAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGc                  >  2:109297/1‑237 (MQ=255)
                                                                                                                                                                                                                                     aTAAAGCTCGCCAGCAGGCTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGc                  <  1:109297/237‑1 (MQ=255)
                                                                                                                                                                                                                                       aaaGCTCGCCAGCAGGCTCCCGTCGTCAGGATTCGCGATGTTGGTTTGCCTGATACCAGCGGCTTTGAATT‑ATGCCGCCAGTTACTGGCGCTCCAGCCTGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGTCTGCTTGGGCTGTAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCggttcaataagtt  <  2:394031/251‑12 (MQ=255)
                                                                                                                                                                                                                                                       ||                                                                                                                                                                                                                                          
AACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATT‑GTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTT  >  NC_000913/4635899‑4636380

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