Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,652,812 A→T 6.0% I163F (ATC→TTC)  ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,652,8120AT 6.0% 1310.5 / 52.4 356I163F (ATC→TTC) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (182/151);  new base T (15/7);  total (198/158)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 366 total aligned reads displayed.

ACCCAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGA  >  CP009273/1652666‑1652960
                                                                                                                                                  |                                                                                                                                                    
 cccAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCg                                                                                                                                                   <  1:3122542/149‑1 (MQ=255)
  ccAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGt                                                                                                                                                  <  1:4021409/149‑1 (MQ=255)
     aTGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTgg                                                                                                                                               >  1:2298603/1‑149 (MQ=255)
         aaCGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTggcggc                                                                                                                                           >  1:1228607/1‑149 (MQ=255)
           cgcgTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCaa                                                                                                                                         <  1:1096049/149‑1 (MQ=255)
           cgcgTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCaa                                                                                                                                         <  1:1012199/149‑1 (MQ=255)
             cgTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAata                                                                                                                                       <  1:1615995/149‑1 (MQ=255)
              gTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAatat                                                                                                                                      <  2:404293/149‑1 (MQ=255)
              gTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAatat                                                                                                                                      <  1:2040024/149‑1 (MQ=255)
                    aaaCGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCg                                                                                                                                >  1:3835361/1‑149 (MQ=255)
                      aCGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGct                                                                                                                              <  2:1951774/149‑1 (MQ=255)
                      aCGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGct                                                                                                                              >  1:2457023/1‑149 (MQ=255)
                            cgAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGc                                                                                                                        <  2:4728124/149‑1 (MQ=255)
                                ggCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATc                                                                                                                    <  1:36118/149‑1 (MQ=255)
                                   aaaTTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGc                                                                                                                 >  2:1543369/1‑149 (MQ=255)
                                    aaTTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGcc                                                                                                                <  2:1304302/149‑1 (MQ=255)
                                     aTTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCt                                                                                                               <  2:3286957/149‑1 (MQ=255)
                                     aTTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCt                                                                                                               <  2:3129307/149‑1 (MQ=255)
                                     aTTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCt                                                                                                               >  2:300372/1‑149 (MQ=255)
                                         gAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcggc                                                                                                           <  2:2089230/149‑1 (MQ=255)
                                           aCGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcggcgg                                                                                                         <  2:902074/149‑1 (MQ=255)
                                           aCGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcggcgg                                                                                                         >  1:1700120/1‑149 (MQ=255)
                                           aCGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcggcgg                                                                                                         <  2:3101878/149‑1 (MQ=255)
                                               aTTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGTTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCaa                                                                                                     >  1:3992369/1‑149 (MQ=255)
                                               aTTAGCCGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCaa                                                                                                     >  1:156433/1‑149 (MQ=255)
                                                ttAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGGGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCaaa                                                                                                    >  1:2415031/1‑149 (MQ=255)
                                                ttAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCaaa                                                                                                    <  2:433421/149‑1 (MQ=255)
                                                ttAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCaaa                                                                                                    <  1:4987042/149‑1 (MQ=255)
                                                   gCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGc                                                                                                 >  1:3276452/1‑149 (MQ=255)
                                                        gATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtgat                                                                                            >  2:1690533/1‑149 (MQ=255)
                                                          tGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtgatga                                                                                          >  1:3398407/1‑149 (MQ=255)
                                                          tGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtgatga                                                                                          <  1:4582621/149‑1 (MQ=255)
                                                           gAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGaa                                                                                         >  2:1792364/1‑149 (MQ=255)
                                                           gAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGaa                                                                                         <  2:1198195/149‑1 (MQ=255)
                                                           gAAGCCCGGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGaa                                                                                         <  1:3492660/149‑1 (MQ=255)
                                                             aGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACt                                                                                       >  2:2316614/1‑149 (MQ=255)
                                                                ccTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATAGTCAGTACTCTTCGGGGTTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCt                                                                                    >  2:2737286/1‑149 (MQ=255)
                                                                    gATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcgg                                                                                >  1:3584097/1‑149 (MQ=255)
                                                                    gATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcgg                                                                                >  1:987240/1‑149 (MQ=255)
                                                                          aTCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTc                                                                          >  1:877463/1‑149 (MQ=255)
                                                                            cGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCAc                                                                        <  1:4407848/149‑1 (MQ=255)
                                                                              ccAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTc                                                                      >  1:1461750/1‑149 (MQ=255)
                                                                                   aGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCa                                                                 <  2:191172/149‑1 (MQ=255)
                                                                                      cTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTa                                                              >  1:3363594/1‑149 (MQ=255)
                                                                                       tGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTAt                                                             >  2:2118753/1‑149 (MQ=255)
                                                                                       tGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTAt                                                             <  2:727627/149‑1 (MQ=255)
                                                                                        gaagaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATg                                                            <  1:4241885/149‑1 (MQ=255)
                                                                                          agaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGc                                                          >  1:5141081/1‑149 (MQ=255)
                                                                                          agaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGc                                                          >  1:3166214/1‑149 (MQ=255)
                                                                                             aaaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTcc                                                       >  2:2361440/1‑149 (MQ=255)
                                                                                             aaaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTcc                                                       <  1:1837064/149‑1 (MQ=255)
                                                                                               aCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTa                                                     >  1:948919/1‑149 (MQ=255)
                                                                                               aCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTa                                                     <  1:1141083/149‑1 (MQ=255)
                                                                                                   tCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGc                                                 >  1:2603994/1‑149 (MQ=255)
                                                                                                     gAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCAc                                               >  2:1065257/1‑149 (MQ=255)
                                                                                                     gAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCAc                                               <  1:2874724/149‑1 (MQ=255)
                                                                                                      aaCAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACt                                              >  2:4941675/1‑149 (MQ=255)
                                                                                                         aaTATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgcg                                           >  2:4152408/1‑149 (MQ=255)
                                                                                                             tGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCaaa                                       >  2:1291357/1‑149 (MQ=255)
                                                                                                               gCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAtt                                     >  1:2143038/1‑149 (MQ=255)
                                                                                                                cAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAttt                                    <  1:4508854/149‑1 (MQ=255)
                                                                                                                  aTGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTcc                                  <  1:635546/149‑1 (MQ=255)
                                                                                                                  aTGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTcc                                  >  2:3090341/1‑149 (MQ=255)
                                                                                                                  aTGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTcc                                  <  1:3903905/149‑1 (MQ=255)
                                                                                                                   tGAGGCTGTATATATTCAGTACTCTTCGGGGTTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCt                                 >  1:2925412/1‑149 (MQ=255)
                                                                                                                     aGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtgt                               <  2:4664786/149‑1 (MQ=255)
                                                                                                                     aGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtgt                               >  2:3913514/1‑149 (MQ=255)
                                                                                                                       gCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGc                             >  2:2940308/1‑149 (MQ=255)
                                                                                                                       gCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGc                             >  2:4875907/1‑149 (MQ=255)
                                                                                                                       gCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGc                             <  2:2659100/149‑1 (MQ=255)
                                                                                                                       gCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGc                             >  2:4977081/1‑149 (MQ=255)
                                                                                                                        cTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGcc                            <  2:2851693/149‑1 (MQ=255)
                                                                                                                        cTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGcc                            >  1:2306567/1‑149 (MQ=255)
                                                                                                                        cTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGcc                            >  2:1535341/1‑149 (MQ=255)
                                                                                                                         tGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCa                           <  1:4126286/149‑1 (MQ=255)
                                                                                                                         tGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCa                           <  1:4953968/149‑1 (MQ=255)
                                                                                                                         tGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCa                           <  1:4297095/149‑1 (MQ=255)
                                                                                                                           tatataTTCAGTACTCTTCGGGGTTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCAtt                         >  2:4003966/1‑148 (MQ=255)
                                                                                                                            atatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGc                        <  2:2210973/149‑1 (MQ=255)
                                                                                                                               tatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCt                     <  1:4764637/149‑1 (MQ=255)
                                                                                                                               tatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCt                     <  2:3847642/149‑1 (MQ=255)
                                                                                                                                atTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTa                    <  2:3584097/149‑1 (MQ=255)
                                                                                                                                   cAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTacac                 >  2:5084454/1‑149 (MQ=255)
                                                                                                                                    aGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACAcc                >  2:4450823/1‑149 (MQ=255)
                                                                                                                                    aGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACAcc                <  1:1542104/149‑1 (MQ=255)
                                                                                                                                      tACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTa              >  2:2085765/1‑149 (MQ=255)
                                                                                                                                      tACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTa              >  2:1725342/1‑149 (MQ=255)
                                                                                                                                       aCTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTAc             <  2:3398407/149‑1 (MQ=255)
                                                                                                                                        ctctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACg            <  2:948919/149‑1 (MQ=255)
                                                                                                                                          ctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGc          >  2:4215449/1‑149 (MQ=255)
                                                                                                                                           ttCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCa         >  1:57805/1‑149 (MQ=255)
                                                                                                                                           ttCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCa         <  1:2576682/149‑1 (MQ=255)
                                                                                                                                               nggATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTaa     >  1:3694554/2‑149 (MQ=255)
                                                                                                                                               gggATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTaa     <  2:2052791/149‑1 (MQ=255)
                                                                                                                                                 gATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAatg   >  2:1861497/1‑149 (MQ=255)
                                                                                                                                                 gATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAatg   >  2:1426432/1‑149 (MQ=255)
                                                                                                                                                  aTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAatga  >  2:2529018/1‑149 (MQ=255)
                                                                                                                                                  aTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAatga  >  2:519003/1‑149 (MQ=255)
                                                                                                                                                  aTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAatga  <  1:99203/149‑1 (MQ=255)
                                                                                                                                                  aTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAatga  <  1:1424955/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ACCCAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGA  >  CP009273/1652666‑1652960

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