Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,571,9950TG39.1% 240.7 / 43.2 128Y659D (TAT→GAT) paoCPaoABC aldehyde oxidoreductase Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base T (13/65);  new base G (50/0);  total (63/65)
Fisher's exact test for biased strand distribution p-value = 9.35e-24
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

TGTGTGTGCTGCAGGACGCATCCTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGAGCTGGGCCTGTGCGGCGTGAGCGCGGCTATCGCCAACGCGGTGTATAAC  >  DH1/3571847‑3572144
                                                                                                                                                    |                                                                                                                                                      
tgtgtgtgCTGCAGGACGCATCCTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGAt                                                                                                                                                      >  2:1993016/1‑151 (MQ=255)
    tgtgCTGCAGGACGCATCCTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAgg                                                                                                                                                  >  2:2200733/1‑151 (MQ=255)
      tgctgcAGGACGCATCCTGAATCCGAAAACTGCGCGCAGCCACGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATGACGATATGGCGGGGTCGGAgg                                                                                                                                                  >  2:577062/1‑149 (MQ=255)
       gctgcAGGACGCATCCTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTg                                                                                                                                                <  1:647878/150‑1 (MQ=255)
       gctgcAGGACGCATCCTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAGGTGc                                                                                                                                               >  2:426321/1‑151 (MQ=255)
         tgcAGGACGCATCCTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGTTTAGGGGCCg                                                                                                                                             >  1:2423668/1‑151 (MQ=255)
           cAGGACGCATCCTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAGGGGCCggg                                                                                                                                           >  1:1176929/1‑150 (MQ=255)
            aGGACGCATCCTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGtt                                                                                                                                          <  1:1454736/151‑1 (MQ=255)
              gACGCATCCTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCTTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATTAGGGGGCGGttt                                                                                                                                         >  2:894614/1‑149 (MQ=255)
              gACGCATCCTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCa                                                                                                                                        <  1:1094901/151‑1 (MQ=255)
                cGCATCCTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGt                                                                                                                                                        <  2:1469706/133‑1 (MQ=255)
                 gCATCCTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGc                                                                                                                                     <  1:1647405/151‑1 (MQ=255)
                  cATCCTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGAGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCg                                                                                                                                    <  2:1209573/151‑1 (MQ=255)
                   gtCCTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCa                                                                                                                                        <  1:2410339/145‑1 (MQ=255)
                   gtCCTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAGGGGCCGGGTCa                                                                                                                                        >  2:2410350/2‑146 (MQ=255)
                   aTCCTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAGGGGGCGGGTTACGCgg                                                                                                                                   >  1:50165/1‑151 (MQ=255)
                     ccTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGat                                                                                                                                 <  1:1733939/151‑1 (MQ=255)
                     ccTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGa                                                                                                                                  <  2:2619249/150‑1 (MQ=255)
                       tGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATTAGGTGGCGGGTTATggggcgat                                                                                                                               >  2:16792/1‑147 (MQ=255)
                        gAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGatat                                                                                                                               <  2:809677/150‑1 (MQ=255)
                        gAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGatat                                                                                                                               <  1:1756480/150‑1 (MQ=255)
                          aTCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATccc                                                                                                                            <  1:985496/151‑1 (MQ=255)
                            ccGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGGGGGGGCGGGTCATGCgggttttccc                                                                                                                           >  1:1471211/1‑142 (MQ=255)
                             cGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGCTGAGGTGGCCGGTCATTCggtaatcgcca                                                                                                                          >  2:1218071/1‑141 (MQ=255)
                               aaaaCTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAAc                                                                                                                       <  1:426322/151‑1 (MQ=255)
                                aaaCTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACa                                                                                                                      <  1:505566/151‑1 (MQ=255)
                                 aaCTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACag                                                                                                                     <  2:1635668/151‑1 (MQ=255)
                                  aCTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACagg                                                                                                                    <  2:1999504/151‑1 (MQ=255)
                                    tGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACag                                                                                                                     <  2:1028631/148‑1 (MQ=255)
                                     gcgcgcAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACaggagg                                                                                                                 <  2:399315/151‑1 (MQ=255)
                                       gcgcAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTAt                                                                                                                                                      <  2:486806/112‑1 (MQ=255)
                                         gcAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAGGGGGCGGTTCATGCGGATATCCCAAAAAAGGAGGGGAt                                                                                                            >  2:2360975/1‑151 (MQ=255)
                                          cAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGGGCCGGGTTATTCGGGTATCCCCAAACaggggggggtt                                                                                                            >  2:48398/1‑146 (MQ=255)
                                          cAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAGGGGCCGGGTCATGCGGGTATTCCCAAACaagggggggtt                                                                                                            >  2:2449542/1‑141 (MQ=255)
                                           aGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGAttt                                                                                                          <  1:1993008/151‑1 (MQ=255)
                                           aGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGAttt                                                                                                          <  1:2262788/151‑1 (MQ=255)
                                              cAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATc                                                                                                                              <  2:194282/129‑1 (MQ=255)
                                              cAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAGGGGGCGGGTCATGCGGgtatt                                                                                                                              >  1:194283/1‑128 (MQ=255)
                                               aGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCt                                                                                                      <  1:52716/151‑1 (MQ=255)
                                               aGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAGGGGCCGGTTTATGCGGGTATCCCCAAACAGGAGGGGATTTTTCt                                                                                                      >  2:52716/1‑151 (MQ=255)
                                               aGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAGGGGCCGGTTCATGCGGATATCCCCAAACAGGAGGGGATTTTTCt                                                                                                      >  2:1505/1‑151 (MQ=255)
                                                ggTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTg                                                                                                     <  2:1176930/151‑1 (MQ=255)
                                                   cATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGgat                                                                                                  <  2:1728135/151‑1 (MQ=255)
                                                   cATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATTAGGGGCCGGGTCATTCGGATATCCCCAAAAAGGAGGGGAttagaaccgg                                                                                                   >  1:913983/1‑142 (MQ=255)
                                                    aTTGGCGCAATGACTATGGGCATGGGAGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGgatg                                                                                                 <  2:1191845/151‑1 (MQ=255)
                                                      tGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGat                                                                                                                                 <  1:2319410/118‑1 (MQ=255)
                                                      tGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATTAGGGGGCGGGTCCTGCGGAAAATCCCAAACAAGAGGGGAGTTTTCTTggtgat                                                                                               >  2:2278666/1‑151 (MQ=255)
                                                         cgcAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATAcc                                                                                            <  1:16792/151‑1 (MQ=255)
                                                                   atgggcatgggcGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATAcc                                                                                            <  1:2403771/141‑1 (MQ=255)
                                                                   atgggcatgggcGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAGGTGGCGGGTCATGCGGGTATCCCAAAAAAAGGGGGGGGTTTTCTGGGTGATACCCACCCCacat                                                                                  >  1:367125/1‑151 (MQ=255)
                                                                   atgggcatgggcGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAGGGGCCGGGTCATTCGGGTATCCCCAAAAAGGAGGGGATTTTTCTGGGTGATaaccgatcggaa                                                                                   >  2:2403782/1‑139 (MQ=255)
                                                                     gggcatgggcGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATAtcc                                                                                <  2:2423679/151‑1 (MQ=255)
                                                                     gggcatgggcGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATAtcc                                                                                <  1:48398/151‑1 (MQ=255)
                                                                       gcatgggcGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGGGCCGGTTCCTGCGGATATCCCAAAACAGGAGGGGATTTTCCTGGATGAAACCCACCCCCTAtcctc                                                                              >  2:474497/1‑151 (MQ=255)
                                                                       gcatgggcGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAGGGGCCGGTTCATGCGGGTATCCCAAAACAGGGGGGGGTTTTTCTGGATGATACCGACCCCCTAtcctc                                                                              >  2:676841/1‑151 (MQ=255)
                                                                          tgggcGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGgat                                                                                                  <  2:780292/128‑1 (MQ=255)
                                                                           gggcGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTcccc                                                                           <  2:1471216/150‑1 (MQ=255)
                                                                                cgGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCat                                                                                    <  2:2453795/136‑1 (MQ=255)
                                                                                cgGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGaa                                                                     <  1:1189824/151‑1 (MQ=255)
                                                                                  gCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAgg                                                                   <  1:652043/151‑1 (MQ=255)
                                                                                  gCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAgg                                                                   <  1:2278657/151‑1 (MQ=255)
                                                                                      tGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAGGGGCCGGGTCATTCGGATATCCCAAAAAAGGAGGGGATTTTTCTGGGTGATACCGACCCCATATCCTCCCCGATGAAGGGCa                                                                >  2:2086890/1‑150 (MQ=255)
                                                                                       gATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCaa                                                               <  1:1505/150‑1 (MQ=255)
                                                                                        aTGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGGTGAGGGGCCGGGTCCTGGGGGTATCCCCAAAAAGGGGGGGAATTTTCTTGGGGGTAACAGACCCATATCCTCCCCCATGAAGGGCaaa                                                              >  1:2569446/1‑150 (MQ=255)
                                                                                              gagCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCg                                                       <  1:2415497/151‑1 (MQ=255)
                                                                                               agCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCatat                                                                                  <  2:371795/123‑1 (MQ=255)
                                                                                                 cTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAGGGGCCGGTTCATTCGGATATTCCAAAACAGGAGGGGGTTTTCCTGGGTGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCgggg                                                    >  2:2188887/1‑149 (MQ=255)
                                                                                                  tGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCatat                                                                                  <  1:333170/120‑1 (MQ=255)
                                                                                                  tGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATTAAGGGCCGGCTTATGCGGGTATCCCCCAACAGGAGGGGATTTTTCCGGCTGATACCGACCCCCTTACCTCCCCCATGGAAGCCAAAGgtgttcgggg                                                   >  2:864865/1‑145 (MQ=255)
                                                                                                   ggCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGa                                                   <  1:464130/150‑1 (MQ=255)
                                                                                                   ggCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGa                                                   <  2:1759589/150‑1 (MQ=255)
                                                                                                   ggCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGa                                                   <  2:283018/150‑1 (MQ=255)
                                                                                                   ggCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAGGGGCCGGTTTATGCGGATATCCCCAAAAAGGGGGGGATTTTTCTGGATGATACCGACCCCCTATCCTCCCCGCTGAAGGGCAAAGGTGTCGGTg                                                    >  2:2605301/1‑149 (MQ=255)
                                                                                                     cGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCCAAACAGGAGGGGGTTTTCCTGGGTGATATCGACCCCATATCCTCCCCGATGAAAGCCAAAGGTGTCGGTGAGCt                                                >  1:2142995/1‑151 (MQ=255)
                                                                                                     cGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAGGTGCCGGGTCCTTCGGATATCCCAAAAAAGGAGGGGGTTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAAGGCAAAGGTGTCGGGGAGCt                                                >  2:552776/1‑151 (MQ=255)
                                                                                                          gATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAGGTGCCGGTTCATTCGGATGTCCCAAAAAAGGGGGGGATTTTCCTCGGGGGTACCGACCCCATATTCCCCCCCAATAAGGCCAAAGGTGTCGGTGAGCTGGGcc                                           >  2:2537504/1‑151 (MQ=255)
                                                                                                           aTGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGAGCTGGGCCt                                          <  2:2081081/151‑1 (MQ=255)
                                                                                                            tGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGAGCTggg                                             <  1:2478953/147‑1 (MQ=255)
                                                                                                              aCCTTTTGTGCTACTTCGTTAATCACGATATAGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTTGATGATACCGTCCACATATCCTCCCCGATGAAGGCCAAAGGTTTCGGTGAGCGGGGCCtgtg                                       <  1:474497/151‑1 (MQ=255)
                                                                                                              aCCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGGCGGGTCATGCGGGTATCCCCAAAAAGGAGGGGGTTTTTCTGGGTGATACCGAcccc                                                                                      >  1:1155914/1‑104 (MQ=255)
                                                                                                              aCCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGAcccc                                                                                      <  2:1155915/104‑1 (MQ=255)
                                                                                                              aCCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGAGCTGGGCCtgtg                                       <  2:2487035/151‑1 (MQ=255)
                                                                                                                 gTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGAGCTGGGc                                            <  1:318413/143‑1 (MQ=255)
                                                                                                                 gTTTGGGCTACTTCGTTAATCACGATATGGCGGGGGATGAGGGGCCGGTTCATGCGGATATCCCAAAACAGGAGGGGATTTTTCTGGATGATACCGACCCCATATCCTCCCCCATGAAGGCCAAAGGTGTCGGTGAGCTGGGc                                            >  2:318412/1‑143 (MQ=255)
                                                                                                                    tGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGAGCTGGGCCTGTgcggcg                                  <  1:676840/150‑1 (MQ=255)
                                                                                                                     gggCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGAGCTGGGCCTGTGCGGCGTg                                <  1:2188878/151‑1 (MQ=255)
                                                                                                                      ggCTACTTCGTTAATCACGATATGGCGGGGGATGAGGGGCCGGGTCATGCGGATATCCCAAAACAGGAGGGGATTTTCCTGGGTGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGGGAGCTGGGCCTGTGCGGCGTg                                >  2:1162367/1‑150 (MQ=255)
                                                                                                                        cTACTTCGTTAATCACGATATGGCGGGGGATGAGGGGGCGGGTCATGCGGATCTCCCAAAACAGGGGGTGATTTTCCTGGATGATAACGACCCCATATCCCCCCCCCTGCAGGCCAAAGGGGTCGGGGAGCTGGGCCTgcgtggcgtgggg                             >  2:559175/1‑148 (MQ=255)
                                                                                                                         tACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAgg                                                                   <  1:2186752/112‑1 (MQ=255)
                                                                                                                         tACTTCGTTAATCACGATATGGCGGGGGATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGGGATTTTCCTGGGTGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGGGTCGGGGAGCTGGGCCTGTGCGGCGTGAgcg                            >  2:2190240/1‑151 (MQ=255)
                                                                                                                             tCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGAGCTGGGCCTGTGCGGCGTGAGCGCgg                         <  1:2190231/150‑1 (MQ=255)
                                                                                                                              cGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGAGCTGGGCCTGTGCGGCGTGAGCGCGGCt                       <  2:367124/151‑1 (MQ=255)
                                                                                                                                 tAATCACGATATGGCGGGGGATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGGGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGAGCTGGGCCTGTGCGGCGTGAGCGCGGCTATc                    >  2:114837/1‑151 (MQ=255)
                                                                                                                                       cGATATGGCGGGGGATGAGGTGCCGGTTCATGCGGATATCCCCAAACAGGAGGTGATTTTTCTGGGTGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGAGCTGGGCCTGTGCGGCGTGAGCGCGGCTATCGCCAAc              >  2:1819168/1‑151 (MQ=255)
                                                                                                                                         atatGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGCGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGAGCTGGGCCTGTGCGGCGTGAGCGCGGCTATCGCCAAcgc            >  1:1794194/1‑151 (MQ=255)
                                                                                                                                           atGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATATCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGAGCTGGGCCTGTGCGGCGTGAGCGCGGCTATCGCCAAcgcg           <  2:2143003/150‑1 (MQ=255)
                                                                                                                                            tGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGAGCTGGGCCTGTGCGGCGTGAGCGCGGCTATCGCCAAc              >  1:2543523/1‑146 (MQ=255)
                                                                                                                                            tGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGAGCTGGGCCTGTGCGGCGTGAGCGCGGCTATCGCCAAc              <  2:2543538/146‑1 (MQ=255)
                                                                                                                                            tGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGAGCTGGGCCTGTGCGGCGTGAGCGCGGCTATCGCCAACGCGgt         >  2:695728/1‑151 (MQ=255)
                                                                                                                                                   gTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGAGCTGGGCCTGTGCGGCGTGAGCGCGGCTATCGCCAACGCGGTGTATAAc  >  1:624085/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TGTGTGTGCTGCAGGACGCATCCTGAATCCGAAAACTGCGCGCAGCCAGGTCATTGGCGCAATGACTATGGGCATGGGCGCGGCACTGATGGAGGAGCTGGCGGTGGATGACCGTTTGGGCTACTTCGTTAATCACGATATGGCGGGGTATGAGGTGCCGGTTCATGCGGATATCCCAAAACAGGAGGTGATTTTCCTGGATGATACCGACCCCATATCCTCCCCGATGAAGGCCAAAGGTGTCGGTGAGCTGGGCCTGTGCGGCGTGAGCGCGGCTATCGCCAACGCGGTGTATAAC  >  DH1/3571847‑3572144

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