Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 2,688,963 A→C K404T (AAG→ACG)  purB → adenylosuccinate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,688,9630AC96.2% 1393.1 / 26.2 371K404T (AAG→ACG) purBadenylosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base A (5/9);  new base C (178/177);  total (183/188)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

Only 100 of 371 total aligned reads displayed.

GCATATCAATCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGAAGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCACGATGGTTGA  >  DH1/2688818‑2689110
                                                                                                                                                 |                                                                                                                                                   
gCATATCAATCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGa                                                                                                                                                <  7:289717/151‑1 (MQ=255)
    atCAATCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGaa                                                                                                                                               >  5:231214/1‑148 (MQ=255)
    atCAATCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAaga                                                                                                                                            <  3:989425/151‑1 (MQ=255)
     tCAATCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGaa                                                                                                                                               >  1:973895/1‑147 (MQ=255)
      cAATCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAagag                                                                                                                                           <  4:677135/150‑1 (MQ=255)
       aaTCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCt                                                                                                                                         >  3:1424124/1‑151 (MQ=255)
          ccaccaTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATTCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGAc                                                                                                                                      <  1:1368224/151‑1 (MQ=255)
          ccACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGAAGCTGAAAGAGCTGa                                                                                                                                       <  2:757408/150‑1 (MQ=255)
          ccACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCAGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGAc                                                                                                                                      <  8:376336/151‑1 (MQ=255)
             cccTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGAAGCTGAAAGAGCTGACTcg                                                                                                                                   <  3:751997/151‑1 (MQ=255)
                tGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCgg                                                                                                                                >  5:1347715/1‑151 (MQ=255)
                      gCGTGAGCAAAATGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAACACCGTACGAGACGCTGAAAGAGCGGACTCGCGGTAAgcg                                                                                                                          >  6:1074953/1‑151 (MQ=255)
                        gTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGcgcg                                                                                                                        >  7:1127283/1‑151 (MQ=255)
                          gAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGtt                                                                                                                      <  5:808750/151‑1 (MQ=255)
                          gAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGtt                                                                                                                      <  5:662547/151‑1 (MQ=255)
                          gAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGtt                                                                                                                      <  2:1230041/151‑1 (MQ=255)
                           aGCAAACTGGACGTGAACCGTGACCAGCTGCTGGATGAACTGGAGCACAAGTGGGACGTGCTGGCTGAACCAATCCATACAGTTATGCGTCGCTATGTCATCGAAAAACCGTACGAGACTCTGAAAGAGCTGACTCGCGGTAAGCGCGTTg                                                                                                                     <  7:77158/151‑1 (MQ=255)
                                 cTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCg                                                                                                               >  8:1199617/1‑151 (MQ=255)
                                  tGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCg                                                                                                               >  7:390855/1‑150 (MQ=255)
                                   ggAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGaa                                                                                                             <  3:648675/151‑1 (MQ=255)
                                   ggAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGa                                                                                                              <  2:240719/150‑1 (MQ=255)
                                   ggAAGTGAACCGGGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGaa                                                                                                             <  6:398378/151‑1 (MQ=255)
                                    gAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGa                                                                                                              >  6:179677/1‑149 (MQ=255)
                                      aGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGc                                                                                                          <  5:587977/151‑1 (MQ=255)
                                         gAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATg                                                                                                       >  4:1197038/1‑151 (MQ=255)
                                           aCCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGaa                                                                                                     <  1:317044/151‑1 (MQ=255)
                                           aCCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGa                                                                                                      >  3:267430/1‑150 (MQ=255)
                                            ccGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACTGTTATGCGTCGCTATGGCATCGAAAAACCGTACGGGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAg                                                                                                    <  6:836150/151‑1 (MQ=255)
                                              gTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCa                                                                                                  <  1:658180/151‑1 (MQ=255)
                                                gACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGt                                                                                                <  7:142333/151‑1 (MQ=255)
                                                gACAATCTGATGGATGAACTGGATCACAACAGGGAACTGCTCGCTGAACCAATCAAGACAGTTATTCGTCGCTAAGGCATCGAAAAACCGTACGAGACGAAAAAAGAGATCACTCGCGCTAAGAGCGTAGACGAAGAATGCATGAAGc                                                                                                   >  2:1066436/1‑148 (MQ=255)
                                                 aCCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGtt                                                                                               >  4:189047/1‑151 (MQ=255)
                                                   cATCTGCTGGATGAACTGGATGACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGGCGCTATGGCATCGAAAAAGCGTACGAGAGGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTa                                                                                             <  8:1314172/151‑1 (MQ=255)
                                                   cATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGttt                                                                                              >  1:1424860/1‑150 (MQ=255)
                                                   cATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTa                                                                                             <  7:495352/151‑1 (MQ=255)
                                                       tgctgGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGATGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCg                                                                                          >  4:932679/1‑150 (MQ=255)
                                                        gctgGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAACCCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTAt                                                                                            <  3:860166/147‑1 (MQ=255)
                                                        gctgGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGa                                                                                         <  6:335965/150‑1 (MQ=255)
                                                        gctgGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGAt                                                                                        >  2:879439/1‑151 (MQ=255)
                                                         ctgGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATg                                                                                       >  8:1534875/1‑151 (MQ=255)
                                                          tgGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATACGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGACGCAGTTTATCGATgg                                                                                      >  6:695242/1‑151 (MQ=255)
                                                           ggATGAACTGGATCACAACTGGGAAGTGCTTGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATgg                                                                                      <  8:746070/150‑1 (MQ=255)
                                                           ggATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATg                                                                                       <  7:181058/149‑1 (MQ=255)
                                                           ggATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATg                                                                                       <  6:983485/149‑1 (MQ=255)
                                                           ggATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGt                                                                                     <  6:1402067/151‑1 (MQ=255)
                                                               gAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTgg                                                                                 >  8:62709/1‑151 (MQ=255)
                                                               gAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTgg                                                                                 <  2:294790/151‑1 (MQ=255)
                                                                 aCTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCg                                                                               >  5:383822/1‑151 (MQ=255)
                                                                  cTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGt                                                                              >  2:341854/1‑151 (MQ=255)
                                                                   tGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGtt                                                                             >  5:470503/1‑151 (MQ=255)
                                                                    ggATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTg                                                                            >  2:959221/1‑151 (MQ=255)
                                                                    ggATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTg                                                                            <  3:1082261/151‑1 (MQ=255)
                                                                        cacaACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCa                                                                         <  8:1473813/150‑1 (MQ=255)
                                                                            aCTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCagaaga                                                                    <  2:1516816/151‑1 (MQ=255)
                                                                             cTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCagaaga                                                                    >  1:1052367/1‑150 (MQ=255)
                                                                             cTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCagaaga                                                                    >  4:1512598/1‑150 (MQ=255)
                                                                               gggAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACCCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCAAGAAGCAGTTTATCGATGGTCTGGCGTTGCCagaaga                                                                    >  6:515424/1‑148 (MQ=255)
                                                                                ggAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGAAGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAagag                                                                <  1:154075/151‑1 (MQ=255)
                                                                                   aGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGaaa                                                             <  7:839591/151‑1 (MQ=255)
                                                                                   aGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGaaa                                                             <  7:1258420/151‑1 (MQ=255)
                                                                                   aGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGAAGCTGAAATAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGaaa                                                             <  4:1333973/151‑1 (MQ=255)
                                                                                        tggctgAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGCAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCg                                                        >  1:236321/1‑151 (MQ=255)
                                                                                        tggctgAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCg                                                        >  4:54239/1‑151 (MQ=255)
                                                                                        tggctgAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCCTTGCCAGAAGAAGAGAAAGCCCg                                                        >  8:269286/1‑151 (MQ=255)
                                                                                        tggctgAACCAATCCAGACAGTTAAGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCg                                                        >  8:1402043/1‑151 (MQ=255)
                                                                                         ggctgAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCg                                                        <  6:896638/150‑1 (MQ=255)
                                                                                            tgAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTg                                                    >  2:1396430/1‑151 (MQ=255)
                                                                                             gAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGa                                                   <  3:419808/151‑1 (MQ=255)
                                                                                              aaCCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGaa                                                  <  1:696185/151‑1 (MQ=255)
                                                                                              aaCCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGaa                                                  <  7:120439/151‑1 (MQ=255)
                                                                                               aCCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCCAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGACAAAGCCCGCCTGaa                                                  >  1:53518/1‑150 (MQ=255)
                                                                                                 cAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGc                                               <  6:60783/151‑1 (MQ=255)
                                                                                                 cAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGc                                               <  3:580191/151‑1 (MQ=255)
                                                                                                 cAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGc                                               <  1:1022500/151‑1 (MQ=255)
                                                                                                     ccAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATg                                           <  4:82166/151‑1 (MQ=255)
                                                                                                           aGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCg                                     <  2:636181/151‑1 (MQ=255)
                                                                                                            gTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCgg                                    >  7:576125/1‑151 (MQ=255)
                                                                                                             ttATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGc                                   >  5:1101521/1‑151 (MQ=255)
                                                                                                               aTGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTa                                 <  4:460536/151‑1 (MQ=255)
                                                                                                                 gcgtttatATGGCATCGAAAAACCGTAGGAGAGGCTGAAAGAGCTTACTCTCGTTAATCTCGTTGACTAAGAATTAATGAATCATTTTATGGATGTTCTGGCGTTTACATAAGAAGAGAAAGCCCTCATGAAAGCGATGACGCCGGCTaa                                <  5:1228062/143‑1 (MQ=255)
                                                                                                                 gCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAAc                               >  5:1549591/1‑151 (MQ=255)
                                                                                                                   gTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACta                             <  7:885442/151‑1 (MQ=255)
                                                                                                                   gTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACt                              >  2:459775/1‑150 (MQ=255)
                                                                                                                     ggCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACtata                           <  8:573302/150‑1 (MQ=255)
                                                                                                                        tATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTg                        >  5:400888/1‑151 (MQ=255)
                                                                                                                         aTGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTgg                       >  8:20370/1‑151 (MQ=255)
                                                                                                                           ggCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTc                     <  1:1230837/151‑1 (MQ=255)
                                                                                                                              aTCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCg                    >  8:440889/1‑149 (MQ=255)
                                                                                                                                   aaaaCCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATc             >  1:56311/1‑151 (MQ=255)
                                                                                                                                    aaaCCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTACCTATATTGGTCGAGCTagcc            >  4:436808/1‑148 (MQ=255)
                                                                                                                                     aaCCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCAc           >  3:229652/1‑151 (MQ=255)
                                                                                                                                      aCCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCACg          <  4:389561/151‑1 (MQ=255)
                                                                                                                                      aCCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCACg          >  6:1049722/1‑151 (MQ=255)
                                                                                                                                         gTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCACGATg       <  5:548930/151‑1 (MQ=255)
                                                                                                                                           aCGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCACGATGGt     <  5:333284/151‑1 (MQ=255)
                                                                                                                                           aCGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCACGATGGt     <  7:922709/151‑1 (MQ=255)
                                                                                                                                           aCGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCACGATGGt     >  3:1359905/1‑151 (MQ=255)
                                                                                                                                            cGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCAAGAATGAAAGAGATGACGCCGGCTAAATATATAGGTAGAGCTATCACTATGGtt    >  5:145859/1‑151 (MQ=255)
                                                                                                                                             gagaCGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCACGATGGTtg   >  3:866746/1‑151 (MQ=255)
                                                                                                                                              agaCGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCACGATGGTtga  <  6:31325/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCATATCAATCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGAAGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCACGATGGTTGA  >  DH1/2688818‑2689110

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