Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,652,812 A→T 7.2% 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 7.2% 1200.8 / 57.4 343I163F (ATC→TTC) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (180/137);  new base T (15/10);  total (196/147)
Fisher's exact test for biased strand distribution p-value = 8.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 354 total aligned reads displayed.

ATTACCCAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGAT  >  CP009273/1652663‑1652961
                                                                                                                                                     |                                                                                                                                                     
aTTACCCAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGa                                                                                                                                                       >  1:338355/1‑150 (MQ=255)
aTTACCCAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGa                                                                                                                                                       <  1:5047323/150‑1 (MQ=255)
aTTACCCAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGa                                                                                                                                                       <  1:2333178/150‑1 (MQ=255)
aTTACCCAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATAGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGa                                                                                                                                                       >  2:4655068/1‑150 (MQ=255)
 ttACCCAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGAt                                                                                                                                                      >  2:1976363/1‑150 (MQ=255)
 ttACCCAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGAt                                                                                                                                                      <  1:4355844/150‑1 (MQ=255)
 ttACCCAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGAt                                                                                                                                                      >  2:178436/1‑150 (MQ=255)
 ttACCCAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGAt                                                                                                                                                      >  2:3492780/1‑150 (MQ=255)
  tACCCAATGTATCGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATc                                                                                                                                                     <  2:3593786/150‑1 (MQ=255)
  tACCCAATGAAACGCGTTGTCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGTATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATTGCAATGTGGCTGTATATATTCAGTAATCTTCGGGGATc                                                                                                                                                     <  2:1553007/150‑1 (MQ=255)
   aCCCAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCg                                                                                                                                                    <  1:856214/150‑1 (MQ=255)
    cccAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGt                                                                                                                                                   <  1:5225132/150‑1 (MQ=255)
    cccAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGt                                                                                                                                                   >  2:4806478/1‑150 (MQ=255)
       aaTGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATCTATTCAGTACTCTTCGGGGATCGTTgg                                                                                                                                                >  2:3368152/1‑150 (MQ=255)
       aaTGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTgg                                                                                                                                                >  1:1579187/1‑150 (MQ=255)
              cgcgTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAat                                                                                                                                         >  2:2094670/1‑150 (MQ=255)
                 gTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATg                                                                                                                                      <  2:6534/150‑1 (MQ=255)
                  tGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGa                                                                                                                                     >  1:4969296/1‑150 (MQ=255)
                     gCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGAccc                                                                                                                                  >  1:3869504/1‑150 (MQ=255)
                      cAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCg                                                                                                                                 >  2:4735050/1‑150 (MQ=255)
                        aaCGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGct                                                                                                                               <  2:763841/150‑1 (MQ=255)
                         aCGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGTTCGTTGGCGGCAATATGACCCGctc                                                                                                                              <  2:239382/150‑1 (MQ=255)
                               cgAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGcc                                                                                                                        <  1:3663057/150‑1 (MQ=255)
                                   ggCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCa                                                                                                                    >  2:2174780/1‑150 (MQ=255)
                                     cAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGTTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGc                                                                                                                  >  1:1459327/1‑150 (MQ=255)
                                      aaaTTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGcc                                                                                                                 >  1:1711854/1‑150 (MQ=255)
                                        aTTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTc                                                                                                               <  1:4677901/150‑1 (MQ=255)
                                         ttCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcg                                                                                                              <  2:4940911/150‑1 (MQ=255)
                                          tCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcgg                                                                                                             >  2:4763407/1‑150 (MQ=255)
                                            gAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATCTATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcggcg                                                                                                           >  2:2887917/1‑150 (MQ=255)
                                            gAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGTTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcggcg                                                                                                           >  2:113224/1‑150 (MQ=255)
                                            gAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcggcg                                                                                                           <  1:5138269/150‑1 (MQ=255)
                                            gAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACAGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcggcg                                                                                                           >  2:46396/1‑150 (MQ=255)
                                               cGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTc                                                                                                        >  1:1500627/1‑150 (MQ=255)
                                                        tGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCCGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtg                                                                                               <  1:4811924/150‑1 (MQ=255)
                                                            aTGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtgatgt                                                                                           >  1:4457208/1‑149 (MQ=255)
                                                            aTGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtgatga                                                                                           >  2:3256234/1‑150 (MQ=255)
                                                            aTGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtgatga                                                                                           <  2:1237506/150‑1 (MQ=255)
                                                              gAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAAc                                                                                         <  1:3558313/150‑1 (MQ=255)
                                                              gAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAAc                                                                                         <  2:536202/150‑1 (MQ=255)
                                                               aaGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGTTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACt                                                                                        >  1:3491123/1‑150 (MQ=255)
                                                               aaGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACt                                                                                        <  1:5018655/150‑1 (MQ=255)
                                                                  cccTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCt                                                                                     <  1:2482567/150‑1 (MQ=255)
                                                                    cTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAc                                                                                   <  1:1242919/150‑1 (MQ=255)
                                                                    cTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAc                                                                                   >  1:2925595/1‑150 (MQ=255)
                                                                    cTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATAATCAGTACTCATCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCAACAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAc                                                                                   >  1:3176100/1‑150 (MQ=255)
                                                                     tGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcg                                                                                  >  1:4596698/1‑150 (MQ=255)
                                                                     aggATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcg                                                                                  >  1:2494309/2‑150 (MQ=255)
                                                                       gATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcggc                                                                                <  1:268605/150‑1 (MQ=255)
                                                                            tATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGTTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTc                                                                           <  2:1584909/150‑1 (MQ=255)
                                                                            tATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTc                                                                           >  1:4200210/1‑150 (MQ=255)
                                                                              tCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCAc                                                                         >  1:3004321/1‑150 (MQ=255)
                                                                                gCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTCCGGCGGTTCACTc                                                                       >  2:2917190/1‑150 (MQ=255)
                                                                                      aGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAg                                                                 <  2:1885429/150‑1 (MQ=255)
                                                                                            aagaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGc                                                           >  1:4483044/1‑150 (MQ=255)
                                                                                             agaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCt                                                          >  1:4239720/1‑150 (MQ=255)
                                                                                              gaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGGCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTc                                                         >  2:1437549/1‑150 (MQ=255)
                                                                                               aaaaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTcc                                                        >  2:3330058/1‑150 (MQ=255)
                                                                                                aaaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCt                                                       >  1:3430819/1‑150 (MQ=255)
                                                                                                aaaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCt                                                       <  2:4398400/150‑1 (MQ=255)
                                                                                                 aaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTa                                                      >  1:1135833/1‑150 (MQ=255)
                                                                                                 aaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTa                                                      <  2:3047960/150‑1 (MQ=255)
                                                                                                   ccGTCGAACAATATGGCACTGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACa                                                    <  2:3237047/150‑1 (MQ=255)
                                                                                                       cGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCAc                                                >  1:946120/1‑150 (MQ=255)
                                                                                                        gAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGTTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACt                                               <  2:3538270/150‑1 (MQ=255)
                                                                                                        gAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACt                                               >  2:3752951/1‑150 (MQ=255)
                                                                                                          aCAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgc                                             >  1:2742523/1‑150 (MQ=255)
                                                                                                            aaTATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgcgc                                           >  1:1673227/1‑150 (MQ=255)
                                                                                                             atatGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCa                                          >  1:4913815/1‑150 (MQ=255)
                                                                                                                tGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAt                                       <  1:1065731/150‑1 (MQ=255)
                                                                                                                 ggCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAtt                                      >  2:4118738/1‑150 (MQ=255)
                                                                                                                  gCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAtgt                                     >  2:3666792/1‑148 (MQ=255)
                                                                                                                    aaTGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTcc                                   >  1:5159571/1‑150 (MQ=255)
                                                                                                                       gAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtgt                                >  1:216075/1‑150 (MQ=255)
                                                                                                                        aGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtgtg                               <  1:3275346/150‑1 (MQ=255)
                                                                                                                          gCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGcc                             <  2:4149168/150‑1 (MQ=255)
                                                                                                                           cTGTATATATTCAGTACTCTTCGGGGTTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCa                            <  2:385264/150‑1 (MQ=255)
                                                                                                                           cTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCa                            <  1:2831947/150‑1 (MQ=255)
                                                                                                                            tGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCAt                           <  2:455096/150‑1 (MQ=255)
                                                                                                                            tGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCAt                           >  2:1837005/1‑150 (MQ=255)
                                                                                                                             gTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATg                          >  1:2013816/1‑150 (MQ=255)
                                                                                                                             gTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATg                          <  2:2055565/150‑1 (MQ=255)
                                                                                                                             gTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATATCCAGTGCCaag                          >  2:2662077/1‑148 (MQ=255)
                                                                                                                              tatataTTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGc                         <  1:558201/150‑1 (MQ=255)
                                                                                                                                tatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGccc                       >  1:2052443/1‑150 (MQ=255)
                                                                                                                                 atatTCAGTACTCTTCGGGGTTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCt                      >  2:4429785/1‑150 (MQ=255)
                                                                                                                                  tatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTa                     >  2:5243922/1‑150 (MQ=255)
                                                                                                                                   atTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTac                    >  2:228671/1‑150 (MQ=255)
                                                                                                                                    ttCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTaca                   <  2:404639/150‑1 (MQ=255)
                                                                                                                                    ttCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTaca                   <  1:4332940/150‑1 (MQ=255)
                                                                                                                                        gTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTa               <  2:2253399/150‑1 (MQ=255)
                                                                                                                                         tACTCTTCGGGGATCGTTGGCGGCAATATGACCCTCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTAc              <  1:2642852/150‑1 (MQ=255)
                                                                                                                                         tACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTAc              >  2:2467403/1‑150 (MQ=255)
                                                                                                                                             ctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCa          >  1:1003316/1‑150 (MQ=255)
                                                                                                                                               tCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCACACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGt        <  1:156356/150‑1 (MQ=255)
                                                                                                                                               tCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGt        <  2:2036777/150‑1 (MQ=255)
                                                                                                                                                 ggggATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTaa      >  1:74798/1‑150 (MQ=255)
                                                                                                                                                 ggggATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTaa      >  1:2207230/1‑150 (MQ=255)
                                                                                                                                                    gATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAatga   <  1:2566657/150‑1 (MQ=255)
                                                                                                                                                    gATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAatga   >  2:151364/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTACCCAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGAT  >  CP009273/1652663‑1652961

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: