Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 809,527 T→C R17G (AGA→GGA)  unknown ← NADH pyrophosphatase (EC 3.6.1.22), decaps 5'‑NAD modified RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1809,5270TC100.0% 175.7 / NA 59R17G (AGA→GGA) unknownNADH pyrophosphatase (EC 3.6.1.22), decaps 5'‑NAD modified RNA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/30);  total (29/30)

GAATTGATGGGTCTGATAGAAATGATTTAACTCCACGCCGCAGCTCAACAGATTAAATTGCGCTTGCGGCAGCGGCAGTTGATCGCGTAAAGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCTCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGGCCTTGTTTTTCGATAAAATTTTAGTCTTTCACAAACTCGCTGAAGGTGAGTTTGTACTTATCGCCGTCGCGATTGTA  >  1/809282‑809770
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
gAATTGATGGGTCTGATAGAAATGATTTAACTCCACGCCGCAGCTCAACAGATTAAATTGCGCTTGCGGCAGCGGCAGTTGATCGCGTAAAGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAg                                                                                                                                                                                                                                                >  1:891312/1‑251 (MQ=255)
gAATTGATGGGTCTGATAGAAATGATTTAACTCCACGCCGCAGCTCAACAGATTAAATTGCGCTTGCGGCAGCGGCAGTTGATCGCGTAAAGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAg                                                                                                                                                                                                                                                >  1:225328/1‑251 (MQ=255)
gAATTGATGGGTCTGATAGAAATGATTTAACTCCACGCCGCAGCTCAACAGATTAAATTGCGCTTGCGGCAGCGGCAGTTGATCGCGTAAAGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAg                                                                                                                                                                                                                                                >  1:250646/1‑251 (MQ=255)
gAATTGATGGGTCTGATAGAAATGATTTAACTCCACGCCGCAGCTCAACAGATTAAATTGCGCTTGCGGCAGCGGCAGTTGATCGCGTAAAGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAACATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAg                                                                                                                                                                                                                                                >  1:427538/1‑251 (MQ=255)
          gTCTGATAGAAATGATTTAACTCCACGCCGCAGCTCAACAGATTAAATTGCGCTTGCGGCAGCGGCAGTTGATCGCGTAAAGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCa                                                                                                                                                                                                                                      >  1:801053/1‑251 (MQ=255)
          gTCTGATAGAAATGATTTAACTCCACGCCGCAGCTCAACAGATTAAATTGCGCTTGCGGCAGCGGCAGTTGATCGCGTAAAGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCa                                                                                                                                                                                                                                      >  1:634957/1‑251 (MQ=255)
              gATAGAAATGATTTAACTCCACGCCGCAGCTCAACAGATTAAATTGCGCTTGCGGCAGCGGCAGTTGATCGCGTAAAGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAg                                                                                                                                                                                                                                  >  1:105496/1‑251 (MQ=255)
               aTAGAAATGATTTAACTCCACGCCGCAGCTCAACAGATTAAATTGCGCTTGCGGCAGCGGCAGTTGATCGCGTAAAGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGc                                                                                                                                                                                                                                 <  1:695971/251‑1 (MQ=255)
                  gAAATGATTTAACTCCACGCCGCAGCTCAACAGATTAAATTGCGCTTGCGGCAGCGGCAGTTGATCGCGTAAAGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGct                                                                                                                                                                                                                                >  1:569323/1‑249 (MQ=255)
                        aTTTAACTCCACGCCGCAGCTCAACAGATTAAATTGCGCTTGCGGCAGCGGCAGTTGATCGCGTAAAGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCtt                                                                                                                                                                                                                        >  1:395801/1‑251 (MQ=255)
                             aCTCCACGCCGCAGCTCAACAGATTAAATTGCGCTTGCGGCAGCGGCAGTTGATCGCGTAAAGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGc                                                                                                                                                                                                                                 <  1:892492/237‑1 (MQ=255)
                                     ccgcAGCTCAACAGATTAAATTGCGCTTGCGGCAGCGGCAGTTGATCGCGTAAAGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGCCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAg                                                                                                                                                                                                            >  1:48351/1‑250 (MQ=255)
                                                  gATTAAATTGCGCTTGCGGCAGCGGCAGTTGATCGCGTAAAGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCAttt                                                                                                                                                                                                   >  1:52469/1‑246 (MQ=255)
                                                               ttGCGGCAGCGGCAGTTGATCGCGTAAAGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTg                                                                                                                                                                                 <  1:254348/251‑1 (MQ=255)
                                                               ttGCGGCAGCGGCAGTTGATCGCGTAAAGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTg                                                                                                                                                                                 <  1:770861/251‑1 (MQ=255)
                                                                                         aaGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAGAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAAc                                                                                                                                                                                   <  1:63450/223‑1 (MQ=255)
                                                                                         aaGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTc                                                                                                                                                       <  1:419185/251‑1 (MQ=255)
                                                                                           gAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCa                                                                                                                                                      >  1:678646/1‑250 (MQ=255)
                                                                                                 tAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAg                                                                                                                                                                                                                                                >  1:751202/1‑154 (MQ=255)
                                                                                                         gTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGc                                                                                                                                                                                                                   <  1:830184/175‑1 (MQ=255)
                                                                                                         gTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGc                                                                                                                                                                                                                   <  1:836496/175‑1 (MQ=255)
                                                                                                                                  aaaTATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTg                                                                                                                                                                                 >  1:860595/1‑184 (MQ=255)
                                                                                                                                  aaaTATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTg                                                                                                                                                                                 >  1:865750/1‑184 (MQ=255)
                                                                                                                                    atatAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGc                                                                                                                                                                                                                   <  1:896260/148‑1 (MQ=255)
                                                                                                                                    atatAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAAc                                                                                                                                                                                   <  1:366786/180‑1 (MQ=255)
                                                                                                                                    atatAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAAc                                                                                                                                                                                   <  1:323095/180‑1 (MQ=255)
                                                                                                                                    atatAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTg                                                                                                                                                                                 <  1:655783/182‑1 (MQ=255)
                                                                                                                                    atatAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCa                                                                                                            >  1:109822/1‑251 (MQ=255)
                                                                                                                                    atatAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCa                                                                                                            >  1:897758/1‑251 (MQ=255)
                                                                                                                                    atatAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCa                                                                                                            >  1:734598/1‑251 (MQ=255)
                                                                                                                                              tGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGt                                                                                                   >  1:598094/1‑250 (MQ=255)
                                                                                                                                                         ttCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGt                                                                                                   <  1:739534/239‑1 (MQ=255)
                                                                                                                                                             ccGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATc                                                                                                                                                                                                                          <  1:62179/116‑1 (MQ=255)
                                                                                                                                                             ccGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTg                                                                                                                                                                                 <  1:454908/157‑1 (MQ=255)
                                                                                                                                                             ccGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAAtttttt                                                                                          <  1:390576/244‑1 (MQ=255)
                                                                                                                                                                    gcATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAAc                                                                                                                                                                                   <  1:343585/148‑1 (MQ=255)
                                                                                                                                                                    gcATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTaaaaa                                                                                                                                         <  1:606197/190‑1 (MQ=255)
                                                                                                                                                                       tAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAAc                                                                                                                                                                                   <  1:763803/145‑1 (MQ=255)
                                                                                                                                                                         gCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGTCCTTGttt                                                                       >  1:823719/1‑251 (MQ=255)
                                                                                                                                                                          cATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGTCCTTGtttt                                                                      >  1:301616/1‑251 (MQ=255)
                                                                                                                                                                          cATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGTCCTTGttt                                                                       <  1:724585/250‑1 (MQ=255)
                                                                                                                                                                                     acgaAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGTCCTTGTTTTTCGATAAAAtt                                                           >  1:207217/1‑251 (MQ=255)
                                                                                                                                                                                     acgaAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGTCCTTGTTTTTCGATAAAAtt                                                           >  1:785847/1‑251 (MQ=255)
                                                                                                                                                                                          aCCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGaa                                                                                                                                                                         <  1:301952/136‑1 (MQ=255)
                                                                                                                                                                                          aCCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTc                                                                                                                                                       <  1:158822/154‑1 (MQ=255)
                                                                                                                                                                                           cccAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGTCCTTGTTTTTCGATAAAATTTTAGt                                                      >  1:2017/1‑250 (MQ=255)
                                                                                                                                                                                           cccAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGTCCTTGTTTTTCGATAAAATTTTAGTc                                                     <  1:615579/251‑1 (MQ=255)
                                                                                                                                                                                                     tGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGTCCTTGTTTTTCGATAAAATTTTAGTCTTTCACAAAc                                           <  1:282300/251‑1 (MQ=255)
                                                                                                                                                                                                            gCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGTCCTTGTTTTTCGATAAAATTTTAGTCTTTCACAAACTCGCTGa                                    >  1:102402/1‑251 (MQ=255)
                                                                                                                                                                                                                aaaGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGTCCTTGTTTTTCGATAAAATTTTAGTCTTTCACAAACTCGCTGAAgg                                 <  1:247296/250‑1 (MQ=255)
                                                                                                                                                                                                                aaaGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGTCCTTGTTTTTCGATAAAATTTTAGTCTTTCACAAACTCGCTGAAgg                                 <  1:139470/250‑1 (MQ=255)
                                                                                                                                                                                                                         ttGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGTCCTTGTTTTTCGATAAAATTTTAGTCTTTCACAAACTCGCTGAAGGTGAGTTTGTa                       <  1:452961/251‑1 (MQ=255)
                                                                                                                                                                                                                          tGCCGTCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGTCCTTGTTTTTCGATAAAATTTTAGTCTTTCACAAACTCGCTGAAGGTGAGTTTGTAc                      <  1:515910/251‑1 (MQ=255)
                                                                                                                                                                                                                               tCAATTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTc                                                                                                                                                                                                                       <  1:679086/53‑1 (MQ=255)
                                                                                                                                                                                                                                  aTTAAATGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGTCCTTGTTTTTCGATAAAATTTTAGTCTTTCACAAACTCGCTGAAGGTGAGTTTGTACTTATCGc               >  1:95106/1‑250 (MQ=255)
                                                                                                                                                                                                                                       aTGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGTCCTTGTTTTTCGATAAAATTTTAGTCTTTCACAAACTCGCTGAAGGTGAGTTTGTAc                      >  1:545166/1‑238 (MQ=255)
                                                                                                                                                                                                                                        tGTATTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGTCCTTGTTTTTCGATAAAATTTTAGTCTTTCACAAACTCGCTGAAGGTGAGTTTGTACTTATCGCCGTcgcg        <  1:711970/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aTTTTTGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTc                                                                                                             >  1:496558/1‑147 (MQ=255)
                                                                                                                                                                                                                                              tttGATCCCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGTCCTTGTTTTTCGATAAAATTTTAGTCTTTCACAAACTCGCTGAAGGTGAGTTTGTACTTATCGCCGTCGCGATTGTa  >  1:554376/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
GAATTGATGGGTCTGATAGAAATGATTTAACTCCACGCCGCAGCTCAACAGATTAAATTGCGCTTGCGGCAGCGGCAGTTGATCGCGTAAAGAAAAATAGGTGCGGTTGTCGTTCGGTTGGGTTTCCACCAAATATAACGGTTGCTCCAGCCATTCGCCGATGCGCATAGCATGAAGCCCGACGAAACCCAGTTCGGTGGCATTGCCGAAAGGAAGCTTGCCGTCAATTAAATGTATTTTTGATCTCTGAGTGAAAAGCCAAAAGCTCTTATCTTCAGGCTGAATAGCGTTCATTTTTAGCACCTGTTTAACTGCGTGAGAAATACATAAATTTTTATTCACCCGATTTAAAAAAACGTTGTTTTTTGTGATCGCACTTTTCAATAAAAAGTCAATTTTTTTATTCAGGCGGCCTTGTTTTTCGATAAAATTTTAGTCTTTCACAAACTCGCTGAAGGTGAGTTTGTACTTATCGCCGTCGCGATTGTA  >  1/809282‑809770

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