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,4660GA94.4% 143.5 / 0.3 57A154T (GCG→ACG) creBDNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (34/17);  total (36/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.73e-01

GCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCAT  >  NC_000913/4636219‑4636709
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
gcTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAacgc                                                                                                                                                                                                                                                  <  1:445279/251‑1 (MQ=255)
       ccggggttccctgtacttgcccTGACGGCCCGAAGTGAAGAGGTCGATCGCCGGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGACCCCGCGGCGCAGATCAGC‑GGGTTTGACACGCCATTAACGCTGAcgcg                                                                                                                                                                                                                                           <  1:149506/231‑4 (MQ=255)
           cGTTACCTGTACTGTTCCTGACGGCCAGAAGTGAAGAGGTCGATCGCCTGATTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAACCCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGCAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGACCAGC‑TGGATTGACACGCCTTTAACGCTGACGCGGTAt                                                                                                                                                                                                                                       >  2:71973/1‑251 (MQ=255)
            gTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATg                                                                                                                                                                                                                                      >  1:482054/1‑251 (MQ=255)
             ttACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATgg                                                                                                                                                                                                                                     >  1:321587/1‑250 (MQ=255)
                ccTGTGCTGTTCCTGACGGCCCGGAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGGGCTGTCGACTAGGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCTTCTCCCGTCATCCGTTTTGGTCATTTTGAATTGAAGGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGtt                                                                                                                                                                                                                                  <  2:583153/251‑1 (MQ=255)
                       tGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTAttg                                                                                                                                                                                                                           >  1:397432/1‑251 (MQ=255)
                             tGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGCAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCAATTAACGCTGCCGCGGTATGAGTTTTTATTGCTGa                                                                                                                                                                                                                       >  2:311061/1‑249 (MQ=255)
                                      gAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACCATGTGGCTAAACCGTTTTACCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTc                                                                                                                                                                                                            <  1:334850/251‑1 (MQ=255)
                                        ggTGAAGAGGTCGATCGCCTGCTTGGGCGGGAAATTGGTGCTGACGACGATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCTCTTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCaa                                                                                                                                                                                                          <  2:36546/250‑1 (MQ=255)
                                                   cGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGCGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGcc                                                                                                                                                                                               >  1:327344/1‑251 (MQ=255)
                                                    gATCGCCTGCTTGGGCTGGAAATTGGTGCTGAAGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAAACGTTCTCGACGCCATCTCCCGTCATCCGTATTGGACATTTTGCATTGACTGAACCCGCGGCGCAGATCACC‑TGGTTTGCCCCACCATTAGCGATGACGCGGGATGAGTTTTTATTGTTGAAGTCGTTACTCAAGTCCCCGGGCcg                                                                                                                                                                                              >  2:482146/1‑251 (MQ=255)
                                                       cGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTTGGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCcgcg                                                                                                                                                                                            <  1:344495/250‑1 (MQ=255)
                                                       cGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCGGAAGACGTTTCTCCAGTCGCCAGGCCGCGt                                                                                                                                                                                           >  1:596986/1‑251 (MQ=255)
                                                          cTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCCTTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTg                                                                                                                                                                                        >  2:624645/1‑251 (MQ=255)
                                                            gCTTGGGCTGGAAATGGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAATTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTgg                                                                                                                                                                                       <  1:457204/250‑1 (MQ=255)
                                                                gggCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTccc                                                                                                                                                                                   <  1:311061/250‑1 (MQ=255)
                                                                         aTTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCCACt                                                                                                                                                                         >  2:100626/1‑251 (MQ=255)
                                                                             gTGCTGACGACTATGTGGCTAAACCGTTTTAACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTACGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATg                                                                                                                                                                     >  1:214944/1‑251 (MQ=255)
                                                                                         atgTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTgg                                                                                                                                                         >  2:466021/3‑251 (MQ=255)
                                                                                               cTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGa                                                                                                                                                       >  1:379891/1‑247 (MQ=255)
                                                                                                 aaaCCGTTTTCACCTCGCGAAGTGAGCGACAGGGTGCGCACCTTACTGCGTCGGGTAAAGAAGTTCTCGTCGCCGTCTCACGTCATACGTATGGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGGCAGGCACTTAACGCTGACGCGGAATGAGTTTTTATTTATTAAGACGTTACTCATGTCGCCAGTCCGCGACCGGTCCCGCCAGCCACTACTGGCAAGCGCGTGGGAAGATgc                                                                                                                                                 >  2:60120/1‑251 (MQ=255)
                                                                                                   aCCGTTTTCACCCCGCGAAGTGAGCGCCAGGGTGCGCACCTTACTGCGTCGGGCGAAACAGTTCTCGTCGCCGTCCCCCGTCCTCCGTATTGGTCCTTTTGTATTGAATGAACCCGCAGCGCAGATCCGC‑TGGTTTGACACGCAATTAACGCTGACGCGGTATGAGTTTTTCTTGCTGCAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGACGTTgcgc                                                                                                                                               >  1:443676/1‑251 (MQ=255)
                                                                                                       ttttCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAgg                                                                                                                                            >  1:159866/1‑250 (MQ=255)
                                                                                                       ttttCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAgg                                                                                                                                            >  1:542961/1‑250 (MQ=255)
                                                                                                                  cgaagttagcgcCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTa                                                                                                                                                                                                                                        >  2:503482/9‑147 (MQ=255)
                                                                                                                  cgAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCTTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑GGGTTTGACACGCCATTAACGCTGACGCGGTa                                                                                                                                                                                                                                        <  1:503482/147‑1 (MQ=255)
                                                                                                                       tgtgCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCAcc                                                                                                                           >  1:557482/1‑251 (MQ=255)
                                                                                                                       ggtgCGCCAGGGTGCTCCCCTTACTGCGTGGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCACCCGTATTGGTCATTTTGAATTGAATGACCCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCAcc                                                                                                                           <  1:4903/250‑1 (MQ=255)
                                                                                                                        gtgCGCCAGGGGGCGCTCCTTCCTGCGTCGGGTAAAGAATTTCTCGTCGCCGTCTCCCTTCACCCGTATTGGTCATTTTAAATGGAATGACCCCGCGCCGCAGATCAGC‑TGTTTTGCAACGCCATTACCGCGGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGTCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCg                                                                                                                          <  1:70950/251‑1 (MQ=255)
                                                                                                                                gggTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGGCGCGGTATGAGTTTTTATTGCTGAAGCCGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCCACTAATGGATAGCGTGTGGAAAGAGGCGGAGGGCACCTCCGATCGCACCGttgctacc                                                                                                                  >  2:21529/1‑244 (MQ=255)
                                                                                                                                        aCCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTa                                                                                                           >  1:336388/1‑250 (MQ=255)
                                                                                                                                              cTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTACCGTCATCCGGATTGGTCATTTTGGATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTACTGGATAGCGTGTGGGAAGGTGCGCAGGAAACCTACGATCGCACCGTTGATCCCCACCTTAAAACgtt                                                                                                    >  2:221253/1‑249 (MQ=255)
                                                                                                                                                           ccgccttcccggtcccggctcCCGTCATCCGTATTGGCCATTTTAAATTAAATGAACCCGGGGCGAAGACCATC‑TGGTTGTACACGCCATTAACGCTGCCGCGTTATGAGTTTTTATTGCTGAAGACGTGACTCAAGTCGCAAGGCCGGGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTg                                                                                       <  2:641338/233‑1 (MQ=255)
                                                                                                                                                                        gCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGATGCGTGCCAAGCTGCGCGCCATCAAcc                                                                          >  2:465711/1‑251 (MQ=255)
                                                                                                                                                                          cctccccctccacctttatggggttttttGAATGGAATGACCCCGGGGGGCAGA‑CAGGGTGGTTTTACCCGCCATTACCGCGGACGCGGTATGAGTTTTTATTGCGGAAGACGTTCCTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAAcccc                                                                        <  1:624125/226‑1 (MQ=255)
                                                                                                                                                                          cGTCTCCTGCCATCTGTATTGGTCATTTGGTATTGATGGACCCCGCGGCTCAGACCAGC‑TTGTTCTACACGCCATTACCTCTGCCGCGGTATGAGTGTTTCTTCCTGAAGACGTTACACAAGTCGCCTGGCCGCTTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCACcccc                                                                        <  2:214944/251‑1 (MQ=255)
                                                                                                                                                                          cGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTg                                                                                       >  1:641338/1‑236 (MQ=255)
                                                                                                                                                                               cccGTCATCCGTATGGGTCATTTTGAATTAAATGACCCCGGGGCGCAGACCAGC‑TGTTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCt                                                                   <  1:466021/251‑1 (MQ=255)
                                                                                                                                                                                    cATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCAc                                                              >  1:333231/1‑251 (MQ=255)
                                                                                                                                                                                        cGTATTGGTCATTTTGAATTGAATGACCCCGCGGCGCAGATCAGC‑TGTTTTGACACGCCATTAACGCTGACGCGGTATGATTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGAt                                                          <  1:465711/251‑1 (MQ=255)
                                                                                                                                                                                         gTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTg                                                                                       >  1:610809/1‑221 (MQ=255)
                                                                                                                                                                                                 cATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTa                                                                                                                                     >  2:267647/1‑167 (MQ=255)
                                                                                                                                                                                                 cATTTTAAATTAAATGACCCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTa                                                                                                                                     <  1:267647/167‑1 (MQ=255)
                                                                                                                                                                                                     ttgaattgaaTGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGCCGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGCTAGAGTGTGGGAAGATGCGCAGGACCCCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCACGCCATCAACCCCTATCTTTCCCCGATTAATACTCATccc                                             >  2:20422/1‑249 (MQ=255)
                                                                                                                                                                                                       aaatgaaataaaccgcggaggaaGATCAGT‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTCCTCAATTGGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCgg                                           <  1:285371/229‑1 (MQ=255)
                                                                                                                                                                                                        aattgaatGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGTTCTTTCACCGCTTAATACTCATCGCGGc                                          >  2:117721/1‑251 (MQ=255)
                                                                                                                                                                                                         attgaatGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTactaataatc                                               >  1:137616/1‑236 (MQ=255)
                                                                                                                                                                                                            aaatgaaCCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATgg                                      <  1:100626/250‑1 (MQ=255)
                                                                                                                                                                                                                      gcggcgCAGATCGGT‑GGGTTTGACACGCCATTAACGCTGACGCGGTAGGAGTTTTTATTGCTGAAGACGTTACTAAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCt                            <  1:401724/251‑1 (MQ=255)
                                                                                                                                                                                                                        tgcgcAGATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGAAATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTg                           >  2:178946/2‑250 (MQ=255)
                                                                                                                                                                                                                              gATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGcc                    >  1:297082/1‑251 (MQ=255)
                                                                                                                                                                                                                              gATCAGC‑TGGTTTGACACGCCATTAACGCTGACGCGGGATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGCTCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTAATCGCGGCATGGGATATAGCCTGAGGGGcc                    >  2:492521/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gatggtagaacacGACATTAACGCTGGCGCGGTAGGAGTTTTTCGTGAGAAAGACGTAACTCAAGTCGCAAGGCCGCGTCTGGACCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTaa               <  2:139903/242‑1 (MQ=255)
                                                                                                                                                                                                                                             acacGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCg      >  1:485272/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCAACATTAAAACGTTGCGTGCCAAGCTGCGCGACATCAACCCCGAACTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCGGTAATGCGTaacgtc    >  2:408000/1‑245 (MQ=255)
                                                                                                                                                                                                                                                 gCCATTAACGCTGACGAGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGACAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGGGCCAAGCAGCGCGCCATCGACCCCGCTCTTTCACCGATTAATACTCATCGCGGCATGGGTAATCGCCTGAGGGGCCTGAAATGCGTATCGGcg   >  2:223812/1‑249 (MQ=255)
                                                                                                                                                                                                                                                 cccATTAACGCTGACGCGTTATGAGTTTTTATTGCTGAAGACTTTACTCAGTTCGCCAGGCCTCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGTAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCAt  <  1:221253/250‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
GCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGC‑TGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCAT  >  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 < 28 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑