Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 2,404,019 +G coding (239/861 nt) ydbC ← oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,404,0181.G98.4% 1439.5 / 10.2 378coding (240/861 nt)ydbCoxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (1/5);  new base G (188/184);  total (189/189)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 403 total aligned reads displayed.

CCCCATCATAACGCGCAGGTTAACCACATCCAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACGCAGCACGGTTATAGCGACGTGGCGATCTCGT  >  DH1/2403868‑2404162
                                                                                                                                                       |                                                                                                                                                
ccccATCATAACGCGCAGGTTAACCACATCCAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAAtt                                                                                                                                                   >  8:1514302/1‑151 (MQ=255)
 cccATCATAACGCGCAGGTTAACCACATCCAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTg                                                                                                                                                 >  4:1126405/1‑151 (MQ=255)
 cccATCATAACGCGCAGGTTAACCACATCCAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTg                                                                                                                                                 >  4:836692/1‑151 (MQ=255)
   catcatAACGCGCAGGTTAACCACATCCAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTggt                                                                                                                                               <  7:402929/151‑3 (MQ=255)
    atcatAACGCGCAGGTTAACCACATCCAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTggtc                                                                                                                                              >  5:348258/1‑148 (MQ=255)
    atcatAACGCGCAGGTTAACCACATCCAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACAGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTggtc                                                                                                                                              >  5:195068/1‑148 (MQ=255)
     tcatAACGCGCAGGTTAACCACATCCAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTggtc                                                                                                                                              >  7:537171/1‑147 (MQ=255)
      catAACGCGCAGGTTAACCACATCCAGCACGTCCAGCCCGAGATTACGTAGTTTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGAGCACCAATTTTAGTGACAATTggtc                                                                                                                                              >  3:673985/1‑146 (MQ=255)
      catAACGCGCAGGTTAACCACATCCAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTggtca                                                                                                                                             <  1:1257757/150‑5 (MQ=255)
         aaCGCGCAGGTTAACCACATCCAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGTCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAgggc                                                                                                                                          <  7:738058/151‑3 (MQ=255)
                aGGTTAACCACATCCAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCagag                                                                                                                                   <  2:1826086/151‑1 (MQ=255)
                aGGTTAACCACATCCAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCag                                                                                                                                     >  4:1727175/1‑149 (MQ=255)
                    tAACCACATCCAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTa                                                                                                                                 >  7:822143/1‑149 (MQ=255)
                     aaCCACATCCAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAgg                                                                                                                              >  6:1871156/1‑151 (MQ=255)
                        cacaTCCAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGAt                                                                                                                            >  4:694383/1‑150 (MQ=255)
                        cacaTCCAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGAt                                                                                                                            >  2:127358/1‑150 (MQ=255)
                         acaTCCAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATaa                                                                                                                          <  3:233636/151‑1 (MQ=255)
                             ccAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAAGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAgcg                                                                                                                      <  8:949745/151‑1 (MQ=255)
                               aGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCt                                                                                                                    <  7:1718092/151‑1 (MQ=255)
                               aGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCt                                                                                                                    >  5:1258501/1‑151 (MQ=255)
                                    gTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTcgcg                                                                                                               <  7:1167234/151‑1 (MQ=255)
                                    gTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCCCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTcgcg                                                                                                               >  2:976001/1‑151 (MQ=255)
                                            cGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCt                                                                                                       <  4:1340326/151‑1 (MQ=255)
                                             gagaTTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTg                                                                                                      <  7:1551962/151‑1 (MQ=255)
                                               gaTTACGTAGATTATCGTACACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGAt                                                                                                    <  4:1922292/151‑1 (MQ=255)
                                                aTTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGAtt                                                                                                   <  1:1023230/151‑1 (MQ=255)
                                                aTTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGAt                                                                                                    <  8:1582263/150‑1 (MQ=255)
                                                 ttACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTg                                                                                                  >  6:970027/1‑151 (MQ=255)
                                                  tACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTg                                                                                                  <  1:1937472/150‑1 (MQ=255)
                                                  tACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTg                                                                                                  <  8:1831683/150‑1 (MQ=255)
                                                   aCGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGt                                                                                                >  5:409630/1‑151 (MQ=255)
                                                    cGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTg                                                                                               <  4:1513986/151‑1 (MQ=255)
                                                      tAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGAc                                                                                             >  8:502150/1‑151 (MQ=255)
                                                       aGACTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGAc                                                                                             <  6:830989/150‑1 (MQ=255)
                                                         aTTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTg                                                                                          <  7:330476/151‑1 (MQ=255)
                                                          ttATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGc                                                                                         >  2:228007/1‑151 (MQ=255)
                                                           tATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCATCGCGGATAATCTGATTGGTGACGTGCg                                                                                        >  4:589312/1‑151 (MQ=255)
                                                            aTCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCgg                                                                                       <  1:299040/151‑1 (MQ=255)
                                                            aTCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCg                                                                                        <  6:784677/150‑1 (MQ=255)
                                                                tGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAAATTTAGTGACAATTGGTCACGTCGGCACAGTAAGGATAAAGCGCTTCGCGGATAAGCTCATTGGTGACGTGCGGAcc                                                                                    >  6:694316/1‑150 (MQ=255)
                                                                 gCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCAt                                                                                  >  2:740491/1‑151 (MQ=255)
                                                                 gCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCAt                                                                                  >  5:659426/1‑151 (MQ=255)
                                                                    ccGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATaaa                                                                               <  1:668836/151‑1 (MQ=255)
                                                                     cGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGCCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATaaaa                                                                              <  5:1894954/151‑1 (MQ=255)
                                                                     cGCTTTTTGCAGCTCTACCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATaaaa                                                                              <  2:1733435/151‑1 (MQ=255)
                                                                      gCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAg                                                                             <  1:938742/151‑1 (MQ=255)
                                                                        tttttGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTc                                                                           >  2:1645393/1‑151 (MQ=255)
                                                                        tttttGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTc                                                                           >  3:1785419/1‑151 (MQ=255)
                                                                         ttttGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTc                                                                           <  4:242340/150‑1 (MQ=255)
                                                                         ttttGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCg                                                                          >  3:829807/1‑151 (MQ=255)
                                                                            tGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCAt                                                                                  <  3:510615/140‑1 (MQ=255)
                                                                            tGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCAt                                                                                  >  4:505798/1‑140 (MQ=255)
                                                                              cAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCCCGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGCCGTGCGGACCATAAAAGTCGCTGGt                                                                     >  5:1175354/1‑151 (MQ=255)
                                                                              cAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGt                                                                     <  6:598107/151‑1 (MQ=255)
                                                                                  gctGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCa                                                                 <  7:19888/150‑1 (MQ=255)
                                                                                    tGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAat                                                               >  2:1091189/1‑151 (MQ=255)
                                                                                    tGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAat                                                               >  1:418340/1‑151 (MQ=255)
                                                                                      ccGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAatat                                                             >  2:1421425/1‑151 (MQ=255)
                                                                                      ccGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAatat                                                             <  1:1513358/151‑1 (MQ=255)
                                                                                      ccGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAatat                                                             <  8:728156/151‑1 (MQ=255)
                                                                                        ggAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATg                                                            >  8:898693/1‑150 (MQ=255)
                                                                                         gagaAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGAt                                                          >  6:1470787/1‑151 (MQ=255)
                                                                                           gaAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTg                                                        <  5:1578806/151‑1 (MQ=255)
                                                                                           gaAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTg                                                        <  6:210117/151‑1 (MQ=255)
                                                                                            aaaaTGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGa                                                       <  8:120579/151‑1 (MQ=255)
                                                                                             aaaTGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGAc                                                      <  3:1717272/151‑1 (MQ=255)
                                                                                               aTGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGc                                                    <  7:1424961/151‑1 (MQ=255)
                                                                                                  cGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCa                                                 <  3:1613692/151‑1 (MQ=255)
                                                                                                   gggCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCa                                                 >  5:1417252/1‑150 (MQ=255)
                                                                                                      cAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGATTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGc                                             <  4:1590836/151‑1 (MQ=255)
                                                                                                       aaCCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGcc                                            >  6:173609/1‑151 (MQ=255)
                                                                                                        aCCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGcc                                            >  7:1621064/1‑150 (MQ=255)
                                                                                                        aCCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCa                                           <  8:389409/151‑1 (MQ=255)
                                                                                                        aCCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCa                                           <  1:1636527/151‑1 (MQ=255)
                                                                                                         ccAGGATGCGTCCTCTCCACGCCGCGCAGCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGcc                                            >  5:872918/1‑149 (MQ=255)
                                                                                                          cAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAgc                                         <  3:1120180/151‑1 (MQ=255)
                                                                                                           aGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAgcg                                        <  2:377818/151‑1 (MQ=255)
                                                                                                           aGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAgcg                                        >  8:1311815/1‑151 (MQ=255)
                                                                                                                gCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTc                                    <  8:1670231/150‑1 (MQ=255)
                                                                                                                 cGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTccc                                  >  2:631531/1‑149 (MQ=255)
                                                                                                                    ccTCTCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACgc                                >  4:988727/1‑150 (MQ=255)
                                                                                                                     ctctCCACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACgc                                >  7:1868001/1‑149 (MQ=255)
                                                                                                                         ccACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAACGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACGCAGCACg                          <  5:754646/151‑1 (MQ=255)
                                                                                                                         ccACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACGCAGCAc                           <  4:339792/150‑1 (MQ=255)
                                                                                                                         ccACGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACGCAGCACg                          >  3:1876161/1‑151 (MQ=255)
                                                                                                                           aCGCCGCGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACGCAGCACGGt                        >  6:1355181/1‑151 (MQ=255)
                                                                                                                              ccgcGCACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACGCAGCACGGTtat                     >  6:1444910/1‑151 (MQ=255)
                                                                                                                                gcgcACCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACGCAGCACGGTTATAg                   <  1:1765166/151‑1 (MQ=255)
                                                                                                                                    aCCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGGCGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACGCAGCACGGTTATAGCGAc               <  1:1709676/151‑1 (MQ=255)
                                                                                                                                    aCCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACGCCGCACGGTTATAGCGAc               >  1:1356914/1‑151 (MQ=255)
                                                                                                                                    aCCAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACGCAGCACGGTTATAGCGAc               >  6:296016/1‑151 (MQ=255)
                                                                                                                                     ccAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGCCGCCCAATGCCAGCGCCTCACGCAGCACGGTTATAGCGACg              <  8:143480/151‑1 (MQ=255)
                                                                                                                                     ccAATTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACGCAGCACGGTTATAGCGACg              >  2:1906649/1‑151 (MQ=255)
                                                                                                                                       aaTTTTAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACGCAGCACGGTTATAGCGACGTg            <  8:679203/151‑1 (MQ=255)
                                                                                                                                         ttttAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACGCAGCACGGTTATAGCGACGTGGc          >  6:1203197/1‑151 (MQ=255)
                                                                                                                                            tAGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACGCAGCACGGTTATAGCGACGTGGCGa        <  3:40008/150‑1 (MQ=255)
                                                                                                                                             aGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACGCAGCACGGTTATAGCGACGTGGCGAtc      >  2:590301/1‑151 (MQ=255)
                                                                                                                                             aGTGACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACGCAGCACGGTTATAGCGACGTGGCGAtc      >  7:1277468/1‑151 (MQ=255)
                                                                                                                                                gACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCCGCGCCTCACGCAGCACGGTTATAGCGACGTGGCGATCTCg   >  1:813222/1‑151 (MQ=255)
                                                                                                                                                gACAATTGGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACGCAGCACGGTTATAGCGACGTGGCGAtctc    <  3:1296183/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                
CCCCATCATAACGCGCAGGTTAACCACATCCAGCACGTCCAGCCCGAGATTACGTAGATTATCGTGCACCGCTTTTTGCAGCTCTGCCGGAGAAAATGCGGGCAACCAGGATGCGTCCTCTCCACGCCGCGCACCAATTTTAGTGACAATTGTCAGGTCGTCAGAGTAAGGATAAAGCGCTTCGCGGATAATCTGATTGGTGACGTGCGGACCATAAAAGTCGCTGGTATCAATATGATTGACGCCCAATGCCAGCGCCTCACGCAGCACGGTTATAGCGACGTGGCGATCTCGT  >  DH1/2403868‑2404162

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: