Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,495,194 A→T F339Y (TTT→TAT)  ahr ← NADPH‑dependent aldehyde reductase Ahr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,495,1940AT98.6% 226.6 / NA 69F339Y (TTT→TAT) ahrNADPH‑dependent aldehyde reductase Ahr
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (37/31);  minor base G (0/1);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAACTATTCTCGCGCGAATGCAGCGTCGGGCTGGGCATTTTATGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAAAAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACC  >  NC_000913/4494950‑4495430
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
gAAACTATTCTCGCGCGAATGCAGCGTCGGGCTGGGCATTTTATGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCgg                                                                                                                                                                                                                                        >  1:447618/1‑251 (MQ=255)
   aCTATTCTCGCGCGAATGCAGCGTCGGGCTGGGCATTTTATGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCtt                                                                                                                                                                                                                                     <  1:379567/251‑1 (MQ=255)
                 aaTGCAGCGTCGGGCTGGGCATTTTATGCCGGTAGCGTTGCTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGAGCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCAcgcgt                                                                                                                                                                                                                       <  2:344658/251‑2 (MQ=255)
                    gCAGCGTCGGGCTGGGCATTTTATGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTaa                                                                                                                                                                                                                    <  2:26399/251‑1 (MQ=255)
                    gCAGCGTCGGGCTGGGCATTTTATGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTa                                                                                                                                                                                                                     <  1:217552/250‑1 (MQ=255)
                    gCAGCGGCGGGCTGGGCATTTTATGCCGGTAGCGTTACTAAAAAGTCATTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTa                                                                                                                                                                                                                     <  2:529317/250‑1 (MQ=255)
                     cAGCGTCGGGCTGGGCATTTTATGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAAc                                                                                                                                                                                                                   <  1:602556/251‑1 (MQ=255)
                     cAGCGTCGGGCTGGGCATTTTATGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAAc                                                                                                                                                                                                                   >  2:410408/1‑251 (MQ=255)
                      aGCGTCGGGCTGGGCATTTTATGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGAACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGCAAGAGGGCAGCGCCTCCGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAAcc                                                                                                                                                                                                                  >  2:120144/1‑250 (MQ=255)
                        cGTCGGGCTGGGCATTTTATGCCGGTAGCGTTACAAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAAcgcg                                                                                                                                                                                                                >  1:181410/1‑251 (MQ=255)
                                          aTGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACAt                                                                                                                                                                                              <  2:8739/251‑1 (MQ=255)
                                                gTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTgg                                                                                                                                                                                         >  2:22070/1‑250 (MQ=255)
                                                gTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGa                                                                                                                                                                                        <  2:47395/251‑1 (MQ=255)
                                                 tAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGCAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGAt                                                                                                                                                                                       >  2:175845/1‑251 (MQ=255)
                                                    cGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGc                                                                                                                                                                                    <  2:255374/251‑1 (MQ=255)
                                                      ttACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGcgt                                                                                                                                                                                  >  2:325042/1‑251 (MQ=255)
                                                       tACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGcgt                                                                                                                                                                                  <  2:447618/250‑1 (MQ=255)
                                                        aCTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGcgtc                                                                                                                                                                                 <  1:441984/250‑1 (MQ=255)
                                                         cTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGcgtcgt                                                                                                                                                                               >  1:127738/1‑251 (MQ=255)
                                                                gacAGTTTAAGGCACTCGAGCGTCCCCACGCAGATGTACAGGATCTTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACACTGTCGGCAGGCTGTGCTGGCGCTACGACAAAACAGAATATGTGCGAAAGAGCGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTAAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAAtttt                                                                                                                                                                        >  2:374018/3‑251 (MQ=255)
                                                                gTCAGTTTGAGGTACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAAtttt                                                                                                                                                                        <  2:20816/251‑1 (MQ=255)
                                                                         aGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCgg                                                                                                                                                               >  2:569580/1‑251 (MQ=255)
                                                                           gCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGa                                                                                                                                                             <  1:8476/251‑1 (MQ=255)
                                                                               tGGAGCGTCCACAAGCAGCTGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAAc                                                                                                                                                           >  2:539998/1‑249 (MQ=255)
                                                                                   gCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCg                                                                                                                                                     >  1:168729/1‑251 (MQ=255)
                                                                                    cGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCgg                                                                                                                                                    <  1:539998/251‑1 (MQ=255)
                                                                                       ccACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCggtgg                                                                                                                                                 <  1:268995/251‑1 (MQ=255)
                                                                                       ccACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCggtgg                                                                                                                                                 <  1:566477/251‑1 (MQ=255)
                                                                                             gCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGgcg                                                                                                                                           <  1:86838/251‑1 (MQ=255)
                                                                                               aGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGgcgc                                                                                                                                          >  1:319886/1‑250 (MQ=255)
                                                                                                   gAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAAcc                                                                                                                                      <  1:452722/250‑1 (MQ=255)
                                                                                                       aGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTtgct                                                                                                                                 <  1:233218/251‑1 (MQ=255)
                                                                                                                 gcgcATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCa                                                                                                                       <  1:228353/251‑1 (MQ=255)
                                                                                                                   gcATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCaaa                                                                                                                     <  1:119086/251‑1 (MQ=255)
                                                                                                                         atatCAACCAAGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACGGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCa                                                                                                                >  1:552583/1‑250 (MQ=255)
                                                                                                                            tcaacagtGATATTGCTAATGACACAGACCAATGCCAGCAGGTTGTGCTGGCGATACAACAATAAAGCATATGTGCGAAAGTGGGCAGCGCTCCAGATCAGCGCTGCGAATGCTTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGAATGGCGTCGTTAATTTTCGACATCGGGAACAGATCGGTCGTCGACGCAACCTTGCTGCGGGCGACAAAACGTATc                                                                                                              >  2:287797/8‑249 (MQ=255)
                                                                                                                             caaccctgagagtGCCAATGTCACCGAACATTGGCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGTGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGACTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGc                                                                                                           <  1:354814/239‑1 (MQ=255)
                                                                                                                             caaacATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGc                                                                                                           <  1:431204/247‑1 (MQ=255)
                                                                                                                                        ttGCCAATGTCACCGAACACTGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGACTATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATAGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAACCGCATCAGCTTACGCAGCTc                                                                                                >  2:528406/1‑251 (MQ=255)
                                                                                                                                          gCCAATGTCACCGCACAATCTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAACGAGGGCAGCGCCTCAGATCAGAGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCCAAACGCATCAGCCTACGCAGCTcc                                                                                               >  2:405674/1‑249 (MQ=255)
                                                                                                                                             aaTGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATa                                                                                             >  1:41287/1‑249 (MQ=255)
                                                                                                                                                gTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCg                                                                                        >  1:230518/1‑251 (MQ=255)
                                                                                                                                                  cACCGCACAATGTCGGCAGGCTGTGCTGGCGATACGACCAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCACCACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGAGGCAAAACGCATCAGCTTACGAAGCTCATAAGGCGTg                                                                                      >  2:333207/1‑251 (MQ=255)
                                                                                                                                                      gAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATACTCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCCGCTCCTAAGGCGTGCCgg                                                                                  >  2:250522/1‑251 (MQ=255)
                                                                                                                                                          aaTGTCGGCAGGCTGTGCTGGCGATACGCCAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGc                                                                              >  2:310754/1‑251 (MQ=255)
                                                                                                                                                           atagcGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCGCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGATCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCa                                                                             <  2:459121/247‑1 (MQ=255)
                                                                                                                                                                ggCAGGCTGTGCTGGCGATACGACAAACCAGAATATGCCAGAAAGAGGGCAGCGCCTCCGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACACTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCCTCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAAcc                                                                        >  2:564368/1‑251 (MQ=255)
                                                                                                                                                                ggCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAAcc                                                                        <  2:220132/251‑1 (MQ=255)
                                                                                                                                                                   aGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAAcc                                                                        >  1:265870/1‑248 (MQ=255)
                                                                                                                                                                           cTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGc                                                             >  1:422186/1‑251 (MQ=255)
                                                                                                                                                                           atGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGCTGGCGACGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCACCCTTGCTGCGGGCGGCAAACCGCATCAGCTTACGCAGCTCATAAGGCGTGCAGGTAGCAGAACCAGAGCCGcggc                                                             >  2:350913/2‑248 (MQ=255)
                                                                                                                                                                                gATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCg                                                        >  1:386320/1‑251 (MQ=255)
                                                                                                                                                                                          aaCAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATc                                                                                                                                                                 <  1:3823/136‑1 (MQ=255)
                                                                                                                                                                                          aaCAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATc                                                                                                                                                                 >  2:3823/1‑136 (MQ=255)
                                                                                                                                                                                            cAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAAc                                            <  1:325042/251‑1 (MQ=255)
                                                                                                                                                                                                    gtgCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGc                                    >  2:292104/1‑251 (MQ=255)
                                                                                                                                                                                                    gtgCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAAGAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAgg                                     <  1:365710/250‑1 (MQ=255)
                                                                                                                                                                                                       cGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCgg                                 >  1:188265/1‑251 (MQ=255)
                                                                                                                                                                                                         aaaGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGa                                >  1:194602/1‑250 (MQ=255)
                                                                                                                                                                                                          aaGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAac                              >  1:80243/1‑251 (MQ=255)
                                                                                                                                                                                                           agagGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAac                              >  2:389056/1‑250 (MQ=255)
                                                                                                                                                                                                            gagGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAacag                            <  1:569580/251‑1 (MQ=255)
                                                                                                                                                                                                                        tCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAgcggcgggg                >  2:581550/1‑248 (MQ=255)
                                                                                                                                                                                                                           gATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGaa             <  1:408757/251‑1 (MQ=255)
                                                                                                                                                                                                                             tCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAAcc           >  1:154203/1‑251 (MQ=255)
                                                                                                                                                                                                                                      cGAATGATTTTCAATAATCGGCTTTCAACACCTCGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCaccgacc  <  1:186434/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGTaccgacc  >  2:144559/1‑251 (MQ=255)
                                                                                                                                                                                                                                          tGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAacagacag                        <  2:337932/225‑1 (MQ=255)
                                                                                                                                                                                                                                          tGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAacagacag                        >  1:337932/1‑225 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
GAAACTATTCTCGCGCGAATGCAGCGTCGGGCTGGGCATTTTATGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAAAAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACC  >  NC_000913/4494950‑4495430

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