Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP012306472,8680TG25.6% 552.7 / 13.1 203D435A (GAC→GCC) cusSsensory histidine kinase in two‑component regulatory system with CusR, senses copper ions
Reads supporting (aligned to +/- strand):  ref base T (29/122);  new base G (52/0);  total (81/122)
Fisher's exact test for biased strand distribution p-value = 8.70e-28
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 272 total aligned reads displayed.

TGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGTTC  >  AP012306/472719‑473017
                                                                                                                                                      |                                                                                                                                                     
   tGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCaa                                                                                                                                                    <  2:3394965/150‑1 (MQ=255)
   tGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCaa                                                                                                                                                    >  1:3157796/1‑150 (MQ=255)
    gAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAAc                                                                                                                                                   >  1:3694738/1‑150 (MQ=255)
     aaGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAAcg                                                                                                                                                  >  1:2478332/1‑150 (MQ=255)
      aGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAAcgc                                                                                                                                                 <  1:3819007/150‑1 (MQ=255)
          ttAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATa                                                                                                                                             <  1:4105969/150‑1 (MQ=255)
          ttAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATa                                                                                                                                             <  1:698534/150‑1 (MQ=255)
             aGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGaa                                                                                                                                          >  1:3442638/1‑150 (MQ=255)
             aGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGaa                                                                                                                                          >  2:2563527/1‑150 (MQ=255)
              gCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGaaa                                                                                                                                         <  1:3580491/150‑1 (MQ=255)
              gCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGaaa                                                                                                                                         >  2:1494259/1‑150 (MQ=255)
                    aaTGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAAAGGGCa                                                                                                                                   >  2:1431371/1‑150 (MQ=255)
                     aTGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCaa                                                                                                                                  >  2:1177982/1‑150 (MQ=255)
                     aTGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACCGTCaa                                                                                                                                  >  2:3123992/1‑150 (MQ=255)
                        tgATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACa                                                                                                                               >  1:3733851/1‑150 (MQ=255)
                        tgATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACa                                                                                                                               >  1:2697997/1‑150 (MQ=255)
                         gATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACaa                                                                                                                              <  1:3366768/150‑1 (MQ=255)
                            aaCAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATcg                                                                                                                           <  1:1019892/150‑1 (MQ=255)
                                  ccTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTaa                                                                                                                     <  2:969028/150‑1 (MQ=255)
                                   cTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGGTAGAAACGGGCAAACAAACGCGGGaaa                                                                                                                    >  2:845813/1‑149 (MQ=255)
                                      gTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGc                                                                                                                 <  1:3310799/150‑1 (MQ=255)
                                        cccccGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTTCCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTc                                                                                                               <  2:3132952/150‑1 (MQ=255)
                                         ccccGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCg                                                                                                              <  1:875546/150‑1 (MQ=255)
                                         ccccGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCg                                                                                                              <  1:1241569/150‑1 (MQ=255)
                                         ccccGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCg                                                                                                              <  2:2150756/150‑1 (MQ=255)
                                         ccccACGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCg                                                                                                              >  1:2626738/1‑150 (MQ=255)
                                           ccGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCggg                                                                                                            >  1:589993/1‑150 (MQ=255)
                                            cgcgCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGgc                                                                                                           <  1:1542326/150‑1 (MQ=255)
                                              cgcATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGgcgc                                                                                                         <  2:2694162/150‑1 (MQ=255)
                                                  tCCGACGTTACCGCAGCCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATg                                                                                                     <  1:1920678/150‑1 (MQ=255)
                                                  tCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATg                                                                                                     <  1:1494259/150‑1 (MQ=255)
                                                   ccGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATgg                                                                                                    <  1:2072368/150‑1 (MQ=255)
                                                     gACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGc                                                                                                  <  1:3258107/150‑1 (MQ=255)
                                                     gACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGGGGTCGGGGGCAATGGGc                                                                                                  >  2:552538/1‑150 (MQ=255)
                                                      aCGTTACCGCAACCGTCCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGAGGTGTAACACCAGAGAAAAAATCAAAAAATCTCAGTCAGAGGTCGGTCGCGATGGGCg                                                                                                 >  2:1254780/1‑150 (MQ=255)
                                                       cGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGt                                                                                                <  1:3984285/150‑1 (MQ=255)
                                                            ccGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCgg                                                                                           <  1:1179617/150‑1 (MQ=255)
                                                            ccGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGGACCgg                                                                                           >  1:203710/1‑150 (MQ=255)
                                                                aaCCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCgggggt                                                                                       >  1:2119793/1‑148 (MQ=255)
                                                                aaCCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTACCGGGGtt                                                                                       >  1:4024925/1‑150 (MQ=255)
                                                                 aCCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACCGGCAAACAATCCCGGTGAGAGCTCGGGCGCAAAGGGGGTACCGGGGgtt                                                                                      >  2:2097626/1‑150 (MQ=255)
                                                                 aCCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAAAGGGCAAAAAATCGCGGGAAGTGGCCCGGCGCAATGGGGGGACCGGGGgtt                                                                                      >  2:1974650/1‑150 (MQ=255)
                                                                   cGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTc                                                                                    <  2:1645074/150‑1 (MQ=255)
                                                                      gCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGTCAAAAAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTtgac                                                                                 >  2:2649822/1‑150 (MQ=255)
                                                                       ccTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTcgacg                                                                                <  2:926992/150‑1 (MQ=255)
                                                                         tttATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGAt                                                                              <  1:97069/150‑1 (MQ=255)
                                                                         tttATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAAAGGGCAAACAATCTCGGGAAGGGGTCGGGGGGCCTGGAGGTAACGGGGCTTTTGACCAt                                                                              >  2:238848/1‑150 (MQ=255)
                                                                         tttATGCGAGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTGTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGAt                                                                              <  2:3111584/150‑1 (MQ=255)
                                                                          ttATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGCCAACGCGATAGAAACGGGCAAAAAAACGCCGAAAATGGTCGGGCCCCATGGGCCTTCCGGGGGTTTCGACGATa                                                                             >  2:1201912/1‑150 (MQ=255)
                                                                          ttATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATa                                                                             >  1:3610356/1‑150 (MQ=255)
                                                                            aTGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAAACGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAAc                                                                           >  2:2028928/1‑150 (MQ=255)
                                                                             tGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCCACCCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGGGCACCGGGGGTTTCGACGATAACt                                                                          >  2:728454/1‑150 (MQ=255)
                                                                              gcgcGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACtt                                                                         <  1:830421/150‑1 (MQ=255)
                                                                               cgcgACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAAACGCGGGAAGTGCTCGGGCGCAAAGGGCGTACCGGGGGTTTCGACGATAACTTg                                                                        >  2:2264427/1‑150 (MQ=255)
                                                                                gcgACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGGCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGc                                                                       >  1:1829392/1‑150 (MQ=255)
                                                                                 cgACAACGATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAAAGGTCAAACAATCGCGGTAAGTGGTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCa                                                                      >  2:433944/1‑150 (MQ=255)
                                                                                    caacaaTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCa                                                                   >  1:3129606/1‑150 (MQ=255)
                                                                                     aacaaTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAAAGGTCAAACAATCGCGGTAAATGCTCGGGCGCAATGGGCGGACCGGGGTTTTTGACGATAACTTTCACCAg                                                                  >  2:4081327/1‑150 (MQ=255)
                                                                                      acaaTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGGCAAACAATCGCGGGAAGTGGTCGGGCGCGATGGGCGTAACTGGGTTTTCGACGATAAGTTGCAcgctg                                                                 >  2:3926982/1‑146 (MQ=255)
                                                                                         aTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGTGa                                                              >  1:1359868/1‑150 (MQ=255)
                                                                                           cGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATc                                                            <  1:1331857/150‑1 (MQ=255)
                                                                                             aTTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCCACGCGATAGAAACGGTCAAACAAACGCGGGAAGTGCTCGGGCGCAAAGGGGGGACCGGGGGTTTTGACGATAAATTGCACCCAGGGATCGa                                                          >  2:1961200/1‑150 (MQ=255)
                                                                                             aTTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGGACCGGGGGTTTCGACGATAACTTGCACCAGGGGATCGa                                                          >  2:3342571/1‑150 (MQ=255)
                                                                                             aTTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGGGATCGa                                                          >  1:2914687/1‑150 (MQ=255)
                                                                                              tttCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACACCATAGAAACGGGCAAAAAATCGCGGGAAGTGCTCGGGGGCAAAGGGCGGAACGGGGTTTTGGACGGTAACTTGCACCAGGGGATCggc                                                         >  2:2917776/1‑148 (MQ=255)
                                                                                              tttCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGGTAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCCATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGGGATCGAc                                                         >  2:412005/1‑150 (MQ=255)
                                                                                               ttCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGAcc                                                        >  1:2523811/1‑150 (MQ=255)
                                                                                                   cTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCt                                                    <  1:1639931/150‑1 (MQ=255)
                                                                                                    tATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCCATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTg                                                   >  1:3900116/1‑150 (MQ=255)
                                                                                                    tATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTTGACGATAACTTGCACCAGGGGATCGACCGTCTg                                                   >  1:2576588/1‑150 (MQ=255)
                                                                                                     aTCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGGTAGAAAAGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAAAGGGGGTACCGGGGGTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGa                                                  >  2:2115089/1‑150 (MQ=255)
                                                                                                     aTCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGGCAAAAAATCGCGGTAAGGGCTCGGGCGCAATGGGCGGACCGGGGGTTTCGACGATAACTTGCACCAGGGGATCGACCGTCTgg                                                  >  2:3224069/1‑149 (MQ=255)
                                                                                                       cGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCCGTTGAAACGGTCCAACAATCGCGGTAAGTGGTCGGGCGGCATGGGGGTAACCGGGGTTTCGGCGATAACCTCCACCAGGGGATCCACCGTGTGACa                                                >  2:3919366/1‑150 (MQ=255)
                                                                                                       cGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGGCAAACAATCGCGGTAAGTGCTCGGGCGCCATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGGGATCGACCGGCCGACa                                                >  2:876035/1‑150 (MQ=255)
                                                                                                          cAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACCCGATAGAAACGGACCAACAATCGCGGGAAAAGCTCGGGGGCACTGGGTGGAACGGGGTTTTCGACGATAACTTGCACCAGGTGAACGACCATCTGTCAgcc                                             >  2:2585076/1‑149 (MQ=255)
                                                                                                           aGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGAAAGAAACGGCCCAACAATCCCCGTAAGTGCTCGGGCGCAATGGGGGGACCGGGGCTTTCGACGCTAACTTTGCCCCCGAGATCGACCGTCTGACAgcgc                                            >  2:412848/1‑150 (MQ=255)
                                                                                                            ggCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGTCAACGCGATAGAAACGGGCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCa                                           >  1:2295357/1‑150 (MQ=255)
                                                                                                            ggCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCa                                           >  1:49207/1‑150 (MQ=255)
                                                                                                              ccAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGGCAAACAATCGCGGGAAGGGCTCGGGGGCAATGGGGGGACCGGGGGTTTTGGCGATAAATTGCACCAGGGGGTCGGCCGGCTGACAGGGCCCt                                         >  2:810255/1‑150 (MQ=255)
                                                                                                                aaTACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACCGTCAAACAATCGCGGTAAGTGCTCGGGCGCGATGGGGGTACCGGGGGTTTCGACCATAACTTTCACCAGGTGATCGACCGTCTTACAGCGCACTAc                                       >  2:1894787/1‑150 (MQ=255)
                                                                                                                aaTACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCGCGGGCGCGATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTAc                                       >  2:3739573/1‑150 (MQ=255)
                                                                                                                 aTACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACa                                      <  1:310044/150‑1 (MQ=255)
                                                                                                                  tACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACaa                                     <  2:101294/150‑1 (MQ=255)
                                                                                                                    ccGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAAtt                                   <  1:1042140/150‑1 (MQ=255)
                                                                                                                     cGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGGAGTGGTCGGGGGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTg                                  >  2:1379990/1‑150 (MQ=255)
                                                                                                                     cGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCCAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTg                                  >  1:1629990/1‑150 (MQ=255)
                                                                                                                      gCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGt                                 <  1:1720562/150‑1 (MQ=255)
                                                                                                                        tACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGtct                               <  1:926867/150‑1 (MQ=255)
                                                                                                                        tACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAACGCACTACAATTGtct                               >  2:2854563/1‑150 (MQ=255)
                                                                                                                           cTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAAAAATCGCGGTAAGGGCTCGGGCGCAATGGGCGGAACGGGGGTTTCGCCGATAAATTGCACCAGGGGATCGACCGTCCGACAACGCACTACAATTGtctctc                            >  2:474717/1‑150 (MQ=255)
                                                                                                                                 cTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGCTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTg                      <  1:2978414/150‑1 (MQ=255)
                                                                                                                                      cGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCCGTTGGCGta                 >  2:2334888/1‑150 (MQ=255)
                                                                                                                                        cTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGtata               <  1:3019664/150‑1 (MQ=255)
                                                                                                                                        cTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTAGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGtata               >  2:3530739/1‑150 (MQ=255)
                                                                                                                                          ggCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAAc             >  1:1405526/1‑150 (MQ=255)
                                                                                                                                                gAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAgg        >  2:4033477/1‑150 (MQ=255)
                                                                                                                                                gAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAgg        <  1:2150526/150‑1 (MQ=255)
                                                                                                                                                 aGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGCGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAggg       <  2:2782418/150‑1 (MQ=255)
                                                                                                                                                    ggTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGt    <  1:3851154/150‑1 (MQ=255)
                                                                                                                                                     gTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGtt   >  1:1107055/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGTTC  >  AP012306/472719‑473017

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