Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,646,970 G→T Q224H (CAG→CAT)  gor → glutathione reductase (NADPH)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,646,9700GT93.2% 205.9 / 1.1 74Q224H (CAG→CAT) gorglutathione reductase (NADPH)
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base T (30/39);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

GAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAG‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATAACACGGGT  >  NC_000913/3646727‑3647216
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        
gAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGcaca                                                                                                                                                                                                                                                  <  1:168662/251‑1 (MQ=255)
     aCCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AAc                                                                                                                                                                                                                                            <  1:316805/251‑1 (MQ=255)
       ccGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGAATAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGGAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGGTGAACGCCGAAAGCCCGCAT‑CTGCAAACC‑AACGc                                                                                                                                                                                                                                          >  1:541739/1‑251 (MQ=255)
           atatTCCGGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGGGATTAACGGGCTCGGAGCG‑CAAACGCATCTGGTTGTGAGTAAACATGCGCCCCTGCGCAGCTTCGACCCTGTGATTTCCGAAACGCCGGGCGGAGGGATGAAAACCGCAGGCCCGCCT‑CTGCACACC‑AACGccctc                                                                                                                                                                                                                                      >  1:516020/1‑248 (MQ=255)
                  gggggTGGAATACGTTATTGATTCTTGTGGCTTCCTCGCCCCTCCTGTTTGGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTGTGTGCGTAATCATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AAGGCCATCCCGAAAg                                                                                                                                                                                                                               <  2:568622/247‑1 (MQ=255)
                    gCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGGGTGGTTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGGAAACATGCGCCGCTGCGCAGCTTCGGCCCGGTGATTTCCGAAACGCTGGGCGGAGGGATGAAAGCCGAAGGCCCGCAT‑TCGCCCACC‑CACGCCCTCCCGAAAGCg                                                                                                                                                                                                                             >  1:47305/1‑251 (MQ=255)
                              cGGTATTGATGTTGAGGGCTTCTTCGCCCTGCCTGCTTTGGCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGaaaa                                                                                                                                                                                                                   <  2:200511/251‑1 (MQ=255)
                                gTATTGATTCTGATGGCTTCGTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAAt                                                                                                                                                                                                                 <  1:95909/251‑1 (MQ=255)
                                     gATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCg                                                                                                                                                                                                             <  1:63912/250‑1 (MQ=255)
                                      aTTCTGATGGCTTCGTCGCCCTTCCGGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTGGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCGCGAAAGCGGTAGTGAAAAATACCGAt                                                                                                                                                                                                           <  1:208228/251‑1 (MQ=255)
                                               gCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTg                                                                                                                                                                                                  <  1:603947/251‑1 (MQ=255)
                                                 ttcttcGCTCTTCCGGGTGTGCCAGAGTGGGTTGCGGTTGTTGTGGTGGGTTAAATCGCCGTTGAGCTGGCGGGCGTGATGAACGGACTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGTAGCTTCGACCCGATGATTTCCGAATCGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCTAAAGCGGTAGTGAAAAATTCCGATGGTAGCCtta                                                                                                                                                                                                 <  2:559540/250‑3 (MQ=255)
                                                  tcttcGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGAc                                                                                                                                                                                                <  1:93852/250‑1 (MQ=255)
                                                   tttcGCACTTCCTGCTGTGTCAGGGCGCGGTGCGGTTGGTGGCTCGGGTTTAATCGCCGTTGAGCTGGCGGGCTTGATTTACGGCCTCTGCGCG‑AAAACGTATCTGTTTGTGCGTAAACATGCGCCGCTGCGCTGCTTCGACCCGATGAGTTCCGAAACGCGGGGCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCTTGACg                                                                                                                                                                                               <  2:39139/249‑1 (MQ=255)
                                                     tcGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAG‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACg                                                                                                                                                                                               <  1:194991/248‑1 (MQ=255)
                                                              cTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCCCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACCGCCTCGGCGCGGAAAAC‑CCTCTGTTTGTGCGGAAACATGGGCCGCTTGGCCGCTTCCGCCCGCTGTTTTCCCGAACGCTGGTCGGGGTGGTGAGCCCCCAAAGCCCGCAT‑CTGC‑CACCCAACGCCATCCCCGAAGAGGGAGGGGAAAATACCGCTGGGAGCCTGAAggttgagctggg                                                                                                                                                                                    >  2:106103/1‑249 (MQ=255)
                                                                      ccAGAGCGCGTGGCGGTTGTTGGCGCTGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAG‑CTGCACACC‑AACGCCATCCCGaaa                                                                                                                                                                                                                                <  2:38889/198‑1 (MQ=255)
                                                                      ccAGAGCGCGTGGCGGTTGTTGGCGCTGGTTACATCGCCGTTGAGCTGGCGGGCGGGATTAACGGCCTCGGCGCG‑AAAAAGCATCTGTTTGTGCGTAAACATGCGCCGCCGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAG‑CTGCACACC‑AACGCCATCCCGaaa                                                                                                                                                                                                                                >  1:38889/1‑198 (MQ=255)
                                                                          agCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGGGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAAAGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCAAACC‑AACGCCATCCCGCAAGCGGGAGTGAAAAATAACGATGGTAGCCTGGCGCTGGAGCTGGGAGATGGGCGCAGt                                                                                                                                                                       >  1:412888/1‑251 (MQ=255)
                                                                               ggggCGGTTTTTGGAGCGGGTTACATCGCCGTTTAGCGGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAAc                                                                                                                                                                  <  2:459975/249‑1 (MQ=255)
                                                                                 ggCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACg                                                                                                                                                                 <  1:579082/250‑1 (MQ=255)
                                                                                     gttgttGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCt                                                                                                                                                                                                   >  1:46399/1‑212 (MQ=255)
                                                                                           ggCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGcc                                                                                                                                                       <  1:115376/250‑1 (MQ=255)
                                                                                               cgGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGAtt                                                                                                                                                  >  1:187854/1‑251 (MQ=255)
                                                                                                 ggTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTg                                                                                                                                                <  1:396018/251‑1 (MQ=255)
                                                                                                       aTCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGTCCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCTCACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGAt                                                                                                                                          <  1:469155/251‑1 (MQ=255)
                                                                                                            cGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTc                                                                                                                                     <  1:106103/251‑1 (MQ=255)
                                                                                                             gTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTcg                                                                                                                                    <  1:461794/251‑1 (MQ=255)
                                                                                                             gTTGAGCTGGCGGGCGTGATGAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGGGCGGAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTcg                                                                                                                                    <  2:367951/251‑1 (MQ=255)
                                                                                                                        gggcgTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGcc                                                                                                                          >  1:306670/1‑250 (MQ=255)
                                                                                                                         ggcgTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACg                                                                                                                                                                                                                                           >  1:226421/1‑136 (MQ=255)
                                                                                                                         ggcacgATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACg                                                                                                                                                                                                                                           <  2:226421/131‑1 (MQ=255)
                                                                                                                                 taacgaactaggCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGATGGTCGAAGTGATGAACGCAGAAGGCCCGCCT‑CTGCACCCC‑AACGCCATCCCGAAAGCGGGAGTGAACAATACCGATGGTCGCCTGGCGCTGGAGCTGGCAGATGGTCGCAGAGACACGGTGGATTGCCTGATTTGGGGGATTGGTAGAGGGCCTGCCaagtaaaac                                                                                                                 >  2:620154/11‑243 (MQ=255)
                                                                                                                                 tAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGCTTGCCTGATTTGGGCGATTGGTCGCGAGCCTGACAATGacaac                                                                                                                 >  2:333657/1‑250 (MQ=255)
                                                                                                                                           ggCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGTTTGGTCGCGAGCCTGCCAATGACAACATCAACCTgg                                                                                                       >  1:19747/1‑250 (MQ=255)
                                                                                                                                               cg‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGcc                                                                                                  >  1:468325/1‑251 (MQ=255)
                                                                                                                                                           tggTTTGTGCGCAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAATTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCGCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTa                                                                                        <  1:261574/248‑1 (MQ=255)
                                                                                                                                                              tttGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGCAACGCTGGTCGCAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAAc                                                                                    >  1:534953/1‑251 (MQ=255)
                                                                                                                                                               ttgggcttaaccgtgccccgcTGCGCACCTTCGACCCGATGTTTCCCGAAACGCTGGGTGAAGGGAGTAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGGGGTAGTGAAAAATCCCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAAc                                                                                    <  2:137432/233‑1 (MQ=255)
                                                                                                                                                                          cATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGTTGGTAGCCTGACgctggagctgga                                                                                                                                                                                    >  2:296521/1‑143 (MQ=255)
                                                                                                                                                                          cATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACgctggagctgga                                                                                                                                                                                    <  1:296521/143‑1 (MQ=255)
                                                                                                                                                                              cgccgaagcgcaacatCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCG‑AGTCTGCACACA‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGAAGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGGGGATTGCCTTATTTGGGCGATTGGTCGCGAGCCTGCTAATGACAACATCAACCTGGAATCCGATGTCGTTAAAAAGACCGAAAAAGGCta                                                                     >  2:205743/16‑250 (MQ=255)
                                                                                                                                                                               gccgctggtcagtttctgccatatgatTTCGGAAACGCTGGTAGAGTTGATGAACGCCTGAGGACCGCAT‑CTGCACACC‑AATGCCATGCCGAAAGCGGGAGTGAAAAATACCGATGGTTGCCTGACGCTGGAGCTGGAAGAGGGTCGCAGTGAAACGGTGGATTGCTGGATTTGGGCGATTGGTCGCGAGCCTGACAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCtata                                                                   <  2:313284/229‑1 (MQ=255)
                                                                                                                                                                               gccgcTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCtata                                                                   <  1:256861/251‑1 (MQ=255)
                                                                                                                                                                                 cgcTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATAtc                                                                 <  1:351832/251‑1 (MQ=255)
                                                                                                                                                                                  gcTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGCGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTatgtcggaagagcggcgggtagggg                                                                >  2:230960/1‑227 (MQ=255)
                                                                                                                                                                                  gcTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGtt                                                                                         <  1:230960/226‑1 (MQ=255)
                                                                                                                                                                                   cTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATAtcgt                                                               >  1:535841/1‑251 (MQ=255)
                                                                                                                                                                                    tGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGCTGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATAtcgtc                                                              >  2:410716/1‑251 (MQ=255)
                                                                                                                                                                                     gcgaagatacGACCCGATGAGTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATAtcgtcg                                                             <  2:188777/242‑1 (MQ=255)
                                                                                                                                                                                          gCTTCGACCCGATGTGTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATAtcgtcgtcg                                                          <  2:586067/249‑1 (MQ=255)
                                                                                                                                                                                            ttCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCAGCAT‑CTGCACACC‑AACGACATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGCGCTGGACGATGGTCGCAGTGAAACGGTGGATTGCATGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAACCTAACGAAAAAGGCTATAACGTCATCGATa                                                       >  2:568893/1‑250 (MQ=255)
                                                                                                                                                                                               gACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAat                                                    >  1:578343/1‑250 (MQ=255)
                                                                                                                                                                                               gACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAat                                                    >  1:212579/1‑250 (MQ=255)
                                                                                                                                                                                                aCCCGATGATTTCCGAAACGCTGGTCGCAGAGATGAACGCCGCAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAatat                                                  >  1:525782/1‑251 (MQ=255)
                                                                                                                                                                                                aCCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAatat                                                  <  1:532953/251‑1 (MQ=255)
                                                                                                                                                                                                 ccagatgatTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAACAAGCTATATCGTCGTCGATAAATATc                                                 >  1:337104/4‑251 (MQ=255)
                                                                                                                                                                                                       gatTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGAGGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAaca                                           <  1:60188/251‑1 (MQ=255)
                                                                                                                                                                                                        atTTCCGAAACGCTGGTAGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAacac                                          >  1:149559/1‑251 (MQ=255)
                                                                                                                                                                                                                  cgcgtggcGAAGTGATGAACGCCGAAGTCCCGGAT‑CGGCACCCC‑AACGCCATCCGGAAAGCGGTAGTGAAAAATACCGATGTTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGaa                                <  2:468325/244‑1 (MQ=255)
                                                                                                                                                                                                                     tGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGt                             >  1:614430/1‑251 (MQ=255)
                                                                                                                                                                                                                     tGGTCGAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGAAGGTCGCAGTGAAACGGTGGCTTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATCAATATCAAACAACCAATATTTAAGGt                             >  2:56082/1‑251 (MQ=255)
                                                                                                                                                                                                                      ggTCGAAGTGATGACCGACGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGGCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAgg                              <  1:3546/249‑1 (MQ=255)
                                                                                                                                                                                                                       gTCGAAGTGTTGAACGCCGTAGGCCCGCAG‑CTGCACACC‑AACGCCATCCCGCTAGCGGTAGTGAAACATACCGATGGAAGCCTGACGATGGAGCTGGAAGATGGTAGCAGTGAAACGGTGGCATGCCTGATTTGGGCGATTGGTCGTGTGCCTGCCAATACCAACATCAACATGTAAGCAGCTGGCGTTGAAACTACCATAAAAGACTATTACGTAGTCGATAATTTTCTAAACACCACTATTGAAGGttt                           >  2:569667/1‑249 (MQ=255)
                                                                                                                                                                                                                          gAAGTGATGAACGCCGAAGGCCCGCAT‑CTGCACACC‑ACCGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGACTATATCGTCGTCGATAAATATCACAACACCAATATTGAAGGTAttt                         >  2:288074/1‑250 (MQ=255)
                                                                                                                                                                                                                          aaattgataaACGCCGAAGGCACGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAACTCCGATGTTAGCCTGCCGCTGGAGCGGGAATATGGTCGCATTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTa                        <  2:363849/247‑1 (MQ=255)
                                                                                                                                                                                                                            aGTGATGAACGCGGAGGCCCCGCGT‑CTGCACACC‑AACGCCATCCCGAAAGTTGAATTGAAAAATACCGAGGGTAGCCGTACGCTGGAGCTGGAAGAGGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATTAACCTGGAAGCAGCTGGCGTTAAAACGAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTAcg                      <  2:575757/251‑1 (MQ=255)
                                                                                                                                                                                                                              tgatgaACGCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTAcgcg                    >  1:133440/1‑251 (MQ=255)
                                                                                                                                                                                                                                atgaACGCCTAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTAAAAAATACCGATGGTAGCCGGACGCTGGAGCGTGAAGTTGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGt                  <  2:477897/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aaCGCCGAAGGCCCGCAT‑ATGCACACC‑AACGCCATCCCGCAAGCGGTAGAGAAAACTACCGATGGTAGCCTGACGCTGGCGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGCTTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGACGCCGCTGGCGTTAACACTACCGAAAAAGGCTATACCGTCGTCGATAAATATAAAAACACCACTATTGAGGGTCTATACGCGGTggg               >  2:63236/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cagCCGAAGGCCCGAAT‑CTGCACACC‑AACGCCATCCCGAAGGCGGTAGTAAAAAATACCGATGGTAGCGTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTTATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAATGTATTTACGCGGTGGGc              <  2:614430/249‑1 (MQ=255)
                                                                                                                                                                                                                                      gCCGAAGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCg             >  1:315833/1‑250 (MQ=255)
                                                                                                                                                                                                                                           aGGCCCGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATAaca       <  1:328380/251‑1 (MQ=255)
                                                                                                                                                                                                                                                cGCAT‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATAACACGGGt  <  1:372857/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        
GAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCG‑AAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAG‑CTGCACACC‑AACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATAACACGGGT  >  NC_000913/3646727‑3647216

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