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.6% 127.8 / NA 47A377V (GCG→GTG) gorglutathione reductase (NADPH)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/24);  minor base G (1/2);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GTATTTACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGC  >  NC_000913/3647188‑3647663
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           
gTATTTACGCGGTGGGCGATAACACGGCTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATAcc                                                                                                                                                                                                                                   <  2:346298/251‑1 (MQ=255)
  aTTTACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCTCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTCTAAATCCTCTTTCACCGTGATGTATAccgc                                                                                                                                                                                                                                 >  1:28123/1‑251 (MQ=255)
  aTTTACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATAccgc                                                                                                                                                                                                                                 >  1:301860/1‑251 (MQ=255)
       cgcgGTGGGCGATAACACGGGTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTATTCACACGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTACCGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGGGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTc                                                                                                                                                                                                                             >  2:212478/1‑250 (MQ=255)
        gcgGTGTGCGATAACCAGGGGGCAGTGGAGATGCCGCCGGTGGCATTGTCAGCGGTTTCCCGTCTCTTGGAACGCCTGTTAAATAACAAGCCGGATGAGCATCTGTATTACACCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGATCCGCAGTCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTcac                                                                                                                                                                                                                           <  2:40027/251‑1 (MQ=255)
                             gCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGACGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCg                                                                                                                                                                                                      >  1:130714/1‑251 (MQ=255)
                                            ccGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGtgtg                                                                                                                                                                                       >  2:26700/1‑251 (MQ=255)
                                                   gcagtggcggggggTCCCCGTCTCCCTGACGCCCTGTTAAATAACAAGCCGGTGGGGCTCTTGGTTTACAGCAACATCCCGACCGTGGTCTTAACCCATCCGCCGATTGGTACTGTTGGTGTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGCCGATCAGGTGAAAGTGTATAAATCCTCTTCCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTgg                                                                                                                                                                                 <  2:130714/240‑1 (MQ=255)
                                                        tgcggggggTCCCCGTCCCTCTGACCGCTTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGCCCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCg                                                                                                                                                                             <  2:234009/244‑1 (MQ=255)
                                                           gggggttCGCCGTCTCTTTGAACCCCTTTTAAAAAACAGGCCGGATGGCAATCTGGATTACAGAAACATCCCGACCGGGGTCTTCAGCCTCCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCCCCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaag                                                                                                                                                                         <  2:447136/245‑1 (MQ=255)
                                                           agcgtgtccctctctcccTGAACGCCTTTTTAATAACAAGCCGTATGAGAATCTGTATTACAGCAACATTCCGACCTTGGTTTTCAGCCACCCGCCGATTGTTACGGTTGGTTTAACGGAACCGCAGGCGCTCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTGCACCGTGATGTATACCGCCGTCACCACTCACCGCTAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaag                                                                                                                                                                         <  1:285889/239‑1 (MQ=255)
                                                                  cGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGATAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgtagaga                                                                                                                                                                      >  2:401275/1‑247 (MQ=255)
                                                                              gAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCAGTCACCGCTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGCAGAGACGATTGGAGGTATTc                                                                                                                                                       >  2:498963/1‑250 (MQ=255)
                                                                               aaCGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCAc                                                                                                                                                     >  1:292664/1‑251 (MQ=255)
                                                                                  gCCTGTTTAATAACAAGCCGGATGAGCAGCTGGATTAGAGGAACATTCTGACCTTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGTGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGGATACCGCCGTCGCCTCGCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGc                                                                                                                                                  <  1:72117/251‑1 (MQ=255)
                                                                                              aaaaaGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTgg                                                                                                                                       <  1:330008/248‑1 (MQ=255)
                                                                                                   gCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGCTGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGTAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATgg                                                                                                                                  >  1:65986/1‑251 (MQ=255)
                                                                                                                 gttggtttaaacaaccatccggacggTGGTCTTCACCCATCCGCCTAGTGGTACGGGTGGTTGAACGGAACCGCAGGCTCCCGAGCAGTATGGCGACGATCAGGGGAAAGGGTATAAATCCTCTTTCACCGGGATGTATCCCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGc                                                                                                                      <  2:298233/226‑1 (MQ=255)
                                                                                                                    ggATTACAGCAACATCCCGACCGGGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTg                                                                                                                       <  2:209068/247‑1 (MQ=255)
                                                                                                                    ggATTACACCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCTCCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAg                                                                                                                    <  1:157314/250‑1 (MQ=255)
                                                                                                                         aCAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTc                                                                                                              <  1:228504/251‑1 (MQ=255)
                                                                                                                              aaCATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCggg                                                                                                         >  1:138101/1‑250 (MQ=255)
                                                                                                                                      gACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTg                                                                                                  >  1:310905/1‑250 (MQ=255)
                                                                                                                                              gtggcagccagcCGCCGATTGGTACTGTTGGTTTAACGGAACCGGAGGCGCGCGAGCAGTATGACGACGATCAGGTGAAAGTGTATAAATCCTCTTGCACCGTGAGGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATgggg                                                                                         <  1:465679/247‑1 (MQ=255)
                                                                                                                                                 tCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATCCCGCCGTCACCCCTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTTAAGATGGGGGc                                                                                       <  1:211591/250‑1 (MQ=255)
                                                                                                                                                  cacccgcccgcctgttGGTACTGTTGGTTTAACGGACCCGCAGGCGCGCGATCAGTATGGCGACGATCAGGTGAAAGTGTAAAAATCCTCTTCCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCaa                                                                                     <  2:202542/237‑1 (MQ=255)
                                                                                                                                                  cAGCCAACCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCaa                                                                                     <  1:505587/251‑1 (MQ=255)
                                                                                                                                                     ccATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCa                                                                                  <  1:9511/251‑1 (MQ=255)
                                                                                                                                                        tCCGCCGATTGGTCCTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCaaa                                                                                <  1:208446/250‑1 (MQ=255)
                                                                                                                                                              gATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACt                                                                         >  1:389458/1‑251 (MQ=255)
                                                                                                                                                                   gTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCCGTATGGCGACGCTCAGGTGAAAGTGTCTAAATCCTCTTTCACCGGGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGAATTCACGTCCTTAGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCTCTGAAGATGTGGGCaccccacaaggaattcgac                                                                    >  2:129131/1‑233 (MQ=255)
                                                                                                                                                                        gTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATa                                                                <  1:219102/250‑1 (MQ=255)
                                                                                                                                                                        gTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTAGGGCGACGATCAGGTGAAATTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCCCTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATa                                                                <  1:143856/250‑1 (MQ=255)
                                                                                                                                                                                       ccgaaGGCGCGAGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACTGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGTCGCTGAAGATGGGGGCACCCAAAAAACTCTTCGACAATACCGTCGCCATTCAccc                                                <  1:373518/247‑1 (MQ=255)
                                                                                                                                                                                            ggCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCATCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAAcgg                                           <  1:449139/251‑1 (MQ=255)
                                                                                                                                                                                                    aGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCCCCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTCGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCaga                                    >  2:447275/1‑250 (MQ=255)
                                                                                                                                                                                                      cAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGGGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCaga                                    >  2:439011/1‑248 (MQ=255)
                                                                                                                                                                                                      cAGTATGGCGAAGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGCTGTATACAGCCGTCACAACTCACCGCAAGCCGTGCCGCATGAAGCTGGTATGCGTTGGATCGGAAGAGACGATTGTCGGCATTCACGGCATTGGCGTTGGTATGGTCGCAATGTTGAAGGGGTTCGCGGTGGCGCTGAAGATGGGGGCAACCCAACAAGACTTCGACCATACCGTCGCCATTCACACAAcggagacaac                                    >  2:485980/1‑241 (MQ=255)
                                                                                                                                                                                                               gacgaTCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGAc                        <  1:44237/251‑1 (MQ=255)
                                                                                                                                                                                                                   aTCAGGTGAAAGTGTATAAACCCTCTTTCACCGGGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCATCGGCGGCAGAAGAGTTCGTGACAAtt                    <  1:157016/251‑2 (MQ=255)
                                                                                                                                                                                                                      aGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGt                 >  1:503529/1‑251 (MQ=255)
                                                                                                                                                                                                                          gAAAGTGTTTAAATCCTCTTTCGCCGTGATGTATACCGCCGTCACCCCTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTaa              <  1:403316/250‑1 (MQ=255)
                                                                                                                                                                                                                           aaaGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGGCGCCATTCCCCCAACGGCGGCAGAAGAGTTCGTGACAAGGCGTTaat             >  2:270075/1‑249 (MQ=255)
                                                                                                                                                                                                                             aGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTCTGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGGTAAATGt          >  2:458917/1‑251 (MQ=255)
                                                                                                                                                                                                                               tgtATAAATCCTCTTTCACCGTGATGTATCCCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGAc                                                                    <  1:109431/191‑1 (MQ=255)
                                                                                                                                                                                                                               tgtATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTAACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGAc                                                                    >  2:109431/1‑191 (MQ=255)
                                                                                                                                                                                                                                     aaTCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATTTTAAAGGGc  >  2:431827/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           
GTATTTACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGC  >  NC_000913/3647188‑3647663

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