Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,647,428 C→T A377V (GCG→GTG)  gor → glutathione reductase (NADPH)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,647,4280CT93.9% 139.7 / NA 49A377V (GCG→GTG) gorglutathione reductase (NADPH)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/26);  minor base G (1/2);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTATTTACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAA  >  NC_000913/3647188‑3647651
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                               
gTATTTACTCGGTGGGCGATAACACGGGTGCATTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAAGAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCTCCGATTGGTACTGTTGGTTTAACGGAACAGCAGGCGCGCGATCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATAcc                                                                                                                                                                                                                       <  2:72554/251‑1 (MQ=255)
    ttACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATAccgccg                                                                                                                                                                                                                   >  1:405329/1‑251 (MQ=255)
      aCGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCAGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCTTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTc                                                                                                                                                                                                                 <  1:218111/251‑1 (MQ=255)
       cgcgGTGGGCGATAACACGGGTGCATGGGAGCCGACACCGTTGGCAGTTGCAGCGGGTCGCCGTCTCTCTTAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGCTTACTGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTca                                                                                                                                                                                                                <  2:120589/251‑1 (MQ=255)
                 gatacaacGGGTGCAGGGAGGCTGACACCGGTGGCAGTGCCCGGGGGCGGCCGTCTCTCTAAACCCTGGTTTAATAACAGGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCGTCGGCCACCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTAGGGCGACGATCAGGTGAAAGTGTATAAAGCCTCTTTCACCGTGAGTTATACCGCCGTCACCACTCACCg                                                                                                                                                                                                      <  2:42052/245‑1 (MQ=255)
                 gATAACACGGGTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCCGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTccccc                                                                                                                                                                                                      >  2:127610/1‑247 (MQ=255)
                          ggTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCGGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTg                                                                                                                                                                                             <  1:97562/251‑1 (MQ=255)
                               aGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGACATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCATCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGc                                                                                                                                                                                         >  1:266735/1‑250 (MQ=255)
                                  ggAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGCGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGGGCCGCATg                                                                                                                                                                                      >  2:150208/1‑250 (MQ=255)
                                   gAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGa                                                                                                                                                                                     >  1:322025/1‑250 (MQ=255)
                                       cggACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCt                                                                                                                                                                                 <  1:396350/249‑1 (MQ=255)
                                        tGACACCGGTGGCAGCTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGCATAAATCCTCTTTCACCGGGATGAATACCGCCGTCACCACTCAACGCCAGCCGTGCCGCATGAGGCTg                                                                                                                                                                                >  2:213211/1‑251 (MQ=255)
                                                  ggCAGTTGCAGCGGGTGTCCGTCTCTCTGAACGCCTGTTGAATAACAAGCCGGATGAGCATCTGGATTACAGCAACAGTCCGACCGTTGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCCCCGTTATTTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTg                                                                                                                                                                      <  1:308038/251‑1 (MQ=255)
                                                   gaaGTTGCAGCGGGTCGCCTGCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGACCATCTGGATTACAGTATCATCCCGACGGTGGTCTCCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGATCAGTATGGCGACGTTCAGGTTAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCTCCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTgg                                                                                                                                                                     <  1:197817/249‑1 (MQ=255)
                                                     aGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGAt                                                                                                                                                                   <  1:429464/251‑1 (MQ=255)
                                                          gggCGGGGCGCCGGCTCTCTGAACGCCTGTATAAAAAAAAGCCGGATGAGCATCTGGATTTCAGCAAAATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTAGGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaa                                                                                                                                                              <  2:224284/249‑1 (MQ=255)
                                                          cAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGga                                                                                                                                                               >  1:457638/1‑250 (MQ=255)
                                                             cGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaagag                                                                                                                                                           <  1:387638/251‑1 (MQ=255)
                                                                 tCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaagagaagg                                                                                                                                                       >  2:171265/1‑250 (MQ=255)
                                                                  cGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGAt                                                                                                                                                      <  1:129147/251‑1 (MQ=255)
                                                                                  acgctttttAATAAACAAGCCGGTGGAGCATTGGGATTACAGCAACTTTCCGCCCGTGGTCTCCAGCCTTCCGCCGATTGGTACTGTTGGGTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGGGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACg                                                                                                                                        <  2:143500/245‑1 (MQ=255)
                                                                                         ttAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTgg                                                                                                                                 <  1:130962/250‑1 (MQ=255)
                                                                                           aataaCAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCAGCCGCCGATTGGAACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCCGTATGGCGACGATCAGGTGAAAGTGAATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGATGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATACACGGCATTGGCtt                                                                                                                              >  2:425860/1‑251 (MQ=255)
                                                                                            aaaaaacggCCGGATGAGTATCGGGATTAAAGAACCATGCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATGCCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCt                                                                                                                               <  2:186216/242‑1 (MQ=255)
                                                                                                  aaGCCGGATGAGCATCTGGATTACAGCAACATCCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGGGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTAt                                                                                                                        <  1:200930/251‑1 (MQ=255)
                                                                                                    gCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCCCGGCATTGGCTTTGGTATgg                                                                                                                      >  2:423269/1‑251 (MQ=255)
                                                                                                      cGGATGAGCATCTGGATTCCAGCAACATTCCGCCCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGcc                                                                                                                                                                                           <  1:349643/180‑1 (MQ=255)
                                                                                                      cGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCAGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGcc                                                                                                                                                                                           >  2:349643/1‑180 (MQ=255)
                                                                                                         aTGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGa                                                                                                                  >  1:321488/1‑250 (MQ=255)
                                                                                                           ggccttcgggATTAAAGAAACATCCCGTCCGTTGTTTGCAGCCACCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAAt                                                                                                               <  2:321488/243‑1 (MQ=255)
                                                                                                                                aTTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGgc                                                                                          >  1:214531/1‑251 (MQ=255)
                                                                                                                                   ccGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTGTAACGGAGCCGCAGGCGCTCGACCAGTATGGCGACGATCAGGTGAAAGTGTATAAACCCTCTGTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCt                                                                                       <  1:456497/251‑1 (MQ=255)
                                                                                                                                            gTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTCTACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGGTgggg                                                                             >  2:442580/1‑251 (MQ=255)
                                                                                                                                                       ccccccgATGGGTACTGTGGGTTTACCGGAGCGGCAGGCGACGGAGCATTAGGGCGACGATCAGGTAAAAGGGAATAAAGCCTCTTGCCCCGTGATGTATACCGCCGTCACCACCCACCGCCAGCCGTGCCGCATGAAGCTGGTGGGCGTTGGATCGGAAGAGAAGATTGGCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCGGACGATGGGTGCAACCaaaa                                                                   <  2:226611/247‑1 (MQ=255)
                                                                                                                                                            cgATTGGTTCTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCCCCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTGCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGAc                                                              <  1:156075/251‑1 (MQ=255)
                                                                                                                                                               ttGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTc                                                           >  1:117846/1‑251 (MQ=255)
                                                                                                                                                                       gTTGGTTTAACGGAATCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAAt                                                     <  1:334124/249‑1 (MQ=255)
                                                                                                                                                                                    aaCCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCAc                                      >  1:320633/1‑251 (MQ=255)
                                                                                                                                                                                                    aGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGCCGAAATGTTGCAGGGCTTCGCGGGGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCACCGGCGGCaga                        >  2:282714/1‑250 (MQ=255)
                                                                                                                                                                                                     gaattATGGCGACGATCAGGTGAAAGTGTATAAACCCTCTTTCACCTTGATTTATCCCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGTCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCagaag                      <  1:63527/247‑1 (MQ=255)
                                                                                                                                                                                                          aTGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGtt                  <  1:150208/250‑1 (MQ=255)
                                                                                                                                                                                                                  gaTCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAAt         >  1:21275/1‑251 (MQ=255)
                                                                                                                                                                                                                   accgggTGAAAGTTTATAAATCCTTTTCCACCGTGATGTATCCCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATg        <  2:9179/247‑1 (MQ=255)
                                                                                                                                                                                                                      aGGTGAAAGTTTATAATTCCCTTTTCCCCGTGATGTATCCCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGt     <  1:264385/251‑1 (MQ=255)
                                                                                                                                                                                                                      aGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCaa                                                                     <  1:158622/187‑1 (MQ=255)
                                                                                                                                                                                                                      aGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCaa                                                                     >  2:158622/1‑187 (MQ=255)
                                                                                                                                                                                                                        gTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGTGGCAACCAAAAAAGACTTCGACAATACCGTCACCATTCACCCAACGACGGCAGAAGATGTCGAGACAATGCGtt    >  1:250033/1‑250 (MQ=255)
                                                                                                                                                                                                                         tGAAAGTGTATACATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGCCGCCATTCACCCAACGGCGGCAGAAGAGTTCGGGACAATGCggta   >  2:457461/1‑247 (MQ=255)
                                                                                                                                                                                                                          gAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTaa  <  1:9904/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                               
GTATTTACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAA  >  NC_000913/3647188‑3647651

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: