Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,365,856 G→T T32T (ACC→ACA)  unknown ← Purine nucleotide synthesis repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,365,8560GT98.1% 149.7 / ‑3.9 53T32T (ACC→ACA) unknownPurine nucleotide synthesis repressor
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (35/17);  total (35/18)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GCATTATTTTGCGTATTACATAAAAACAGAGAATAGCCTTGGCGATAGCAATGTTCTTCCACCGCATGAATGATTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTGGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAGATTGAATAGGTTTTGGTTAAGGAGAATCGCGGGAAATCCGATGAAAAATTACACTTTTTTACATTTGGCTGGTTTTTTCTGTGAATTTTTTGTTAG  >  1/1365650‑1366081
                                                                                                                                                                                                              |                                                                                                                                                                                                                                 
gCATTAATTTGCGTATTACATAAAAACAGAGAATAGCCTTGGCGATAGCAATGTTCTTCCACCGCATGAATGATTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATa                                                                                                                                                                                        >  1:663703/1‑250 (MQ=255)
     attTTGCGTATTACATAAAAACAGAGAATAGCCTTGGCGATAGCAATGTTCTTCCACCGCATGAATGATTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATAc                                                                                                                                                                                       >  1:133540/1‑246 (MQ=255)
          gCGTATTACATAAAAACAGAGAATAGCCTTGGCGATAGCAATGTTCTTCCACCGCATGAATGATTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTgtgg                                                                                                                                                                                  <  1:282704/246‑1 (MQ=255)
               ttACATAAAAACAGAGAATAGCCTTGGCGATAGCAATGTTCTTCCACCGCATGAATGATTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATa                                                                                                                                                                         >  1:307711/1‑250 (MQ=255)
                  cATAAAAACAGAGAATAGCCTTGGCGATAGCAATGTTCTTCCACCGCATGAATGATTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGc                                                                                                                                                                      >  1:33600/1‑250 (MQ=255)
                         aCAGAGAATAGCCTTGGCGATAGCAATGTTCTTCCACCGCATGAATGATTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCAtt                                                                                                                                                              >  1:881368/1‑251 (MQ=255)
                             agaATAGCCTTGGCGATAGCAATGTTCTTCCACCGCATGAATGATTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTg                                                                                                                                                          >  1:152185/1‑251 (MQ=255)
                              gaATAGCCTTGGCGATAGCAATGTTCTTCCACCGCATGAATGATTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCg                                                                                                                                                                                                                          <  1:777862/186‑1 (MQ=255)
                              gaATAGCCTTGGCGATAGCAATGTTCTTCCACCGCATGAATGATTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTg                                                                                                                                                          >  1:723779/1‑250 (MQ=255)
                              gaATAGCCTTGGCGATAGCAATGTTCTTCCACCGCATGAATGATTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTg                                                                                                                                                          >  1:757727/1‑250 (MQ=255)
                              gaATAGCCTTGGCGATAGCAATGTTCTTCCACCGCATGAATGATTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTg                                                                                                                                                          >  1:843978/1‑250 (MQ=255)
                                             tAGCAATGTTCTTCCACCGCATGAATGATTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGt                                                                                                                                          >  1:371527/1‑251 (MQ=255)
                                                  aTGTTCTTCCACCGCATGAATGATTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTggtggt                                                                                                                                                                              >  1:761359/1‑210 (MQ=255)
                                                    gTTCTTCCACCGCATGAATGATTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGtggtgg                                                                                                                                                                               >  1:336841/1‑207 (MQ=255)
                                                    gTTCTTCCACCGCATGAATGATTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATa                                                                                                                                   >  1:540346/1‑251 (MQ=255)
                                                     ttcttcCACCGCATGAATGATTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGcc                                                                                                                                                                     >  1:211112/1‑216 (MQ=255)
                                                                   gaatgaTTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTc                                                                                                                                                                             >  1:75041/1‑194 (MQ=255)
                                                                            gacagAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTGGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGc                                                                                                            <  1:681165/246‑1 (MQ=255)
                                                                            cGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCt                                                                                                           <  1:827757/251‑1 (MQ=255)
                                                                             ggCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCta                                                                                                          >  1:92401/1‑251 (MQ=255)
                                                                                  aaTAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATAc                                                                                                                                                                                       >  1:870458/1‑169 (MQ=255)
                                                                                   aTAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTgtgg                                                                                                                                                                                  <  1:123797/173‑1 (MQ=255)
                                                                                     aaGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAg                                                                                                  <  1:222507/251‑1 (MQ=255)
                                                                                        gggCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAGAt                                                                                                <  1:90957/250‑1 (MQ=255)
                                                                                             tCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATgg                                                                                                                                           >  1:82900/1‑202 (MQ=255)
                                                                                              cGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATAc                                                                                                                                                                                       >  1:253329/1‑157 (MQ=255)
                                                                                               gCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTtctc                                                                                                                           <  1:335074/216‑1 (MQ=255)
                                                                                               gCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAGATTGAATAgg                                                                                        >  1:9999/1‑251 (MQ=255)
                                                                                                  ggtggtGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAGATTGAATAGGttt                                                                                     >  1:681628/1‑251 (MQ=255)
                                                                                                   gtggtgACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAAttt                                                                                                                               >  1:235749/1‑208 (MQ=255)
                                                                                                          cGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAGATTGAATAGGTTTTGGt                                                                                 >  1:590844/1‑247 (MQ=255)
                                                                                                           gATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTgtggtg                                                                                                                                                                                <  1:49635/151‑1 (MQ=255)
                                                                                                             tCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAGATTGAATAGGTTTTGGTTAAGgag                                                                          <  1:437407/251‑1 (MQ=255)
                                                                                                                        gATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATAc                                                                                                                                                                                       >  1:589785/1‑131 (MQ=255)
                                                                                                                        gATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAg                                                                                                  >  1:132666/1‑216 (MQ=255)
                                                                                                                        gATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAGATTGAATAGGtt                                                                                      >  1:857253/1‑228 (MQ=255)
                                                                                                                            tGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGt                                                                                                                                          >  1:589920/1‑172 (MQ=255)
                                                                                                                                        cTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTtgttttgttt                                                                                                                                                                                                                               >  1:202374/1‑75 (MQ=255)
                                                                                                                                        cTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGAt                                                                                                                                                                                         <  1:153171/113‑1 (MQ=255)
                                                                                                                                        cTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAg                                                                                                  >  1:446674/1‑200 (MQ=255)
                                                                                                                                        cTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAGATTGAATAGGTTTTGGTTAAGGAGAATCGCGGGAAATCCGATGAAAAATTa                                               >  1:468893/1‑251 (MQ=255)
                                                                                                                                           ttAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAGATTGAATAGGTTTTGGTTAAGGAGAATCGCGGGAAATCCGATGAAAAATTacac                                            <  1:252099/251‑1 (MQ=255)
                                                                                                                                                       ccACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAatat                                                                                                                  <  1:727261/169‑1 (MQ=255)
                                                                                                                                                       ccACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAGATTGAATAGGTTTTGGTTAAGGAGAATCGCGGGAAATCCGATGAAAAATTACACTTTTTTACAttt                                >  1:591126/1‑251 (MQ=255)
                                                                                                                                                                 gggTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTgtgg                                                                                                                                                                                  <  1:486375/95‑1 (MQ=255)
                                                                                                                                                                 gggTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTgtgg                                                                                                                                                                                  <  1:501204/95‑1 (MQ=255)
                                                                                                                                                                  ggTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTgtgg                                                                                                                                                                                  <  1:7182/94‑1 (MQ=255)
                                                                                                                                                                        taatTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAGATTGAATAGGTTTTGGTTAAGGAGAATCGCGGGAAATCCGATGAAAAATTACACTTTTTTACATTTGGCTGGTTTTTTCtgtg               >  1:425061/1‑251 (MQ=255)
                                                                                                                                                                               gATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAGATTGAATAGGTTTTGGTTAAgg                                                                            >  1:675043/1‑183 (MQ=255)
                                                                                                                                                                               gATCTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAGATTGAATAGGTTTTGGTTAAGGAGAATCGCGGGAAATCCGATGAAAAATTACACTTTTTTACATTTGGCTGGTTTTTTCTGTGAAttttt        >  1:766381/1‑251 (MQ=255)
                                                                                                                                                                                  cTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAAttt                                                                                                                               <  1:371110/129‑1 (MQ=255)
                                                                                                                                                                                  cTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCtattat                                                                                                      >  1:486734/1‑154 (MQ=255)
                                                                                                                                                                                  cTTTTATTGCTTGCAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAGATTGAATAGGTTTTGGTTAAGGAGAATCGCGGGAAATCCGATGAAAAATTACACTTTTTTACATTTGGCTGGTTTTTTCTGTGAATTTTTTGt     >  1:480253/1‑251 (MQ=255)
                                                                                                                                                                                            ttgcAAAACCTGTTGTTTTGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAGATTGAATAGGTTTTGGTTAAGGAGAATCGCGGGAAATCCGATGAAAAATTACACTTTTTTACATTTGGCTGGTTTTTTCTGTGAATTTTTTGTTAg  <  1:719424/244‑1 (MQ=255)
                                                                                                                                                                                                              |                                                                                                                                                                                                                                 
GCATTATTTTGCGTATTACATAAAAACAGAGAATAGCCTTGGCGATAGCAATGTTCTTCCACCGCATGAATGATTTCGGCGAAATAAGGGGCTTCGCTGGTGGTGACGATCATGCCGATGGATTTGGTAGTGTTGACTTTTAAACTGCGCGCCACCGCACTGGGTGAATAATTCAGATCTTTTATTGCTTGCAAAACCTGTTGTTTGGTTTCAGCGGCAACAAAACGGGTTTTATTGATAACGTGGGATACTGTGGTGGTCGATACGCCTGCCATTTTTGCGACATCTTTGATGGTAGCCATAATTTTCTCCGTAAATATTTTGGCTATTATAAAGATTGAATAGGTTTTGGTTAAGGAGAATCGCGGGAAATCCGATGAAAAATTACACTTTTTTACATTTGGCTGGTTTTTTCTGTGAATTTTTTGTTAG  >  1/1365650‑1366081

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