Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,496,168 G→T G14G (GGC→GGA)  ahr ← NADPH‑dependent aldehyde reductase Ahr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,496,1680GT98.6% 209.5 / ‑3.8 73G14G (GGC→GGA) ahrNADPH‑dependent aldehyde reductase Ahr
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (34/37);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

GCTACGCGCCGTCCAGCCAATCCCGACACGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCG  >  NC_000913/4495925‑4496408
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
gCTACGCGCCGTCCAGCCAAACCCGACACGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTccgcccg                                                                                                                                                                                                                                           <  1:252305/251‑1 (MQ=255)
      cgcCGTCCAGCCAATCCCGCCCCTCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCGCCCGAGTCCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTACCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCC‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGcttc                                                                                                                                                                                                                                       <  2:146259/249‑1 (MQ=255)
        ccctCCACCCAATCCCGACACGCTTACCGACCTGAAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCtttt                                                                                                                                                                                                                                   <  2:14710/248‑1 (MQ=255)
              acccAAACCCGTCACGCCGACCGACCTGCAACCCTTTATCCTGCGCGGCGCACCCGAGTGCCACCACGCGCCCAATCACCCCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCa                                                                                                                                                                                                                             <  2:523384/249‑1 (MQ=255)
               cccAATCCCTACACGCTGGCCGACTTGCAAACCTGTATCCTGCGCGGATATCCCGAGTGACTCCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAATACCCCAGTCGTTATCGATCAGCTACAGATCGGAATGGCAGATCCCGCAGTAAGCCA‑CCTGCACTCCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCTCCCGCTTCTATTGCGGCAt                                                                                                                                                                                                                            <  2:516178/250‑1 (MQ=255)
                     cccGACACGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCttt                                                                                                                                                                                                                      <  1:137251/251‑1 (MQ=255)
                          cacGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTatc                                                                                                                                                                                                                 >  1:46303/1‑251 (MQ=255)
                              cctccccgcctctcttccgatccccTGCGCGGCGCCCCCGAGTGCCACCCCGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGTTCACCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCg                                                                                                                                                                                                             <  2:532818/228‑1 (MQ=255)
                              cTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCg                                                                                                                                                                                                             >  1:92119/1‑251 (MQ=255)
                                gaccgaccTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGAc                                                                                                                                                                                                           <  1:263667/251‑1 (MQ=255)
                                         gCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGTCAACCAGCTGATATTGTGAAAACCCCCATTCGTTTTCGATCATCGACAGATCGGAATGGCAGACCCCGCAGTAATCCC‑CCTGCACTTCAACACC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACAttttta                                                                                                                                                                                                    <  2:110919/249‑2 (MQ=255)
                                                 ttATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTtctc                                                                                                                                                                                          >  1:399115/1‑251 (MQ=255)
                                                 tgtcccgtcgctgtgcgcctgAGTGCCACCACGCGCCCAATCCCCGCATGCCCGTCAACCAGCGGATCTGGGGAACACCCCAAGTCGTTTGCGATCATCGACAGATCGGAATTGCAGACCCCGCAGTAATCCA‑CCTTCACTTCAACAGC‑GTGTGGCCTCAGTTCGGCGGGAGCGTACTCATAA‑CTGTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTTGTCTTtct                                                                                                                                                                                           <  2:324948/230‑1 (MQ=255)
                                                   aaccTGCGCGGCGCTCCCGAGGGCCACCCCGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTg                                                                                                                                                                                        <  2:304420/249‑1 (MQ=255)
                                                     ccTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCg                                                                                                                                                                                                             >  1:532818/1‑228 (MQ=255)
                                                         cgcgGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGtgttg                                                                                                                                                                                  <  1:146868/251‑1 (MQ=255)
                                                               ccTCCCGTGTTCCACCACGCGCCCAATCACCTCTTGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGAGCATCTACAGATCGGAATGGCAGATCCCGCATTAATCCC‑CCTGCCCTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATTGTTCTCATAATCT‑TCCAGTTCTCCTCCCGCTTCTTTTGCTGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTggg                                                                                                                                                                            <  2:490133/250‑1 (MQ=255)
                                                                  ccagAGTGCCACCACGCGCCCAAGCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGTTCACCGGGATAGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTatcatc                                                                                                                                                                                                              <  1:607558/211‑1 (MQ=255)
                                                                           cacaacGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGACGTGTAGAGc                                                                                                                                                                 >  2:567593/5‑250 (MQ=255)
                                                                            accacGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATGGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCAt                                                                                                                                                               <  1:248143/251‑1 (MQ=255)
                                                                               ccccgcCCAATCCCCTCAGGCCCGGCACCCAGCGGATAATGTGAAAACCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCGCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTGTATCATCGACATTTTGTGGTCCTTCTCTGGTGTTGTTTGGGAAGGGTAGAGCATGGc                                                                                                                                                            <  2:46303/248‑1 (MQ=255)
                                                                                  cgcCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACAtt                                                                                                                                                                                                        >  1:414288/1‑204 (MQ=255)
                                                                                   gcCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAggg                                                                                                                                                        <  1:305292/251‑1 (MQ=255)
                                                                                        aTCACCTCATGCCCGGCAACCAGCGGATATTGTGAAACTCCCCATTCGGTATCTATCATCGACAGATCGGAATGGCAGATCTCGCAGTAATCCT‑CCTGCACTTCAACTCC‑TTGTGTCCTCAGCTCTCCGGGATCGTACTCATAATCT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGTGAAGTGTAGAGCATGGCTGGGCgggg                                                                                                                                                   <  1:23835/251‑1 (MQ=255)
                                                                                         tCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTCTCGCTCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGCTCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCggggg                                                                                                                                                  >  1:283235/1‑250 (MQ=255)
                                                                                          ccccTCATGCCCGGCAACCAGCGTATATTGTGAAAATCCCCATTCTTTATCGATCATCGACAGATCGGAAGGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCTCCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATGTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTg                                                                                                                                                 <  2:256443/249‑1 (MQ=255)
                                                                                                  gCCCGGGAAACAGCGGATATTGTGAAAATCCCCATTCGTTTACGACCACCGACAGATCGGAATGGACGATCCCGCAATAATCCC‑CCTGCACTTCAACATT‑TTGTGGCCCCAGCCCACCGGGGTCGTACTCAAAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTAGCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAgt                                                                                                                                         <  2:282840/251‑1 (MQ=255)
                                                                                                  gCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGAAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCTGTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTgggtg                                                                                                                                         >  2:648568/1‑248 (MQ=255)
                                                                                                              gCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATAATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCGTGGCAGGACGGGGTGCCTGGAGTGTGACAAAGGtt                                                                                                                             >  2:468204/1‑251 (MQ=255)
                                                                                                                ggATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCG‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTTCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTac                                                                                                                           <  1:282009/251‑1 (MQ=255)
                                                                                                                ggATATTGAGAAAATCACCATTCGTTATCGGTCATCGACAGATCGGCATGGCTGATACCGCAGTAATCCA‑ACTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACA‑TACATTTCTACGCCGGCTTCTTTTTCGGCATAGCTTTTTATCATCGACATTTTTTGTTCCTTCTCTGTTGTGGTTTGGGAGGAGTAGCGCATGACTGGACGGGGCGCCTgcactgggaataatgttat                                                                                                                           >  2:195131/1‑233 (MQ=255)
                                                                                                                           aaaTCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTg                                                                                                                                                                                                   >  1:240153/1‑168 (MQ=255)
                                                                                                                           aaaTCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTg                                                                                                                                                                                                   <  2:240153/168‑1 (MQ=255)
                                                                                                                                  cATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATg                                                                                                         <  1:239478/251‑1 (MQ=255)
                                                                                                                                   aTTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACCTCACTGTATGCACTGc                                                                                                        >  2:256092/1‑251 (MQ=255)
                                                                                                                                      cGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGa                                                                                                     <  1:183125/251‑1 (MQ=255)
                                                                                                                                           tCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAAtt                                                                                                >  1:564948/1‑251 (MQ=255)
                                                                                                                                             gATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTc                                                                                              >  1:413957/1‑251 (MQ=255)
                                                                                                                                              atcatcGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATTGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTc                                                                                              >  1:183225/1‑250 (MQ=255)
                                                                                                                                                catcGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAg                                                                                            <  1:636659/250‑1 (MQ=255)
                                                                                                                                                 atcGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGc                                                                                           >  1:599090/1‑250 (MQ=255)
                                                                                                                                                 actgaaagagcggaagagCAGATCCCGCAGTAA‑CCAGCCGTCCCTTCAACATCGTTG‑GGCCACAGCTCACCTGGATCGTACTCATAAACT‑TCCATTTCTCCGCGCGCTTCTTTTGCGGCATAGCTTGTGATCATCGACATTTTTTGGTCCTACTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCa                                                                                          <  2:352041/234‑1 (MQ=255)
                                                                                                                                                   cGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCAc                                                                                         <  1:39451/250‑1 (MQ=255)
                                                                                                                                                      cAGATCGGAATGGCAGATCCCTCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAg                                                                                     <  1:563628/251‑1 (MQ=255)
                                                                                                                                                      cAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCCACATC‑TTGTGGCCTCAGCTCCCCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGAAAGGCGGGGGGCCTGGAGTGTGCGAAAGGTTACAAATCGTTGTATGCAATGCTGAAAATTTCAGAACttcg                                                                                     >  2:156602/1‑249 (MQ=255)
                                                                                                                                                       aGATCGGAATGGCAGATCCCGCATTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCCCCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGc                                                                                    <  1:433802/251‑1 (MQ=255)
                                                                                                                                                         aTCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACCCATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCg                                                                                   >  2:543135/1‑250 (MQ=255)
                                                                                                                                                          tCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAAAATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCGTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCg                                                                                   >  1:232781/1‑249 (MQ=255)
                                                                                                                                                                tGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCg                                                                                                                     <  1:251254/209‑1 (MQ=255)
                                                                                                                                                                tGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCg                                                                                                                     >  2:251254/1‑209 (MQ=255)
                                                                                                                                                                       tCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCt                                                                    >  1:548055/1‑251 (MQ=255)
                                                                                                                                                                         ccGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTgg                                                                  >  1:261461/1‑251 (MQ=255)
                                                                                                                                                                             aGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCt                                                              >  1:141656/1‑251 (MQ=255)
                                                                                                                                                                              gTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCtt                                                             >  1:545277/1‑251 (MQ=255)
                                                                                                                                                                               gtaatctcccTGCACGTCAACAGC‑TGTTGGACTCAGCTCACCGGGATCGTACCTAGAAACT‑TTCAGTTCCCCCCACGCTAGTTGTGCGGCATAGCTGTTTATCATCGACATTTTTTTGTTCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCTGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACAGCGCTGTGTGCAATGTTGAAAATTTCAGCACTTAGCGAGGTGCGAGCGAGCTGgcgcgt                                                             <  2:632079/243‑3 (MQ=255)
                                                                                                                                                                                aaTCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGTTCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCGTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGc                                                           <  1:640670/251‑1 (MQ=255)
                                                                                                                                                                                 aTCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCT‑CTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCCATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCAt                                                         >  2:247819/1‑251 (MQ=255)
                                                                                                                                                                                   ccA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCGCCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGAGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGACGCTTGCAtg                                                        >  1:162460/1‑251 (MQ=255)
                                                                                                                                                                                      cccTGCTCTTCACCATC‑TTTGGGCCGCAGTTCTCCGGGATCGTACTACTATCCT‑TCCATTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTGTTGGTCCTTCTCTGGTGTTGTTTGGGAAGGGGAGAGCATGGCAGGGCGTGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAAGTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCAtggt                                                      <  2:400227/250‑1 (MQ=255)
                                                                                                                                                                                              ttCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAGGTGTAGAGCATGGCAGGGAGGGGTGCCTGGAGTATGACATAGGTTACACATAGCTGTATGAAATGCTGAAAATTTCAGCACTTCGCGAGGTGAGCGAAAGCTGGCGCTTGCAtggtggatgggt                                              >  2:25034/1‑245 (MQ=255)
                                                                                                                                                                                                  aCATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGAc                                            >  1:290329/1‑249 (MQ=255)
                                                                                                                                                                                                               cTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCa                                                                                                                                                                <  1:480247/121‑1 (MQ=255)
                                                                                                                                                                                                               cTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCa                                                                                                                                                                >  2:480247/1‑121 (MQ=255)
                                                                                                                                                                                                                 cAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCcaca                              <  1:623176/249‑1 (MQ=255)
                                                                                                                                                                                                                            taTCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGACCTTCTCTGGTGTTGTTTGGGA‑GTTGAGAGCATGGCTGGGCGGGGTGCCTAGAGTGTGGCAGCGGTTACACCTCGCTGTATGCAATGCTTAAACTTTAAGCACTTACCGAGGTGCGGGGCAGCTGGCTCTTGCATAGTGGCGTGCGACAGGTATAATccccacagttttcagaatc                >  2:189765/2‑234 (MQ=255)
                                                                                                                                                                                                                            gATCGTACTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCACTGCTGAAAATTTCAGCACTTAGCGAGGAGCGAGCAAGCTGGCGCTTGCATGGAGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCAt                 >  2:850/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cTCATAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCtt          >  1:106659/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cTCATAAACT‑TCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCtt          >  2:322590/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aTAAACT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGCGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTc         >  1:391258/1‑249 (MQ=255)
                                                                                                                                                                                                                                        aaaCT‑TCCATTTCTCCGCCCGCTTCGTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTg     <  2:141656/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aaaCT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAATCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGt      <  1:495959/250‑1 (MQ=255)
                                                                                                                                                                                                                                        aaaCT‑TCCAGTTCTCCGCCCGCTTCTTGTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTATGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGt      <  1:25034/250‑1 (MQ=255)
                                                                                                                                                                                                                                           cT‑TCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCg  >  1:130044/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
GCTACGCGCCGTCCAGCCAATCCCGACACGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCA‑CCTGCACTTCAACATC‑TTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACT‑TCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCG  >  NC_000913/4495925‑4496408

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