Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 383,475 A→C intergenic (+175/‑126) unknown → / → unknown Myo‑inositol 2‑dehydrogenase 2 (EC 1.1.1.18)/Inositol transport system sugar‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1383,4750AC97.7% 98.9 / NA 43intergenic (+175/‑126)unknown/unknownMyo‑inositol 2‑dehydrogenase 2 (EC 1.1.1.18)/Inositol transport system sugar‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (23/19);  minor base T (1/0);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

ATGCTATTGCGTTTTCATGTGAACAAAAGCGTTGGATAGACGTAGAAGAGGTTGCGTAACGAATTCGTATTGAGCATGAACATAGGAGTCTTTGGACTCCTATTTTTTTGCTAAAAATAAAATATTTAGTTTATTTTTGCAGTAGTTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTCTTTTAATGATATAGATACT‑AAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATGTGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAGGGAGTAGCATTATGAAAAAGACACTACTTGCGTTACTTTGTGGTGGAATCCTATTAACGTCTCAAGCTGTTATTGCTAAAAATGTTGTAATTGGCGCTCCGATGGTTGCCTTCTCCGATAAATGGCAAAC  >  1/383242‑383719
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     
aTGCTATTGCGTTTTCATGTGAACAAAAGCGTTGGATAGACGTAGAAGAGGTTGCGTAACGAATTCGTATTGAGCATGAACATAGGAGTCTTTGGACTCCTATTTTTTTGCTAAAAATAAAATATTTAGTTTATTTTTGCAGTAGTTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTTTTTTTATGATATATATACT‑CAACGGACGATttttttt                                                                                                                                                                                                                                      >  1:501722/1‑250 (MQ=255)
              tCATGTGAACAAAAGCGTTGGATAGACGTAGAAGAGGTTGCGTAACGAATTCGTATTGAGCATGAACATAGGAGTCTTTGGACTCCTATTTTTTTGCTAAAAATAAAATATTTAGTTTATTTTTGCAGTAGTTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTCTTTTAATGATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTgg                                                                                                                                                                                                                       <  1:32871/251‑1 (MQ=255)
                            gCGTTGGATAGACGTAGAAGAGGTTGCGTAACGAATTCGTATTGAGCATGAACATAGGAGTCTTTGGACTCCTATTTTTTTGCTAAAAATAAAATATTTAGTTTATTTTTGCAGTAGTTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTCTTTTAATGATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTg                                                                                                                                                                                                                     <  1:753846/239‑1 (MQ=255)
                                 ggATAGACGTAGAAGAGGTTGCGTAACGAATTCGTATTGAGCATGAACATAGGAGTCTTTGGACTCCTATTTTTTTGCTAAAAATAAAATATTTAGTTTATTTTTGCAGTAGTTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTCTTTTAATGATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTgg                                                                                                                                                                                                                       <  1:129333/232‑1 (MQ=255)
                                       aCGTAGAAGAGGTTGCGTAACGAATTCGTATTGAGCATGAACATAGGAGTCTTTGGACTCCTATTTTTTTGCTAAAAATAAAATATTTAGTTTATTTTTGCAGTAGTTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAAACTTCTTTTTATGATATAGATACT‑CCACGGACGATTTGTTT‑CAAATTGTATATTTGTGAAATTTATATTCTCATATGAttt                                                                                                                                                                                              >  1:792644/1‑250 (MQ=255)
                                          tAGAAGAGGTTGCGTAACGAATTCGTATTGAGCATGAACATAGGAGTCTTTGGACTCCTATTTTTTTGCTAAAAATAAAATATTTAGTTTATTTTTGCAGTAGTTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTCTTTTAATGATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGAttgttg                                                                                                                                                                                           <  1:798291/251‑1 (MQ=255)
                                                        tAACGAATTCGTATTGAGCATGAACATAGGAGTCTTTGGACTCCTATTTTTTTGCTAAAAATAAAATATTTAGTTTATTTTTGCAGTAGTTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTCTTTTAATGATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGc                                                                                                                                                                             <  1:867489/251‑1 (MQ=255)
                                                        tAACGAATTCGTATTGAGCATGAACATAGGAGTCTTTGGACTCCTATTTTTTTGCTAAAAATAAAATATTTAGTTTATTTTTGCAGTAGTTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTCTTTTAATGATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGc                                                                                                                                                                             <  1:846799/251‑1 (MQ=255)
                                                               ttCGTATTGAGCATGAACATAGGAGTCTTTGGACTCCTATTTTTTTGCTAAAAATAAAATATTTAGTTTATTTTTGCAGTAGTTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTTCAGGAAAACTTCTTTTTATGATATAGATACT‑CAAAGGACGATTTTTTT‑CAAATTGTATATTGGTGAAATTTATATTTGCCTATGATTTTTGTGGTGGGGGATAGAAGTAttt                                                                                                                                                                      >  1:812807/1‑251 (MQ=255)
                                                                       gAGCATGAACATAGGAGTCTTTGGACTCCTATTTTTTTGCTAAAAATAAAATATTTAGTTTATTTTTGCAGTAGTTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGGAATCTTCTTTTTATTATATAGATAAT‑CAACGGACGATTTTTTT‑TAAATTGTATATTGGGGAAATTTATATTTGCATATGGttgttg                                                                                                                                                                                           >  1:357227/1‑222 (MQ=255)
                                                                       gAGCATGAACATAGGAGTCTTTGGACTCCTATTTTTTTGCTAAAAATAAAATATTTAGTTTATTTTTGCAGTAGTTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTCTTTTTATGATATAGATACT‑CAACGGACGATTTTTTT‑TAAATTTTATATTTGTGGAATTTATTTTTGCATATTATTGTTGTGGTGGGGGATAGCAGTATTTTAAttttt                                                                                                                                                              >  1:373157/1‑251 (MQ=255)
                                                                             gAACATAGGAGTCTTTGGACTCCTATTTTTTTGCTAAAAATAAAATATTTAGTTTATTTTTGCAGTAGTTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTCTTTTTATGATATAGATACT‑CAAAGGACGATTTGTTT‑CAAATTGTATAtttg                                                                                                                                                                                                                       >  1:788322/1‑186 (MQ=255)
                                                                                                                     taaaataTTTAGTTTATTTTTGCAGTAGTTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTCTTTTAATGATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAGGGAGTAGCATTATGAAAAAg                                                                                                                <  1:78656/251‑1 (MQ=255)
                                                                                                                             ttAGTTTATTTTTGCAGTAGTTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTCTTTTAATGATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAGGGAGTAGCATTATGAAAAAGACactact                                                                                                        <  1:144002/251‑1 (MQ=255)
                                                                                                                                       tttGCAGTAGTTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTCTTTTAATGATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAGGGAGTAGCATTATGAAAAAGACACTACTTGCGTTACtt                                                                                              <  1:522467/251‑1 (MQ=255)
                                                                                                                                               agtTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTTCCGGAAATCTTTTTTTTATTATATAAATACT‑CAACGGACGATTTGTTT‑CAAATTGTTTTTTTGTGGAATTTATATTCTCATATGATTGTTGTGGGGGGGTATAGCAGTATTTTTATTTTTTTTTTTTTTCGTTTAACCGATAAAAAGGGGGAGCATTATGAAAAAAAAACTACTTTCGTTTCTTTTtggtgg                                                                                      >  1:540316/1‑251 (MQ=255)
                                                                                                                                               agtTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGTAGGAAATTTTCTTTTAATGATATAGATACT‑CAACGGACGATTTTTTTCCAAATT‑TATATTGGTGAAATTTATTTTCGCATATGATTGTTGTGGCGGGGTATAGCAGAATTTTAATTTTTGTTTTTTTTTGTTTAATCGATTAAAGGGAGGAGCATTATTAAAAAGAAAATACTTGCGTTTCTTTgtggggg                                                                                      >  1:9863/1‑248 (MQ=255)
                                                                                                                                                  tAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGGAATATTCTTTTTATTATATAGATAATCTAAGGGA‑GATTTTTTT‑CTAAATGTATATTTGGGAAAATTATATTTGCCTATGTTTGTTGTTGGGGGGGGTTGCAGTATTTTTATTTTTTTTGTTTTTc                                                                                                                                                    >  1:539900/1‑186 (MQ=255)
                                                                                                                                                                   aaaTAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTCTTTTAATGATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAAtt                                                                                                                                                                 <  1:674773/156‑1 (MQ=255)
                                                                                                                                                                   aaaTAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTCTTTTAATGATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAGGGAGTAGCATTATGAAAAAGACACTACTTGCGTTAc                                                                                                <  1:755433/221‑1 (MQ=255)
                                                                                                                                                                   aaaTAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAAACTTCTTTTAATTATATAGATACT‑CAACGGACGATTTGTTT‑Caaaa                                                                                                                                                                                                                                 >  1:639233/1‑91 (MQ=255)
                                                                                                                                                                    aaTAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTCTTTTAATGATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAGGGAGTAGCATTATGAAAAAGACACTACTTGCGTTACTTTGTGGTGGAATCCTATTAACGTCTCAAGc                                                                 <  1:222559/251‑1 (MQ=255)
                                                                                                                                                                         gATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTCTTTTAATGATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAg                                                                                                                                                                           <  1:719460/140‑1 (MQ=255)
                                                                                                                                                                            ctatctCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTCTTTTAATGATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAg                                                                                                                                                                           <  1:232114/137‑1 (MQ=255)
                                                                                                                                                                              atctCACACTTTCAAAAAATGGCGGTTTCAGGAAAACTTCTTTTTATGATATAGATACT‑CAACGGAAGGTTTGTTT‑TAAAATGTTTATTGGTGAAAATTATATTTGCATATGATTGTTGTGGTGGGGGATAGCAGTATTTTAATTTTTGTTGTTATTTGTTTTAACGATAAAAGGGAGTAGCATTATGAAAAAGACACTACTTGCGTTACTTTGTGGGGGGATCCTATTAAAGTCTCAAGcctggctcttt                                                       >  1:288336/1‑241 (MQ=255)
                                                                                                                                                                              atctCACACTTTCAAAAAATGGCGGTTGCCGGGAAACTTTTTTTTATGATATTGATAAT‑CACCGGGCCATTTGTTT‑TTAATTGTATATTTGTGAAATTTATATTCGCCTATGATTGTTGTGGGGGGGGATAGCAGAATTTTTATTTTTTTTGTTTTTCGTTTAAACGATAAAAGGGAGTAGCAATATGAAAAAGACACCACTTTGGTTACtt                                                                                              >  1:187939/1‑212 (MQ=255)
                                                                                                                                                                              atctCACACTTTCAAAAAATGGCGGTTGCAGGAAAACTTTTTTTTATTATATAGATAAT‑CAACGGACGATTTGTTT‑TAAATTTTATATTGGGGAAATTTTTATTCGCATATGATTGTTTTGGGGGGGGATAGGAGTATTTTAATTTTTGTTGTTATTTGTTTAAACGATAAAAAGGAGGAGCAATATGAAAAAAACACTACTTTCGTTACTTTTTGGGGGAAACCTATTAACGGCTCAAACTGGTAttt                                                         >  1:868719/1‑248 (MQ=255)
                                                                                                                                                                                                  ggCGGTTGCAGGAAATCTTCTTTTAATGATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAGGGAGTAGCATTATGAAAAAg                                                                                                                >  1:582233/1‑174 (MQ=255)
                                                                                                                                                                                                                           aTGATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAGGGAGTAGCATTATGAAAAAGACACTACTTGCGTTACTTTGTGGTGGAATCCTATTAACGTCTCAAGCTGTTATTGCTAAAAATGTTGTAATTGGCGCTCCGATGGTTGCCTTCTCCGATaaa          >  1:497826/1‑251 (MQ=255)
                                                                                                                                                                                                                             gATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTc                                                                                                                                                                     <  1:53562/94‑1 (MQ=255)
                                                                                                                                                                                                                             gATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAGGGAGTAGCATTATGAAAAAGACACTACTTGCGTTACTTTGTGGTGGAATcc                                                                                 >  1:19916/1‑178 (MQ=255)
                                                                                                                                                                                                                             gATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAGGGAGTAGCATTATGAAAAAGACACTACTTGCGTTACTTTGTGGTGGAATCCTATTAACGTCTCAAGCTGTTATTGCTAAAAAtgttgt                                            <  1:351754/215‑1 (MQ=255)
                                                                                                                                                                                                                             gATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAGGGAGTAGCATTATGAAAAAGACACTACTTGCGTTACTTTGTGGTGGAATCCTATTAACGTCTCAAGCTGTTATTGCTAAAAATGTTGTAATTGGCGCTCCGATGGTTGCCTTCTCCGATAAATg        >  1:657651/1‑251 (MQ=255)
                                                                                                                                                                                                                             gATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAGGGAGTAGCATTATGAAAAAGACACTACTTGCGTTACTTTGTGGTGGAATCCTATTAACGTCTCAAGCTGTTATTGCTAAAAATGTTGTAATTGGCGCTCCGATGGTTGCCTTCTCCGATAAATg        >  1:251399/1‑251 (MQ=255)
                                                                                                                                                                                                                             gATATAGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAGGGAGTAGCATTATGAAAAAGACACTACTTGCGTTACTTTGTGGTGGAATCCTATTAACGTCTCAAGCTGTTATTGCTAAAAATGTTGTAATTGGCGCTCCGATGGTTGCCTTCTCCGATAAATg        >  1:601993/1‑251 (MQ=255)
                                                                                                                                                                                                                               tataGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGAttgttg                                                                                                                                                                                           <  1:270431/70‑1 (MQ=255)
                                                                                                                                                                                                                               tataGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAggg                                                                                                                                  >  1:451216/1‑127 (MQ=255)
                                                                                                                                                                                                                               tataGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAGGGAGTAGCATTATGAAAAAGACACTACTTGCGTTACTTTGTGGTGGAATCCTATTAACGTCTCAAGCTGTTATTGCTAAAAATGTTGTAATTGGCGCTCCGATGGTTGCCTTCTCCGATAAATgg       >  1:265736/1‑250 (MQ=255)
                                                                                                                                                                                                                                 taGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGAtt                                                                                                                                                                                               <  1:321544/64‑1 (MQ=255)
                                                                                                                                                                                                                                 taGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAGGGAGTAGCATTATGAAAAAGACACTACTTGCGtt                                                                                                  >  1:407293/1‑157 (MQ=255)
                                                                                                                                                                                                                                 taGATACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAGGGAGTAGCATTATGAAAAAGACACTACTTGCGTTACTTTGTGGTGGAATCCTATTAACGTCTCAAGCTGTTATTGCTAAAAATGTTGTAATTGGCGCTCCGATGGTTGCCTTCTCCGATAAATGGCaa    >  1:54215/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aTACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAggg                                                                                                                                  <  1:126792/122‑1 (MQ=255)
                                                                                                                                                                                                                                    aTACT‑CAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATATGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAGGGAGTAGCATTATGAAAAAGACACTACTTGCGTTACTTTGTGGTGGAATCCTATTAACGTCTCAAGCTGTTATTGCTAAAAATGTTGTAATTGGCGCTCCGATGGTTGCCTTCTCCGATAAATGGCAAAc  >  1:293807/1‑250 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     
ATGCTATTGCGTTTTCATGTGAACAAAAGCGTTGGATAGACGTAGAAGAGGTTGCGTAACGAATTCGTATTGAGCATGAACATAGGAGTCTTTGGACTCCTATTTTTTTGCTAAAAATAAAATATTTAGTTTATTTTTGCAGTAGTTAAATTTTCCTTCCAAAAAATAAGATCTATCTCACACTTTCAAAAAATGGCGGTTGCAGGAAATCTTCTTTTAATGATATAGATACT‑AAACGGACGATTTGTTT‑CAAATTGTATATTGGTGAAATTTATATTCGCATGTGATTGTTGTGGTGGGGTATAGCAGTATTTCAATTTTTGTTGTTATTCGTTTAACCGATAAAAGGGAGTAGCATTATGAAAAAGACACTACTTGCGTTACTTTGTGGTGGAATCCTATTAACGTCTCAAGCTGTTATTGCTAAAAATGTTGTAATTGGCGCTCCGATGGTTGCCTTCTCCGATAAATGGCAAAC  >  1/383242‑383719

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