Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 3,154,902 A→T I274F (ATC→TTC)  kdpB → potassium‑transporting ATPase B chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,154,9020AT97.9% 1333.3 / 9.0 386I274F (ATC→TTC) kdpBpotassium‑transporting ATPase B chain
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base T (184/193);  total (187/199)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

Only 100 of 386 total aligned reads displayed.

GCTGATTGCCCTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGATCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAGGCGTCGGAGTTTA  >  DH1/3154754‑3155049
                                                                                                                                                    |                                                                                                                                                   
gctgATTGCCCTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTc                                                                                                                                                   <  4:868270/151‑1 (MQ=255)
 ctgATTGCCCTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCCCTATTGGCGGCCTGTTGTCAGCGTTCg                                                                                                                                                  >  3:1457803/1‑151 (MQ=255)
   gagtGCCCTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGc                                                                                                                                                <  1:1538664/148‑1 (MQ=255)
   gATTGCCCTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCg                                                                                                                                                  <  4:1298037/149‑1 (MQ=255)
   gATTGCCCTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCg                                                                                                                                                  <  8:1093776/149‑1 (MQ=255)
    aTTGCCCTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGCTCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCg                                                                                                                                               >  2:1542810/1‑151 (MQ=255)
      tGCCCTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTc                                                                                                                                             >  4:1684140/1‑151 (MQ=255)
       gCCCTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTc                                                                                                                                             <  6:1193988/150‑1 (MQ=255)
          cTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCg                                                                                                                                         >  2:246367/1‑151 (MQ=255)
              cTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGAt                                                                                                                                     >  5:1510596/1‑151 (MQ=255)
               tATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGATCGGCGTCGCCGGGATg                                                                                                                                    >  5:1826179/1‑151 (MQ=255)
                aTCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATg                                                                                                                                    <  1:1732936/150‑1 (MQ=255)
                aTCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATg                                                                                                                                    >  3:1361806/1‑150 (MQ=255)
                aTCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGGGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGAt                                                                                                                                     <  2:969803/149‑1 (MQ=255)
                  cGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGc                                                                                                                                 <  4:1421292/151‑1 (MQ=255)
                   gTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGcc                                                                                                                                >  3:1549843/1‑151 (MQ=255)
                   gTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGcc                                                                                                                                <  8:338071/151‑1 (MQ=255)
                      tttttACTGGCAACCGCCACGCTGTGGCGGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTGTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCgggt                                                                                                                             >  3:645279/1‑150 (MQ=255)
                       ttttACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCg                                                                                                                               <  8:1617125/148‑1 (MQ=255)
                           aCTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCt                                                                                                                         <  6:914010/150‑1 (MQ=255)
                           aCTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCt                                                                                                                         <  4:606872/150‑1 (MQ=255)
                              ggCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGgc                                                                                                                     <  1:1623330/151‑1 (MQ=255)
                                  aCCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGa                                                                                                                 <  5:887936/151‑1 (MQ=255)
                                    cGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAAt                                                                                                               <  5:582714/151‑1 (MQ=255)
                                     gCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGGTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGCAtg                                                                                                              >  7:1551731/1‑151 (MQ=255)
                                     gCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGGCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGATCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAAtg                                                                                                              <  7:1188043/151‑1 (MQ=255)
                                       cACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAAtgtg                                                                                                            >  7:1415655/1‑151 (MQ=255)
                                        aCGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAAtgtg                                                                                                            >  4:312848/1‑150 (MQ=255)
                                        aCGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGAAGAGCCGGATGCTAGGCGCGAAtg                                                                                                              >  4:392129/1‑148 (MQ=255)
                                          gCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGAtt                                                                                                         >  6:1428042/1‑151 (MQ=255)
                                          gCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGAtt                                                                                                         <  4:1189111/151‑1 (MQ=255)
                                               ggCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGccac                                                                                                    >  8:918550/1‑151 (MQ=255)
                                                   gTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGc                                                                                                >  3:251297/1‑151 (MQ=255)
                                                      ttCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGGACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCgg                                                                                              <  8:758554/150‑1 (MQ=255)
                                                       tCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCTTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACAAGCGGAc                                                                                            >  5:1397870/1‑151 (MQ=255)
                                                        ccGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACg                                                                                           <  7:1726267/151‑1 (MQ=255)
                                                         cgcgTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGt                                                                                          <  3:1169009/151‑1 (MQ=255)
                                                         cgcgTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGt                                                                                          <  7:125555/151‑1 (MQ=255)
                                                          gcgTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTg                                                                                         >  8:492468/1‑151 (MQ=255)
                                                            gTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCa                                                                                       <  3:1700700/151‑1 (MQ=255)
                                                             tGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGc                                                                                        >  4:1755419/1‑149 (MQ=255)
                                                              ggggCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAg                                                                                      <  8:65782/150‑1 (MQ=255)
                                                              ggggCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGt                                                                                     <  5:380317/151‑1 (MQ=255)
                                                                ggCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTg                                                                                   <  4:200259/151‑1 (MQ=255)
                                                                ggCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTg                                                                                   <  2:1808723/151‑1 (MQ=255)
                                                                 gCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCTCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGa                                                                                  <  2:1611247/151‑1 (MQ=255)
                                                                  gggTAATGCAGGCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGaa                                                                                 <  5:289686/150‑1 (MQ=255)
                                                                   ggTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGa                                                                                  >  7:1040983/1‑149 (MQ=255)
                                                                    gTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGc                                                                               <  4:661943/151‑1 (MQ=255)
                                                                    gTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGc                                                                               <  3:181270/151‑1 (MQ=255)
                                                                     tAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCg                                                                              >  8:1735162/1‑151 (MQ=255)
                                                                      aaTGCAGTCAGCGGAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGCGATTGCCACCAGCGGACGTGCAGTTGAAGCgg                                                                             <  3:1624188/151‑1 (MQ=255)
                                                                       aTGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGc                                                                            <  5:966291/151‑1 (MQ=255)
                                                                        tGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGc                                                                            >  6:770970/1‑150 (MQ=255)
                                                                        tGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGc                                                                            >  7:997990/1‑150 (MQ=255)
                                                                         gCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCa                                                                           <  8:1414729/150‑1 (MQ=255)
                                                                             tCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGCCGTGCAGTTGAAGCGGCAGGtgc                                                                      >  6:499725/1‑150 (MQ=255)
                                                                              cAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCACAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGtgac                                                                     >  7:1498192/1‑151 (MQ=255)
                                                                               aGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGtgacg                                                                    <  1:795631/151‑1 (MQ=255)
                                                                               aGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGtgacg                                                                    <  7:53089/151‑1 (MQ=255)
                                                                                gCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGtgacgt                                                                   <  7:416559/151‑1 (MQ=255)
                                                                                gCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGtgacgt                                                                   >  1:1149551/1‑151 (MQ=255)
                                                                                    aaCGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGtgacgttgac                                                               >  2:1261923/1‑151 (MQ=255)
                                                                                    aaCGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGtgacgttgac                                                               <  8:256970/151‑1 (MQ=255)
                                                                                      cGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGtgacgttgacgt                                                             >  6:220438/1‑151 (MQ=255)
                                                                                       gggaCTGGTGGCGCTGCTGGTCTGGCTGATCCCAACGACTATTGGCGGCCTGTTGTCAGCGTTCGGGGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTgacgttgacgtt                                                            <  5:1277987/148‑1 (MQ=255)
                                                                                             ggtggCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGGCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTa                                                      <  6:1560628/151‑1 (MQ=255)
                                                                                              gtggCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTAc                                                     <  1:394190/151‑1 (MQ=255)
                                                                                                ggCGCTGCTGGTCTGTCTGCTCCCAACCCCTATTGGCGGCCTGTTGTGAGCGTGCGGGGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATGGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGGGACGTTGACGTTCTGCTACTg                                                   <  7:1365274/151‑1 (MQ=255)
                                                                                                ggCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTg                                                   >  3:127996/1‑151 (MQ=255)
                                                                                                ggCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTg                                                   <  3:619977/151‑1 (MQ=255)
                                                                                                  cgcTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTgg                                                  >  6:1760923/1‑150 (MQ=255)
                                                                                                   gctgctGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGAt                                                <  2:6764/151‑1 (MQ=255)
                                                                                                   gctgctGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGAt                                                <  1:607541/151‑1 (MQ=255)
                                                                                                    ctgctgGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGAt                                                >  1:679277/1‑150 (MQ=255)
                                                                                                              gtctgATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCAc                                     >  5:779141/1‑151 (MQ=255)
                                                                                                                ctgATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCAcc                                    >  6:348867/1‑150 (MQ=255)
                                                                                                                 tgATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCAt                                  >  6:1061691/1‑151 (MQ=255)
                                                                                                                 tgATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCAt                                  >  3:530177/1‑151 (MQ=255)
                                                                                                                  gNTCCCAACCACTATNGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCAc                                     <  3:565914/147‑1 (MQ=255)
                                                                                                                  gATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATc                                 <  5:1690536/151‑1 (MQ=255)
                                                                                                                     cccAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATcac                               <  8:150988/150‑1 (MQ=255)
                                                                                                                      ccAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCacac                             >  6:964720/1‑151 (MQ=255)
                                                                                                                       cAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACt                            <  1:1038328/151‑1 (MQ=255)
                                                                                                                        aaCCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACCCTc                           >  6:1641934/1‑151 (MQ=255)
                                                                                                                          ccACTATTGGCGGCCTGGTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCgg                         <  3:800281/151‑1 (MQ=255)
                                                                                                                            aCTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGt                        >  5:1004361/1‑150 (MQ=255)
                                                                                                                             cTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTa                       >  1:427345/1‑150 (MQ=255)
                                                                                                                              tATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAAc                     >  8:1325932/1‑151 (MQ=255)
                                                                                                                               aTTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAAcc                    >  7:673196/1‑151 (MQ=255)
                                                                                                                               aTTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAAcc                    <  5:1155752/151‑1 (MQ=255)
                                                                                                                                     ggcCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAgg              >  7:580480/1‑151 (MQ=255)
                                                                                                                                       ccTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAGGCg            <  5:476309/151‑1 (MQ=255)
                                                                                                                                        cTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAGGCgg           >  7:1808373/1‑150 (MQ=255)
                                                                                                                                          gttgtCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAGGCGt           <  5:700249/149‑1 (MQ=255)
                                                                                                                                            tgtCAGCGTGCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAGGCGTCgg        <  7:223814/150‑1 (MQ=255)
                                                                                                                                             gtCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAGGCGTCgg        >  3:303304/1‑149 (MQ=255)
                                                                                                                                             gtCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAGGCGTCGGa       <  1:1379434/150‑1 (MQ=255)
                                                                                                                                               cAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAAGCGTCGGAGtt    >  1:1298243/1‑151 (MQ=255)
                                                                                                                                                aGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAGGCGTCGGAGttt   <  2:1169701/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCTGATTGCCCTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGATCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAGGCGTCGGAGTTTA  >  DH1/3154754‑3155049

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