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,1680GT96.6% 167.6 / ‑3.8 61G14G (GGC→GGA) ahrNADPH‑dependent aldehyde reductase Ahr
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (30/27);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

ACAGCTACGCGCCGTCCAGCCAATCCCGACACGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTG  >  NC_000913/4495922‑4496413
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
acaGCTACGCGCCGTCCAGCCAATCCCGACACGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTctccgc                                                                                                                                                                                                                                                   >  1:456164/1‑251 (MQ=255)
  aGCTACGCGCCGTCCAGCCAATCCCGACACGCTGACCGACCTGCACACCTTTCTCCTGCGCGGCGCTCCCGCGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCACCGACAGATCGGAATGGCAGGTCCCGCAGTAATCCACATGCACTTCCACCTCTTGTGGCCTCCGCTCACCGGGATCGTCCTCATAAACTTCCCGTTctccgacc                                                                                                                                                                                                                                                 >  2:469507/1‑248 (MQ=255)
  aGCTACGCGCCGTCCAGCCAATCCCGACACGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTccgccc                                                                                                                                                                                                                                                 >  1:4331/1‑251 (MQ=255)
     tACGCGCCGTCCAGCCAATCCCGACACGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGct                                                                                                                                                                                                                                              >  1:248111/1‑251 (MQ=255)
       cgcgCCGTCCAGCCAATCCCGACACGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGcttc                                                                                                                                                                                                                                            >  1:257707/1‑251 (MQ=255)
         cgccgcCCAGCCAATCCCGACACGCTGACCGACCTGCAAACCTTTAGCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGcttctt                                                                                                                                                                                                                                          <  2:407215/251‑1 (MQ=255)
                           gACACGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATAACCTCATGCACGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTCATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCATCAGCTCAACGGGATCGTATTCATAAACTTCCAGTTCTACGCCCGCTTCTTTTGCGGCATAGCttttt                                                                                                                                                                                                                         >  1:278890/1‑250 (MQ=255)
                              ccGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCCCCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTatca                                                                                                                                                                                                                     <  2:347060/250‑1 (MQ=255)
                               cGCTGACCGACCTGCAAACCTTGATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTatcat                                                                                                                                                                                                                    <  1:99/251‑1 (MQ=255)
                                       gaccTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACAttt                                                                                                                                                                                                            >  2:76403/1‑251 (MQ=255)
                                               aaCCTTTATCCCGAGCGGCGCCCCCGAGGGCCACCACGCGCCCAACCCCCCCCTGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTcc                                                                                                                                                                                                    <  2:394596/251‑1 (MQ=255)
                                                aCCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCt                                                                                                                                                                                                   >  2:413850/1‑251 (MQ=255)
                                                  cTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCCTAGCTTTTTATCATCGACATTTTTTGGTCCTtc                                                                                                                                                                                                 >  2:1210/1‑251 (MQ=255)
                                                      aTCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTg                                                                                                                                                                                             <  1:412610/251‑1 (MQ=255)
                                                           gcgcGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGgtg                                                                                                                                                                                          <  1:76403/249‑1 (MQ=255)
                                                               ggcgCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGCTCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCCCTTCAACATCTTGTGGCCTCAGCTCACCGGGCTCGTACTCCTAAACTTCCCGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATAATCGCCATTTTTTGGTCCTTCTATGGTGTTGttt                                                                                                                                                                                    >  2:294527/1‑251 (MQ=255)
                                                                gcgcTCCCGTGTGCCAACACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGCTATTGTGAAAATCCCCATTCGTTATCGATCTTCGACAGCTCGGAATGTCAGGACCCGCAGTAATCCACCTGCACTTCAAAATCTTGTGGCCTCAGCTCACCGGGATCGTACTACTAAACTTACAGTTCTCCGCCCGCTACTTTTGCGGCATAGCTTTTTATCATCAACATTTTTTGGTCCTTCTTTGGTGTTGTTTg                                                                                                                                                                                   >  2:294185/1‑251 (MQ=255)
                                                                    ccccggggtcccccacGCGCCCAATCCCCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTGAGCGATCATCGACAGATCGGAATGGCATATCCCGCAGTAATCCCCCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGAGCGTACTCATAAACTCCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTGCTCTGGTGTTGTTTGGGaa                                                                                                                                                                               <  2:436851/239‑1 (MQ=255)
                                                                     cccGAGTGCCACCACGCGCCCAATCACCTCAGGCCCTGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCTAGCACCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGGGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGa                                                                                                                                                                                <  2:433384/249‑1 (MQ=255)
                                                                                      gcCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAATTGTAGAGCATGGCAggg                                                                                                                                                             >  2:84087/1‑251 (MQ=255)
                                                                                        cccctcgccTCATGCGCGGCAACCAGCGGATATTGTGAAAAACCCCATTCGTTATCGATCATCGCCATACCGGAAGGGCAGAGCCCGCAGTAAGCCACCTGCAGTCCAACAACTTGTGGCCGCAGCCCACCGGGATCGTACTCATAAACTACCAGTTCCCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCAGCGACAGTTTTTGGGCCTACGCTGGTGTTGTTTTGGAAGTGTAGAGCATGGCAGGGCg                                                                                                                                                           <  2:475884/244‑1 (MQ=255)
                                                                                           aTCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGAgggg                                                                                                                                                        >  1:363477/1‑251 (MQ=255)
                                                                                             ccccccagtcccgtaaaCCAGCGGATATTGTGAAAATCCCCATTCGTTACCGACCATCGCCCGATCGGAATGGCATATCCCGCAGTAATCCACCGGCCCTTCAACATCTTGTGGCTTCAGCTCACCGGGCTCGTATTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTAGCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTg                                                                                                                                                      <  2:274155/236‑1 (MQ=255)
                                                                                                 ccATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTAACGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTg                                                                                                                                                  <  2:4331/250‑1 (MQ=255)
                                                                                                                   ggaatattgaaaaaTCCCCATTCTTTATCTACCATCGACATATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACAGCTTGTGCCCTCAGCTCACCGGGATCGTCCTCAGAATCGTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTGTTGGTCCTTCTCTGGTGTTGTTTTGGAAGTGTAGAGCATGGCAGGGCGGGGTTCCTGGAGGGTGACAAAGGTTac                                                                                                                                <  2:409877/241‑1 (MQ=255)
                                                                                                                            gAAACTCCCCATTCTTTATCGATCATCGACAGATCTGAAGGGCAGATCTCGCAGTAATCCATCTGCACTTCACCACCTTTTGGCCTCAGCTCACCGGGATCGTACTCATAATCTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGTAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGGTACACATCGCTg                                                                                                                       <  1:311992/251‑1 (MQ=255)
                                                                                                                               aaTCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGGTCTCCGCCCGCTTCGTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTTTATAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTAt                                                                                                                    <  1:469507/251‑1 (MQ=255)
                                                                                                                                aTCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTAtt                                                                                                                   >  1:383247/1‑250 (MQ=255)
                                                                                                                                aTCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTATTTGGGAAGTGTAGAGCATGGCAGGGCGGGGGGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATg                                                                                                                   >  1:141454/1‑251 (MQ=255)
                                                                                                                                                gATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACAtttttt                                                                                                                                                                                                         >  1:86250/1‑149 (MQ=255)
                                                                                                                                                   catcGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGc                                                                                                >  1:389492/1‑251 (MQ=255)
                                                                                                                                                      cGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCGGCATTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAGCTGCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCGTTTTATCATCGACATTGTTTGGGCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCAc                                                                                              <  1:22360/250‑1 (MQ=255)
                                                                                                                                                          aGATCGGAATGGCAGACCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGc                                                                                         <  2:196097/251‑1 (MQ=255)
                                                                                                                                                              cGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGCGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAgg                                                                                     >  2:309180/1‑251 (MQ=255)
                                                                                                                                                              cGGAAGGGCAGTTCCCGCAGTAACCCACCTGCACTTCAACAGTTTGTGGCCTCAGCTCAACGGGATCTTACTCATAAACTCCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGat                                                                                      <  2:4231/250‑2 (MQ=255)
                                                                                                                                                                     gCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTTGCGAGGTGCGAGc                                                                              >  2:433800/1‑251 (MQ=255)
                                                                                                                                                                       agaccccGCAGTAATCCACCGGTCCTTCAACATCTTGTGGCCCCAGCGCCCCGGTATCGTATACATAAACTTTCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTAACATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCaa                                                                            <  2:132570/247‑1 (MQ=255)
                                                                                                                                                                       aGATCCCGCAGTAATCCACATGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGGGCCTGGAGTGTGAAAAAGGTTACACATAGCTTGATGCACAGCTGAAAATTTCCGAACTTAGCGAGGGGCGAGCa                                                                             >  2:379439/1‑250 (MQ=255)
                                                                                                                                                                          tCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTTCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCCATGCTGAAAATTTCAGCACTTAGCGAGGGGCGAGCAAGCt                                                                         >  2:112060/1‑251 (MQ=255)
                                                                                                                                                                                   aaTCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGc                                                                >  1:136192/1‑251 (MQ=255)
                                                                                                                                                                                   aaTCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGc                                                                >  1:176860/1‑251 (MQ=255)
                                                                                                                                                                                    aTCCACCTGCACTTCAACATCTTGAGGCCTCAGCTCACCGGGGTCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATAATCGACATTTTTTGGTCCTTCTCTGGTGTTGATTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGTGTACACATTGCTGTATGCAAGGCTGCAAATTTCAGCACTTTGCGCGGTGAGGGCGAGCGGGCGCtttc                                                                >  2:69815/1‑248 (MQ=255)
                                                                                                                                                                                       cACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCAtgg                                                            >  1:247989/1‑251 (MQ=255)
                                                                                                                                                                                            gcccTTCAACATCTTGTGGCCTCAGCTCCCCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGc                                                        <  2:136192/247‑1 (MQ=255)
                                                                                                                                                                                            gCACTTCAACATCTTGTGGACTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGTCAAAGGTTGCACATCGCTGGATGCAATGCTGAAAATTTCAGCACTTAGCGAGGGTCGAGCAAGCTGGCGCTTGAATGGTGGCg                                                       >  2:259266/1‑251 (MQ=255)
                                                                                                                                                                                                        cTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGtat                                           >  2:343804/1‑251 (MQ=255)
                                                                                                                                                                                                           gtgGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAAt                                        <  1:309180/251‑1 (MQ=255)
                                                                                                                                                                                                            tgGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATc                                       >  1:35741/1‑251 (MQ=255)
                                                                                                                                                                                                                  cAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCGTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGAGCCTGGAGTGTGACAACGGGTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGTAAGCTGGCGCTTGCATGGAGGCGTGCTACAGGCATAATCCACAAc                                 >  1:358669/1‑251 (MQ=255)
                                                                                                                                                                                                                    gCTCACCGGGATCGTACTCATAAACTTCCAGGTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGt                               <  2:248111/251‑1 (MQ=255)
                                                                                                                                                                                                                     cTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTg                                                     <  1:425269/228‑1 (MQ=255)
                                                                                                                                                                                                                     cTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTg                                                     >  2:425269/1‑228 (MQ=255)
                                                                                                                                                                                                                     cTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGtt                              >  1:293614/1‑251 (MQ=255)
                                                                                                                                                                                                                      tCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGGGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAAAGttt                             >  1:221013/1‑251 (MQ=255)
                                                                                                                                                                                                                         ccGGGATCGTACTCAAAAACTCCCGGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTcc                          <  2:389492/251‑1 (MQ=255)
                                                                                                                                                                                                                          cGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCg                         <  1:104135/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gTACTCAGAAACTTCCATTTTCCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCAGCGACATTTTTTGGTCCTTCTCTGGTGTTGTGTGGGAAGTGTAGAGCATGGCAGGGCGGGGGGCCTGGAGGGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTGTTCCGCATACct                  <  2:261499/251‑1 (MQ=255)
                                                                                                                                                                                                                                  tACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATa                     <  2:257707/247‑1 (MQ=255)
                                                                                                                                                                                                                                          agcTTCCCGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTATTTACCATCGACATTTTGTTGTCCTTCTCTTGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTGACACATCGCTGTATGCAATGCTGAAAATTTCAGCCCTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGc         <  2:239458/249‑1 (MQ=255)
                                                                                                                                                                                                                                                cAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGt   <  2:456164/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 aGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCa                                                                                                  >  2:423687/1‑155 (MQ=255)
                                                                                                                                                                                                                                                 aGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCa                                                                                                  <  1:423687/155‑1 (MQ=255)
                                                                                                                                                                                                                                                 aGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTg  <  1:475188/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
ACAGCTACGCGCCGTCCAGCCAATCCCGACACGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTG  >  NC_000913/4495922‑4496413

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: ‑