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,9020AT98.0% 1365.1 / 9.5 400I274F (ATC→TTC) kdpBpotassium‑transporting ATPase B chain
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base T (187/204);  total (190/210)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-01

Only 100 of 400 total aligned reads displayed.

GCTGATTGCCCTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGATCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAGGCGTCGGAGTTTA  >  DH1/3154754‑3155049
                                                                                                                                                    |                                                                                                                                                   
  tgATTGCCCTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCgg                                                                                                                                                 >  8:1076690/1‑151 (MQ=255)
   gATTGCCCTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCg                                                                                                                                                  <  7:224763/149‑1 (MQ=255)
   gATTGCCCTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCg                                                                                                                                                  <  3:793580/149‑1 (MQ=255)
     ttGCCCTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGt                                                                                                                                              >  1:104069/1‑151 (MQ=255)
      tGCCCTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGt                                                                                                                                              <  1:213884/150‑1 (MQ=255)
      tGCCCTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGGGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTc                                                                                                                                             >  1:1119470/1‑151 (MQ=255)
        cccTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGc                                                                                                                                           <  8:1144445/151‑1 (MQ=255)
        cccTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGc                                                                                                                                           <  5:581119/151‑1 (MQ=255)
         ccTGACTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGcc                                                                                                                                          >  8:815215/1‑151 (MQ=255)
             aCTATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGa                                                                                                                                      <  8:230831/151‑1 (MQ=255)
               tATCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATg                                                                                                                                    >  3:1106247/1‑151 (MQ=255)
                aTCGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGAt                                                                                                                                     <  1:1249743/149‑1 (MQ=255)
                  cGTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAg                                                                                                                                  <  7:1648046/150‑1 (MQ=255)
                   gTCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGATCGGCGTCGCCGGGATGAGcc                                                                                                                                <  1:310866/151‑1 (MQ=255)
                    tCTTTTTACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCg                                                                                                                               >  2:1148991/1‑151 (MQ=255)
                      tttttACTGGCAACCGCCACGCTGTGGCCGTTTTCCCCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGCTGCCGCTGCTGGTGTGTCTGATCCCAACCACTATTGCCGGCCTGTTGTCAGGGTTCGGCGTCGCCGGGATGAGCCggc                                                                                                                             >  8:1624000/1‑150 (MQ=255)
                       ttttACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGa                                                                                                                             <  5:1413553/150‑1 (MQ=255)
                       ttttACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGGCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCgg                                                                                                                              <  7:532997/149‑1 (MQ=255)
                        tttACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGAt                                                                                                                            <  6:54835/150‑1 (MQ=255)
                        tttACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGAt                                                                                                                            <  5:276283/150‑1 (MQ=255)
                        tttACTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATg                                                                                                                           >  3:1419931/1‑151 (MQ=255)
                           aCTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATTCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTa                                                                                                                        <  4:1351132/151‑1 (MQ=255)
                           aCTGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCt                                                                                                                         >  8:672010/1‑150 (MQ=255)
                           aCTGGCAACCGCCACGCTGTGGCCGGTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCt                                                                                                                         <  4:1063684/150‑1 (MQ=255)
                             tGGCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAgg                                                                                                                      >  4:940570/1‑151 (MQ=255)
                               gCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGgcg                                                                                                                    >  4:1913036/1‑151 (MQ=255)
                               gCAACCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGgcg                                                                                                                    >  2:1632965/1‑151 (MQ=255)
                                 caccgagaCGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGGCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTGCGGCGGCGCCGGGATGAGCCGGATGCTATGcgcg                                                                                                                  <  4:1222424/144‑1 (MQ=255)
                                  aCCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGcgcg                                                                                                                  >  8:1741699/1‑150 (MQ=255)
                                  aCCGCCACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGcgcg                                                                                                                  >  5:1558575/1‑150 (MQ=255)
                                       cACGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAAtgtg                                                                                                            >  1:476196/1‑151 (MQ=255)
                                         cGCTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGAt                                                                                                          >  3:1806282/1‑151 (MQ=255)
                                           gtgtgGACGTTTTCCGCGTGTGGCTGTAATGCAGTCAGCGTAACGGTACTGTTGGCGCTGNNGGTCTNNNTNNNNNNANNNNNNATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTg                                                                                                        <  4:758278/150‑1 (MQ=255)
                                           cTGTGGCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTg                                                                                                        >  5:1514916/1‑151 (MQ=255)
                                              gggCCGTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGcca                                                                                                     <  7:1028365/150‑1 (MQ=255)
                                                   gTTTTCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTATTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGc                                                                                                >  5:1472344/1‑151 (MQ=255)
                                                   gTTTTCCGCGTGGCGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTTGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGc                                                                                                <  4:1866144/151‑1 (MQ=255)
                                                    ttttCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCg                                                                                               >  2:229586/1‑151 (MQ=255)
                                                    ttttCCGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCg                                                                                               >  4:636920/1‑151 (MQ=255)
                                                        ccGCGTGGGGCTTTAATGCAGTCAGCGTAACGGTACTTGTGGAGCTGCTGGTCTTTCTGATCCGAACCACTATTTGCGGCATGTTGTCATCGTTCGGCGTCTCCGGTATGAGCCGGATGCTAGGCTCTAATGTGATTGCCACCAGCGGAc                                                                                            <  6:1785364/150‑1 (MQ=255)
                                                        ccGCGTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGAc                                                                                            <  1:1466556/150‑1 (MQ=255)
                                                         cgcgTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGt                                                                                          >  2:557123/1‑151 (MQ=255)
                                                         cgcgTGGGGCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGt                                                                                          >  7:1858540/1‑151 (MQ=255)
                                                               ggggggTAATGCATTCATCGTAACGGTACTGGTGGCGCTGCTGGTCTTTCTGATGCCAACGACTATTGGGGGCCTGTTGTCAGAGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGTCGCGAATGTGATTGCCACCAGCGGAAGTGCAGt                                                                                     <  5:894073/150‑1 (MQ=255)
                                                               gggCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGtt                                                                                    >  8:580541/1‑151 (MQ=255)
                                                                ggCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTg                                                                                   <  2:1038733/151‑1 (MQ=255)
                                                                 gCGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGa                                                                                  <  4:528505/151‑1 (MQ=255)
                                                                  cGGTAATGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCCGCCTGTTGTCAGCCTTCGGCGTCGCCGGGAAGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGa                                                                                  >  6:1657101/1‑150 (MQ=255)
                                                                      aaTGCAGTCAGCGTACCGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCgg                                                                             <  4:341061/151‑1 (MQ=255)
                                                                      aaTGCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTAATGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCgg                                                                             >  6:150292/1‑151 (MQ=255)
                                                                         gCAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCa                                                                           <  5:713148/150‑1 (MQ=255)
                                                                          cAGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAgg                                                                         >  3:232000/1‑151 (MQ=255)
                                                                           aGTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGATCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGt                                                                        >  2:1424847/1‑151 (MQ=255)
                                                                            gTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGtg                                                                       <  5:1946287/151‑1 (MQ=255)
                                                                            gTCAGCGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGCCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGtg                                                                       >  8:332999/1‑151 (MQ=255)
                                                                                 cGTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGtgacgtt                                                                  >  3:1713530/1‑151 (MQ=255)
                                                                                  gTAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGtgacgttg                                                                 >  5:1878996/1‑151 (MQ=255)
                                                                                   tAACGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGtgacgttg                                                                 >  6:227575/1‑150 (MQ=255)
                                                                                      cGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGtgacgttgacgt                                                             >  3:1507773/1‑151 (MQ=255)
                                                                                      cGGTACTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGtgacgttgacgt                                                             >  6:1109093/1‑151 (MQ=255)
                                                                                          aCTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTg                                                         >  5:128160/1‑151 (MQ=255)
                                                                                          aCTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTg                                                         <  1:1873118/151‑1 (MQ=255)
                                                                                          aCTGGTGGCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTg                                                         <  8:1333238/151‑1 (MQ=255)
                                                                                             ggtggCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCt                                                       >  5:516025/1‑150 (MQ=255)
                                                                                             ggtggCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTa                                                      <  2:1694865/151‑1 (MQ=255)
                                                                                              gtggCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTAc                                                     <  6:1907121/151‑1 (MQ=255)
                                                                                              gtggCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTAc                                                     <  2:1239746/151‑1 (MQ=255)
                                                                                                ggCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACt                                                    <  6:934898/150‑1 (MQ=255)
                                                                                                ggCGCTGCTGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTg                                                   <  6:1738798/151‑1 (MQ=255)
                                                                                                   gctgctGGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGATCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGAt                                                <  7:383144/151‑1 (MQ=255)
                                                                                                      gctgGTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTGCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATa                                               <  6:1353852/149‑1 (MQ=255)
                                                                                                         ggTCTGTCTGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATaaaa                                            <  7:119400/149‑1 (MQ=255)
                                                                                                         ggTCTGTCTGATCCCAACCACTATTGGCGCCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCACCGGACCTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAAcc                                          >  5:41807/1‑151 (MQ=255)
                                                                                                          gtctgtctGATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCg                                         >  7:1568874/1‑151 (MQ=255)
                                                                                                          gtctgtctGATCCCAACCACTATTGGCGGCCTGTTGGCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCg                                         <  1:1611203/151‑1 (MQ=255)
                                                                                                            ctgtctgATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGc                                       >  8:1889371/1‑151 (MQ=255)
                                                                                                              gtctgATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCgg                                        <  2:857438/148‑1 (MQ=255)
                                                                                                                ctgATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAACACCGGCAcc                                    >  3:1284255/1‑150 (MQ=255)
                                                                                                                ctgATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCAcc                                    >  6:827626/1‑150 (MQ=255)
                                                                                                                  gATCCCAACCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCAt                                  <  5:1005871/150‑1 (MQ=255)
                                                                                                                      ccAACCACTATTGGCGGCCTGTTGTGAGCGTTCGGCGTCGCCGGGATGAGCCTGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCacac                             <  2:436579/151‑1 (MQ=255)
                                                                                                                        aaCCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTc                           <  7:813541/151‑1 (MQ=255)
                                                                                                                         aCCACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCg                          <  4:1626578/151‑1 (MQ=255)
                                                                                                                          ccACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCgg                         >  2:1863882/1‑151 (MQ=255)
                                                                                                                           cACTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGt                        <  3:454612/151‑1 (MQ=255)
                                                                                                                             cTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTa                       >  1:888325/1‑150 (MQ=255)
                                                                                                                             cTATTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATAGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGCCACCATCACACTCGGTa                       >  8:121723/1‑150 (MQ=255)
                                                                                                                               aTTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAAcc                    <  3:599974/151‑1 (MQ=255)
                                                                                                                               aTTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAAcc                    <  7:608208/151‑1 (MQ=255)
                                                                                                                               aTTGGCGGCCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAAc                     <  3:665562/150‑1 (MQ=255)
                                                                                                                               aTTGGCGGCCTGTTGTCAGCGATCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAAcc                    <  6:1226648/151‑1 (MQ=255)
                                                                                                                                  ggcggcCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGt                  <  2:752078/150‑1 (MQ=255)
                                                                                                                                     ggcCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAgg              >  2:1350313/1‑151 (MQ=255)
                                                                                                                                      gcCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAGGc             >  2:1543139/1‑151 (MQ=255)
                                                                                                                                      gcCTGTTGTCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAGGc             <  5:1335199/151‑1 (MQ=255)
                                                                                                                                         tgttgtCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAGGCGTc          >  3:270889/1‑151 (MQ=255)
                                                                                                                                          gttgtCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAGGCGTc          <  7:951487/150‑1 (MQ=255)
                                                                                                                                          gttgtCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAGGCGTCg         <  8:29682/151‑1 (MQ=255)
                                                                                                                                           ttgtCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAGGCGTCgg        >  1:689323/1‑151 (MQ=255)
                                                                                                                                             gtCAGCGTTCGGCGTCGCCGGGATGAGCCGGATGCTAGGCGCGAATGTGATTGCCACCAGCGGACGTGCAGTTGAAGCGGCAGGTGACGTTGACGTTCTGCTACTGGATAAAACCGGCACCATCACACTCGGTAACCGTCAGGCGTCGGAg      <  1:1723922/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: