Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH14,194,5020AG54.5% 174.7 / inf 541Q195Q (CAA→CAG) phnIalpha‑D‑ribose 1‑methylphosphonate 5‑triphosphate synthase subunit PhnI
Reads supporting (aligned to +/- strand):  ref base A (118/127);  new base G (134/160);  total (252/289)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-01

Only 100 of 548 total aligned reads displayed.

GACAGCGAACAGCAGCCGTCGCCGCACGTTTTCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGAAGAGCTGGGATTTGCGGTGAACGTC  >  DH1/4194353‑4194652
                                                                                                                                                     |                                                                                                                                                      
 aCAGCGAACAGCAGCCGTCGCCGCACGTTTTCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAAtt                                                                                                                                                      <  4:773046/151‑1 (MQ=255)
 aCAGCGAACAGCAGCCGTCGCCGCACGTTTTCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAAtt                                                                                                                                                      >  7:1617028/1‑151 (MQ=255)
  cAGCGAACAGCAGCCGTCGCCGCACGTTTTCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTg                                                                                                                                                     <  4:1938022/151‑1 (MQ=255)
  cAGCGAACAGCAGCCGTCGCCGCACGTTTTCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTg                                                                                                                                                     >  4:1530300/1‑151 (MQ=255)
  cAGCGAACAGCAGCCGTCGCCGCACGTTTTCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTg                                                                                                                                                     <  1:22867/151‑1 (MQ=255)
    gCGAACAGCAGCCGGCGCCGCACGTTTTCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGAt                                                                                                                                                   <  7:1288975/151‑1 (MQ=255)
       aaCAGCAGCCGTCGCCGCACGTTTTCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATgcg                                                                                                                                                >  7:83095/1‑151 (MQ=255)
          agcagcCGTCGCCGCACGTTTTCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCgcgg                                                                                                                                             >  2:1920380/1‑151 (MQ=255)
           gcagcCGTCGCCGCACGTTTTCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCgcggc                                                                                                                                            <  5:695049/151‑1 (MQ=255)
               ccGTCGCCGCACGTTTTCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGcgacg                                                                                                                                        >  1:483844/1‑151 (MQ=255)
                cGTCGCCGCACGTTTTCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGcgacga                                                                                                                                       <  2:977802/151‑1 (MQ=255)
                   cgccgcACGTTTTCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAg                                                                                                                                     <  7:734738/150‑1 (MQ=255)
                        cACGTTTTCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTAtt                                                                                                                               >  2:1359501/1‑151 (MQ=255)
                         aCGTTTTCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTAttt                                                                                                                              <  6:1716097/151‑1 (MQ=255)
                         aCGTTTTCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTAttt                                                                                                                              >  1:1898770/1‑151 (MQ=255)
                           gTTTTCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTgc                                                                                                                            <  5:1449673/151‑1 (MQ=255)
                             tttCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTgct                                                                                                                           <  4:915150/150‑1 (MQ=255)
                               tCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTgctggc                                                                                                                        >  8:70273/1‑151 (MQ=255)
                                   ccTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTgctggcgctg                                                                                                                    >  3:1372274/1‑151 (MQ=255)
                                   ccTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTgctggcgctg                                                                                                                    <  6:391846/151‑1 (MQ=255)
                                   ccTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTgctggcgctg                                                                                                                    >  6:606742/1‑151 (MQ=255)
                                      gctgGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTgctggcgctggc                                                                                                                  <  5:491242/150‑1 (MQ=255)
                                        gggCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTAGCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTa                                                                                                               <  7:1209425/150‑1 (MQ=255)
                                         ggCGCGTCAGGGGCGGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCt                                                                                                                <  5:281638/149‑1 (MQ=255)
                                            gcgTCAGGGGCTGGCGAAGTTTGAAGAGGATCGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTcc                                                                                                           <  3:1449271/151‑1 (MQ=255)
                                            gcgTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTcc                                                                                                           >  6:1212870/1‑151 (MQ=255)
                                              gTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCAc                                                                                                         >  1:1011269/1‑151 (MQ=255)
                                               tCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCAcc                                                                                                        >  8:883374/1‑151 (MQ=255)
                                                  ggggCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCac                                                                                                     >  4:1192708/1‑150 (MQ=255)
                                                   gggCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAAc                                                                                                    <  8:1644794/151‑1 (MQ=255)
                                                      cTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAAcgcg                                                                                                 >  7:435677/1‑151 (MQ=255)
                                                        ggCTAAGTTTGAAGTGGATAGCGGCGCAGAGCCGGATGACTTCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTGCCGTTTGCAGCAATTGAAGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAAcgcg                                                                                                 <  1:608212/149‑1 (MQ=255)
                                                        ggCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGt                                                                                               <  2:1441196/151‑1 (MQ=255)
                                                         gCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGtt                                                                                              <  6:93105/151‑1 (MQ=255)
                                                          cGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTa                                                                                             <  7:965277/151‑1 (MQ=255)
                                                            aaGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATCGGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCCGTTAc                                                                                            >  6:1975657/1‑150 (MQ=255)
                                                             aGTTTGACGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACCCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCACCGCGGTTACgg                                                                                          <  8:140537/151‑1 (MQ=255)
                                                                ttGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGgc                                                                                        <  4:1637300/150‑1 (MQ=255)
                                                                 tGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGgcg                                                                                       <  3:1913485/150‑1 (MQ=255)
                                                                  gAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGgcgc                                                                                      <  6:527897/150‑1 (MQ=255)
                                                                    agagGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAAt                                                                                   >  8:782630/1‑151 (MQ=255)
                                                                         aTAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCAccc                                                                              >  2:813083/1‑151 (MQ=255)
                                                                           aGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCg                                                                             <  3:183327/150‑1 (MQ=255)
                                                                            gcggcgCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGtt                                                                           <  7:1416868/151‑1 (MQ=255)
                                                                             cggcgCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTc                                                                          <  6:1852645/151‑1 (MQ=255)
                                                                             cggcgCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTc                                                                          >  3:1170897/1‑151 (MQ=255)
                                                                                cgcACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGc                                                                        >  7:544224/1‑150 (MQ=255)
                                                                                cgcACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGc                                                                        >  5:834487/1‑150 (MQ=255)
                                                                                cgcACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCa                                                                       <  3:1531210/151‑1 (MQ=255)
                                                                                  cacaGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAgg                                                                     <  1:995535/151‑1 (MQ=255)
                                                                                   acaGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGc                                                                    >  3:1605113/1‑151 (MQ=255)
                                                                                   acaGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGc                                                                    <  2:434311/151‑1 (MQ=255)
                                                                                     aGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCg                                                                   >  7:1498842/1‑150 (MQ=255)
                                                                                     aGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCg                                                                   >  8:1836867/1‑150 (MQ=255)
                                                                                          gATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATcc                                                             <  5:584334/151‑1 (MQ=255)
                                                                                             gACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGc                                                           >  2:1831946/1‑150 (MQ=255)
                                                                                             gACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGc                                                           >  8:539744/1‑150 (MQ=255)
                                                                                              aCATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAg                                                         <  6:1522950/151‑1 (MQ=255)
                                                                                               cATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGc                                                        <  8:331927/151‑1 (MQ=255)
                                                                                                aTCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCg                                                       <  3:677415/151‑1 (MQ=255)
                                                                                                aTCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCg                                                       <  1:74539/151‑1 (MQ=255)
                                                                                                 tCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCgg                                                      >  2:510521/1‑151 (MQ=255)
                                                                                                  cACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCgg                                                      <  2:1144151/150‑1 (MQ=255)
                                                                                                   aCCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCt                                                    <  4:564939/151‑1 (MQ=255)
                                                                                                   aCCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCt                                                    >  6:276345/1‑151 (MQ=255)
                                                                                                     ccGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCtat                                                  <  1:1180513/151‑1 (MQ=255)
                                                                                                      cGCACGCCGCCGGTTTTCCGCTGCTCACGCTCGTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCtata                                                 <  8:1266629/151‑1 (MQ=255)
                                                                                                       gCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTatat                                                <  3:539006/151‑1 (MQ=255)
                                                                                                       gCACGCCGCCGGTTCACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTatat                                                >  8:1140863/1‑151 (MQ=255)
                                                                                                         aCGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCTCCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCg                                              >  2:1556159/1‑151 (MQ=255)
                                                                                                         aCGCCGCCGGGTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCg                                              <  3:406678/151‑1 (MQ=255)
                                                                                                            ccgccgGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACg                                           >  7:325315/1‑151 (MQ=255)
                                                                                                            ccgccgGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACg                                           >  7:211408/1‑151 (MQ=255)
                                                                                                                 ggTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCg                                      <  7:1038090/151‑1 (MQ=255)
                                                                                                                      aCCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATtgt                                 <  1:1610941/151‑1 (MQ=255)
                                                                                                                        cccTGCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGc                               <  8:1745113/151‑1 (MQ=255)
                                                                                                                         ccTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGcc                              <  5:1178090/151‑1 (MQ=255)
                                                                                                                            gCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCgg                            <  2:320289/150‑1 (MQ=255)
                                                                                                                            gCTCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGa                           <  4:1802590/151‑1 (MQ=255)
                                                                                                                              tCACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGa                           >  5:1291450/1‑149 (MQ=255)
                                                                                                                               cACGCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGAaga                        <  4:1275947/151‑1 (MQ=255)
                                                                                                                                  gCTCCTCCCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGAAGAGCt                     <  5:379498/151‑1 (MQ=255)
                                                                                                                                    tcctccCGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATACGCAGCGGCTATATCGACGTTTCGATTGTGCCGGAAGCTCTgg                   >  2:215048/1‑151 (MQ=255)
                                                                                                                                     cctccCGTTTGCAGCAGTTGATGCGCGGCGACGAAGTCTATTTGCTGGCGCTGGCCTACTCCACACAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGAGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTTCCGGAAGAGCAggg                  >  5:1997070/1‑151 (MQ=255)
                                                                                                                                        cccGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGAAGAGCTGGGAtt               <  2:1699290/151‑1 (MQ=255)
                                                                                                                                        cccGTTTGCAGCAGTTGATCCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGCCGATTGTGCCGGAAGAGCTGGGAt                <  4:295737/150‑1 (MQ=255)
                                                                                                                                         ccGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGAAGAGCTGGGAttt              >  1:131187/1‑151 (MQ=255)
                                                                                                                                          cGTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGTGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGAAGAGCTGGGATTTg             >  6:71377/1‑151 (MQ=255)
                                                                                                                                           gTTTGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGAAGAGCTGGGATTTGc            >  3:860735/1‑151 (MQ=255)
                                                                                                                                           gTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCTACGTGTCGATTGTGCCGGAAGAGCTGGGATTTGc            <  1:378913/151‑1 (MQ=255)
                                                                                                                                              tGCAGCAGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGAAGAGCTGGGATTTGCGGt         >  4:850047/1‑151 (MQ=255)
                                                                                                                                               gcagcaGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGAAGAGCTGGGATTTGCGGt         <  7:1828000/150‑1 (MQ=255)
                                                                                                                                               gcagcaATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGAAGAGCTGGGATTTGCGGTg        <  5:1965423/151‑1 (MQ=255)
                                                                                                                                                 agcaGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGGTGGCCTACTCCAGCCAACGCGGTTTCGGGCGGAATCTCCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGAAGAGCTGGGATTTGCTGTGaa      <  4:383043/151‑1 (MQ=255)
                                                                                                                                                 agcaATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGAAGAGCTGGGATTTGCGGTGa       >  5:1619989/1‑150 (MQ=255)
                                                                                                                                                  gcaGTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGAAGAGCTGGGATTTGCGGTGAAc     <  3:599363/151‑1 (MQ=255)
                                                                                                                                                  gcaATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGAAGAGCTGGGATTTGCGGTGAAc     <  1:1360551/151‑1 (MQ=255)
                                                                                                                                                    aaTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGAAGAGCTGGGATTTGCGGTGAACGt   >  7:1319185/1‑151 (MQ=255)
                                                                                                                                                    aaTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGAAGAGCTGGGATTTGCGGTGAACGt   <  5:1827196/151‑1 (MQ=255)
                                                                                                                                                     aTTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGAAGAGCTGGGATTTGCGGTGAACGTc  <  4:1011726/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GACAGCGAACAGCAGCCGTCGCCGCACGTTTTCAGCCTGCTGGCGCGTCAGGGGCTGGCGAAGTTTGAAGAGGATAGCGGCGCACAGCCGGATGACATCACCCGCACGCCGCCGGTTTACCCCTGCTCACGCTCCTCCCGTTTGCAGCAATTGATGCGCGGCGACGAAGGCTATTTGCTGGCGCTGGCCTACTCCACCCAACGCGGTTACGGGCGCAATCACCCGTTCGCAGGCGAGATCCGCAGCGGCTATATCGACGTGTCGATTGTGCCGGAAGAGCTGGGATTTGCGGTGAACGTC  >  DH1/4194353‑4194652

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: ‑