Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 498,116 C→T P21L (CCT→CTT)  hemH → ferrochelatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913498,1160CT97.1% 231.2 / NA 68P21L (CCT→CTT) hemHferrochelatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/33);  minor base G (0/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACCTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTAC  >  NC_000913/497874‑498360
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
aGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTGTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTATTTCTTCTGCCTCACTTTTCCGCTACAATTATCAACAAGGTGAATCGATAATAGGCGGTAATGCTTCATACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCgg                                                                                                                                                                                                                                              <  1:330092/251‑1 (MQ=255)
 ggggAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGt                                                                                                                                                                                                                                             <  2:53182/251‑1 (MQ=255)
  gggAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGTCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTa                                                                                                                                                                                                                                            <  1:57662/251‑1 (MQ=255)
  gggAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAAGGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTa                                                                                                                                                                                                                                            <  1:78485/251‑1 (MQ=255)
     aaCAGGCACGCACAAACAAACTTATCTGCAATATAAGGCAGGCTATTAATTCGGCCTTTTACCACACCATCAATCTGCGCTATTCCTTTTTATGCTTTGACCTCTCACAGCCATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGAGGCAATGCGTCCGACTACAACCGGTATCCTGCTGGCAAACCTGGGTCCGCACGAAGCCCCCACACTTGAAGCggacaa                                                                                                                                                                                                                                          >  2:278502/1‑246 (MQ=255)
      ccAGGCCCGCCCAAGCCAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGTGCGATTCCTTTTTCCGATTTGACCTCTCCCATCAATTAGTTCTTCTTCCTCCCTTCGCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCCGACTAAACCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCCCACTTGAAGCGGTAAAAc                                                                                                                                                                                                                                        <  2:414645/250‑1 (MQ=255)
           cccGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTAt                                                                                                                                                                                                                                   <  2:403633/251‑1 (MQ=255)
                     cAAACTTCTCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCCGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCCCACTTGAAGCGGTAAAAAGCTATCTGAAACAAt                                                                                                                                                                                                                         >  2:275635/1‑251 (MQ=255)
                      aaaCTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGACAGACTAAAAACGGTATCCTGCTGGCAAACCTGGGTACGCCCGCTGCCCCCACACTTGAAGCGGTCAAACGCTATCTGCAACAAtt                                                                                                                                                                                                                        >  2:620035/1‑251 (MQ=255)
                        aCTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAAtttt                                                                                                                                                                                                                      >  2:559138/1‑251 (MQ=255)
                              cAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAAACTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATACGAGGCGGTAATGCGTCAGACTAAAACCGGTATACTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCg                                                                                                                                                                                                                >  1:651009/1‑251 (MQ=255)
                                gCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGAc                                                                                                                                                                                                              <  1:728483/251‑1 (MQ=255)
                                  aaTCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAg                                                                                                                                                                                                            <  2:526666/251‑1 (MQ=255)
                                  aaTCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAAAAATTTTTAAGCGAc                                                                                                                                                                                                              >  2:385538/1‑249 (MQ=255)
                                  aaTCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAAGCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTTACTTTTCCGCTACAATGATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAATCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAGCGCTATCTGAAACAATTTTTAATCGACAg                                                                                                                                                                                                            <  1:401100/251‑1 (MQ=255)
                                       cAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGAcgc                                                                                                                                                                                                        <  1:347450/250‑1 (MQ=255)
                                          gCCGGATATTCATTCGGCCTTTTACAAACACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCCATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTCTCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGCAGCGGGAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTgg                                                                                                                                                                                                    >  2:466473/1‑251 (MQ=255)
                                             ggATATTCATTCGGCCTTTTCCAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCCATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATTCCCCCACACTTGAAGCGGTAAAACGCTCTCTGAAACAATTTTTAAGCGACAGACGAGTGGTTg                                                                                                                                                                                                 >  1:623832/1‑251 (MQ=255)
                                                 attcattcGGACTTTTACAAACACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACATCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGAAAGACGCGTGGTTGATAc                                                                                                                                                                                             >  2:654822/1‑251 (MQ=255)
                                                          gCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAAGTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCGGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGgtt                                                                                                                                                                                    <  1:66554/251‑1 (MQ=255)
                                                              tttACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGgttgt                                                                                                                                                                                  <  1:548361/249‑1 (MQ=255)
                                                                 aCAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTTtgg                                                                                                                                                                             >  2:348772/1‑251 (MQ=255)
                                                                   aaaaaCCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTAtggtg                                                                                                                                                                           <  1:325961/251‑1 (MQ=255)
                                                                           cAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCTCTCGTGAAGCGTTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCAtt                                                                                                                                                                     <  1:644191/249‑1 (MQ=255)
                                                                                        tCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGAt                                                                                                                                                      >  1:131567/1‑251 (MQ=255)
                                                                                                 tGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCCCTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGAACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATTGTGGCCATTGCTGAGCGGCGTGATTTTGCCGCt                                                                                                                                             >  1:389648/1‑251 (MQ=255)
                                                                                                   aTTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTTGCCATGGCTGCGCGGCGTGATTTTGCCGCTgc                                                                                                                                           <  2:593908/251‑1 (MQ=255)
                                                                                                              cacaGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTcgccgc                                                                                                                                  <  1:559138/249‑1 (MQ=255)
                                                                                                                  gCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTgg                                                                                                                            <  2:297181/251‑1 (MQ=255)
                                                                                                                  gCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTgg                                                                                                                            <  2:60154/251‑1 (MQ=255)
                                                                                                                   cAATTAGTTCTTCTTCCTCGCTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGc                                                                                                                           <  2:63496/251‑1 (MQ=255)
                                                                                                                       tAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGCTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGAAAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCCTTGCTGCGCGGCGAGTTTTTGCCGCTGCGCTCGCCGCGTGTGGCGa                                                                                                                         >  1:410844/1‑249 (MQ=255)
                                                                                                                        aGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAg                                                                                                                       <  1:324788/250‑1 (MQ=255)
                                                                                                                         gTTCTTCTTCCTCACTTTTCCGCTACAATTATCACCAAGTTGAATCGATCCGAGGCGGTACAGCGTAAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCTGGATGTTGCCGCTGCGCTCGCCGCGTGTGGCggaacc                                                                                                                     >  2:341767/1‑246 (MQ=255)
                                                                                                                            cttcttcCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGt                                                                                                                   >  2:145823/1‑250 (MQ=255)
                                                                                                                              tcttcCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGAATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATg                                                                                                                >  1:68800/1‑251 (MQ=255)
                                                                                                                               cttcCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGc                                                                                                               >  1:161862/1‑251 (MQ=255)
                                                                                                                                 tcCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGcc                                                                                                              <  1:385538/250‑1 (MQ=255)
                                                                                                                                  ccTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCctt                                                                                                            >  1:399359/1‑250 (MQ=255)
                                                                                                                                     cACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTACAACCGGCATCCTGCCTGCAAAACTGGGCACGCCCGCTGACCCCACACTTGAAGCGGTAAAACCCTATCTGAAACAATTTTTAAGCGCCAGAAGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGGGTGCTTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGTCtctgt                                                                                                         >  1:48443/1‑251 (MQ=255)
                                                                                                                                     cACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCctct                                                                                                           <  1:499685/249‑1 (MQ=255)
                                                                                                                                                      tATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGTTACGCCCGATGCCCCCACGCTGGAAGCGGTAAAACGCTATCTGAAACAAGTTTTAAGCGACAGACGCGTGGTTGATACCTCACGTTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTGTGCCGCTGCGCACGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAAGTGGc                                                                                         <  2:54359/250‑1 (MQ=255)
                                                                                                                                                      tATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGGGTGATTTTGCCGCTGCGCTCGCCGCGTGGGGCGAAGCTGGATGCCTCTGTCGGGTTGGAAGGTGGcc                                                                                        >  2:603317/1‑250 (MQ=255)
                                                                                                                                                                      cGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGGGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATTCCTCTGTCTGGATGGAAGGTGGCGCTCCGCTGATGGttt                                                                         <  2:389648/250‑1 (MQ=255)
                                                                                                                                                                       gATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTAc                                                                       >  1:54674/1‑251 (MQ=255)
                                                                                                                                                                            gagGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGAAAACCTGGGCACGCCCGATGCCCCCACACTTGAAGCGGTAACACGCTATCTGAAACACTTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAgccg                                                                  >  1:697646/1‑251 (MQ=255)
                                                                                                                                                                              ggCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAgccgcc                                                                <  1:196307/251‑1 (MQ=255)
                                                                                                                                                                                 ggTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGc                                                             <  2:161862/251‑1 (MQ=255)
                                                                                                                                                                                    aaTGCGTCAGACTAAAACCGATATCCAGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGATAAAACGCTAGCAGAAACAATGTTTAAGCGACAGAGGCGTGGTTGATACCTTACGGTTGTTATGGTGGCCATTGCAGCGTGACGGGATTTTGCCGCTGCGTTCGCGGCGTGTGGCGATGATGTATGCTTGTCTGTGTATTTAAGGTGGCTCGCCGCTGATGGTTTACATCCGCCAGCAAc                                                          <  1:463532/251‑1 (MQ=255)
                                                                                                                                                                                     aTGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAAc                                                          >  1:349822/1‑250 (MQ=255)
                                                                                                                                                                                       gCGTCAGAATAAAACCGGTATCCTGCTGGCAAACCTGAGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAAAGCTATATGAAACACTTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGGGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGCCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAAcagc                                                       >  1:240199/1‑251 (MQ=255)
                                                                                                                                                                                         gTCAGACTAAAACCGGTCTCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAAcagc                                                       >  1:221366/1‑249 (MQ=255)
                                                                                                                                                                                           cAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGgc                                                   <  1:90126/251‑1 (MQ=255)
                                                                                                                                                                                              atgaaaaCCGGTATCCTGCTGGCAAACCTGGGTACGCCCTATGCCCCCACACTTGAAGCGGTAAAACTCTAGCTGAAACAATTTTTAATCGACAGACGCGTGGTTGATACCGCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAATCTGTATGCCTCTTCCTGTATGGAAGGTGTCTCGCCGCTGATGGTTGACAGCCGTCAGCAACAGCAGGCGCt                                                <  2:651009/248‑1 (MQ=255)
                                                                                                                                                                                                    aCCGGTATCCTGCTGGCAAACCTGGGTCCGCCCGATGCCCCCACACTTGCAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGcac                                           >  2:114589/1‑250 (MQ=255)
                                                                                                                                                                                                    aCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGcac                                           >  2:552382/1‑250 (MQ=255)
                                                                                                                                                                                                       ggTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCCCGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAAcagca                                                      >  2:41175/1‑236 (MQ=255)
                                                                                                                                                                                                       ggTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAAcagca                                                      <  1:41175/236‑1 (MQ=255)
                                                                                                                                                                                                       ggTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACg                                       >  1:544595/1‑251 (MQ=255)
                                                                                                                                                                                                                    cAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATg                          >  1:525902/1‑251 (MQ=255)
                                                                                                                                                                                                                      aaCCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGcc                        <  1:114589/251‑1 (MQ=255)
                                                                                                                                                                                                                       aCCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGccc                       <  1:348772/251‑1 (MQ=255)
                                                                                                                                                                                                                        ccTGGGTACGCACCACGCCTCCACACTTGAAGCGCAAACCCCCTATCTGAAACAATTTTTAAGCGCCAGACGCGTGGTTGATACCTCACGCTTGCTATGGTGGCCATCGCTGCGCGGCGTGATTTTGCCGCTGCGATCGCCGCGTGTGGCGCAGCTGTATGCCTCTGTCTGGTTGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGCTTACCCGAGATGTCCg                      >  2:120937/1‑251 (MQ=255)
                                                                                                                                                                                                                           gggTACGCCCGATGCCCCCACACTTGAAGCGGTAAACCGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATCCCTCGCGGTTGTTAGGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCTGAGATGCCCGTAg                   <  2:544595/251‑1 (MQ=255)
                                                                                                                                                                                                                                cGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTg              >  1:233842/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cccgaggaccccACACGTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTAGGGTGGCCATTGCTGCGCGGCGTTATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTggg            <  2:221366/243‑1 (MQ=255)
                                                                                                                                                                                                                                    cGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCAGTAGCGCTGGGa           >  2:221695/1‑250 (MQ=255)
                                                                                                                                                                                                                                            ccACACTTGAAGCGGTAAAACGCTATCTGAAACAATTGTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACATCTTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTAc  <  1:664887/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
AGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACCTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTAC  >  NC_000913/497874‑498360

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